hind <- read.table("hinds.csv",sep="\t",header=TRUE)
names(hind)
hind
plot(hind$rank,hind$ratiomales,ylim=c(0,1),col="brown3")
mod <- lm(ratiomales~rank,data=hind)
summary(mod)
mod$coef
lines(c(0,1),mod$coef[1]+c(0,1)*mod$coef[2],col="darkgreen")

qqnorm(mod$residuals)
qqnorm(hind$ratiomales)

plot(mod)