In MoisesExpositoAlonso/popgensim: Wright Fisher simulations: C++ implementation with interactive Shiny app visualization

knitr::opts_chunk$set(echo = TRUE)

n <- 100
x <- rnorm(n)
y <- 2*x + rnorm(n)
out <- lm(y ~ x)
library(xtable)
tab <- xtable(summary(out)$coef, digits=c(0, 2, 2, 1, 2))
print(tab, type="html")

MSE= 1/N * sum( (Y - (X %*% b))^2 )
var_b= rep(MSE,M) %*% solve(t(X) %*% X)

MSEgwa= 1/N * sum( (Y - (X %*% bgwa))^2 )
var_bgwa= rep(MSEgwa,M) %*% solve(t(X) %*% X)

out=matrix(c(MSEgwa,MSE)
rownames(out)<-"MSE"
colnames(out)<-c("bgwa","b")

tab <- xtable(out),digits=c(2,2))
print(tab, type="html")

MoisesExpositoAlonso/popgensim documentation built on May 7, 2019, 8:57 p.m.