Displaying 2 results from an estimated 2 matches for "covtrait".
Did you mean:
contraint
2011 Apr 15
1
no solution yet, please help: extract p-value from mixed model in kinship package
...nd y variables
*********************
> set.seed(3456)
>
> dat <- sample(c(-1,0,1), 10000, replace = TRUE)
>
> snpmat<- data.frame(matrix(dat, ncol = 100))
>
> names(snpmat) <- c(paste ("VR",1:100, sep='' ))
>
> yvar <- rnorm(100, 30, 5)
> covtrait <- rnorm(100, 10, 5)
>
> mydata <- data.frame(id, yvar, covtrait, snpmat)
>
#******************************mixed model in lmekin
*******************************************
>
> fmod <- lmekin(yvar ~ mydata[,3] , data= mydata, random = ~1|id,
> varlist=list(kmat)) $coeff...
2011 Apr 14
0
extract p-value from mixed model in kinship package
...d, ped$momid, ped$dadid)
set.seed(3456)
SNPdata <- c(1:1000)
dat <- sample(c(-1,0,1), 10000, replace = TRUE)
snpmat<- data.frame(matrix(dat, ncol = 100))
names(snpmat) <- c(paste ("VR",1:100, sep='' ))
yvar <- rnorm(100, 30, 5)
mydata <- data.frame(trait1, covtrait, snpmat)
fmod <- lmekin(y~ mydata[,i] , data= snpdata, random = ~1|id,
varlist=list(kmat)) $coefficients[2,4] # does not work
Ultimately I want to put into the loop:
for(i in 3:length(snpdata)) {
P <- vector (mode="numeric", length = 1000)
P[i] <- lmekin(trait1~ snpdata[,i]...