Nothing
gkin <- ibs(ge03d2ex.clean)
gkin[upper.tri(gkin)] <- gkin[lower.tri(gkin)]
diag(gkin) <- 1
gkin <- gkin/2
#a <- ibs(ge03d2ex.clean,snps=sample(autosomal(ge03d2ex.clean),1000,replace=FALSE))
#a[upper.tri(a)] <- gkin[lower.tri(a)]
#diag(a) <- 1
plotKinship2(2*gkin, "hist")
df <- ge03d2ex.clean@phdata
df <- mutate(df,
sex = sex + 1)
mod <- solarPolygenic(height ~ sex + age, df, kinship = gkin)
summary(mod)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.