Displaying 2 results from an estimated 2 matches for "makefamid".
2011 Apr 15
1
no solution yet, please help: extract p-value from mixed model in kinship package
...)
momid <- c(rep(0, 5), rep(2, 5), rep(4, 5), rep(6, 5), rep(8, 10),
rep(10, 10), rep(12, 10), rep(14, 10), rep(16, 10), rep(18, 10),
rep(20, 10), rep(22, 10) )
ped <- data.frame(id, dadid, momid)
# *****************************kmatrix**************************************
> cfam <- makefamid(ped$id,ped$momid, ped$dadid)
>
> kmat <- makekinship(cfam, ped$id, ped$momid, ped$dadid)
>
> #*****************************************x and y variables
*********************
> set.seed(3456)
>
> dat <- sample(c(-1,0,1), 10000, replace = TRUE)
>
> snpmat<- data....
2011 Apr 14
0
extract p-value from mixed model in kinship package
..., 10), rep(11, 10), rep(13, 10), rep(15, 10), rep(17, 10),
rep(19, 10), rep(21, 10))
momid <- c(rep(0, 5), rep(2, 5), rep(4, 5), rep(6, 5), rep(8, 10),
rep(10, 10), rep(12, 10), rep(14, 10), rep(16, 10), rep(18, 10),
rep(20, 10), rep(22, 10) )
ped <- data.frame(id, dadid, momid)
cfam <- makefamid(ped$id,ped$momid, ped$dadid)
kmat <- makekinship(cfam, ped$id, 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='' )...