print(save.em.try20, nlines=20)
print(save.em, nlines=20)
print(save.em.try20, nlines=20)
print(save.em, nlines=20)
print(save.em.try20, nlines=20)
print(save.em, nlines=20)
print(save.em.try20, nlines=20)
names(save.em)
save.em$haplotype
save.em$hap1code
save.em.try20$hap1code
save.em.try20$haplotype
library(rms)
score.gaus.adj
haplo.score
locator.haplo
locator.haplo()
locator.haplo(score.gaus.add)
score.gaus
score.gaus.add
score.gaus.min10
score.gaus.dom
fit.gaus
5/(2*nrow(glm.data))
50*50
residuals.glm
gaussian()
gaussian()$resid
names(gaussian())
gaussian()$variance
q()
n
