# Hardy-Weinberg moth R code

obs<-c(1469,138,5)
Tot<-2*(1469+138+5)
fA<-(2*1469 + 138)/Tot
fa<-(2*5 + 138)/Tot
fa+fA 
expected<-c(fA^2,2*fA*fa,fa^2)
chisq.test(obs,p=expected)