similar to: kinship package: drawing pedigree error

I am making the question clear. Please help. > Dear R experts > > I was using kinship package to fit mixed model with kinship matrix. > The package looks like lme4, but I could find a way to extract p-value > out of it. I need to extract is as I need to analyse large number of > variables (> 10000). > > Please help me: > > require(kinship) > > #Generating

I've been looking at the kinship package which looks as though it might be appropriate for my purposes. What I can't find is any reference to the data that is used in the example code. A dataframe called d10 with column names, upn, dadid, momid, sex and affect is required. One can get an idea of what sort of values should be in most columns from the description in the pedigree function,

I have a data frame containing the Id, Mother, Father and Sex from about 10,000 animals in our colony. I am interested in graphing simple family trees for a given subject or small number of subjects. The basic idea is: start with data frame from entire colony and list of index animals. I need to identify all immediate relatives of these index animals and plot the pedigree for them. We're

Dear R experts I was using kinship package to fit mixed model with kinship matrix. The package looks like lme4, but I could find a way to extract p-value out of it. I need to extract is as I need to analyse large number of variables (> 10000). Please help me: require(kinship) Generating random example data id <- 1:100 dadid <- c(rep(0, 5), rep(1, 5), rep(3, 5), rep(5, 5), rep(7,

Hello, I am using the Kinship package to draw a pedigree plot. I am trying to control the spacing between individuals in the pedigree. Currently, it is drawing a pedigree with too little space, so the text is overlapping. I'd like to increase the distance between individuals. According to the documentation, I would expect the following to be relevant:

Hi, I'm using kinship2 to calculate heritabilty, but I would like calculate in GenABEL too. I trying the code: > require(kinship2) > require(GenABEL) > pedig = with(Dados, pedigree(id=IID, dadid=PAT, momid=MAT, sex=SEX, famid=FID, missid=0)) > kmat = kinship(pedig) > (mod1 = polygenic(altura ~ SEX + idade, data=Dados, kin=kmat)) Erro em intI(i, n = d[1],

I wish to make modifications to the plot.pedigree function in the kinship package. My attempts to contact the maintainer have been unsuccessful, but my question is general, so specifics of the kinship package might not be an issue. My first attempt was to make a new function Plot.pedigree in the .GlobalEnv which mostly achieved what I wanted to. However, I'm sure that's not the tidiest

Hello! I am building a package, which includes also one Fortran subroutine, which works fine if I compile it as a shared library and load it into R via dyn.load(). However, when I launch R CMD check it doesn't stop with checking examples. It's just doing and doing ... I pasted the whole output from R CMD check. Does anyone have any suggestions? I'm still using R 2.0.1.

I'm currently packaging up some of the kinship matrix routines more formally, these are used in coxme when dealing with family correlation structures. One of my test programs exercises error conditions, i.e., it purposely feeds particular types of invalid pedigree data in to see if the right error message results. So there are comment-action pairs # the next line should generate a

I have library in development with a function that works when called from the top level, but fails under R CMD check. The paricular line of failure is rsum <- rowSums(kmat>0) where kmat is a dsCMatrix object. I'm currently stumped and looking for some ideas. I've created a stripped down library "ktest" that has only 3 functions: pedigree.R to create a pedigree or

Hi,I use the 'lmekin' model of the 'kinship' package of R in order to estimate heritability. I want to estimate the confidence interval of the variance coefficient and so I should use a bootstrap simulation. The pedigree file has 1386 subjects so I create a kinship matrix [1386*1386].This is the code of R I use: kfit2 <- lmekin(IT~1+AGE +(1|ID), dati1,

I've finally narrowed down a puzzling problem: here is the short test case. tmt34% R --vanilla R version 2.10.0 (2009-10-26) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 > temp <- matrix(runif(50), ncol=2) > t(temp) %*% temp [,1] [,2] [1,] 7.916016 6.049698 [2,] 6.049698 7.650694 > library(kinship) Loading required package:

Dear useRs: Announcing the release of kinship2, version 1.2.1, to CRAN. kinship2 is a branch from the original kinship package with some key updates for the pedigree and kinship functions, with some additional functions that work with the pedigree object. Highlights: * pedigree and pedigreeList objects implemented with S3 methods print and "[," * pedigree plotting routine

Dear useRs: Announcing the release of kinship2, version 1.2.1, to CRAN. kinship2 is a branch from the original kinship package with some key updates for the pedigree and kinship functions, with some additional functions that work with the pedigree object. Highlights: * pedigree and pedigreeList objects implemented with S3 methods print and "[," * pedigree plotting routine

Hola a todos, A partir de los siguientes datos: d <- list(`1` = structure(list(ped = c(1L, 1L, 1L, 1L, 1L, 1L, 1L), id = 1:7, father = c(2L, 0L, 0L, 2L, 2L, 2L, 2L), mother = c(3L, 0L, 0L, 3L, 3L, 3L, 3L), sex = c(2L, 1L, 2L, 2L, 2L, 1L, 2L), affected = c(1L, 2L, 1L, 1L, 2L, 2L, 2L)), row.names = c("1", "2", "3", "4", "5",

Hello, I have a pedigree file such this: FAMID ID FA ID MO ID SEX STATUS 1 1 0 0 2 0 1 2 3 1 2 2 1 4 3 1 1 2 1 5 3 1 1 0 1 6 3 1 1 0 1 7 3 1 2 2 1 8 3 1 1 0 1 9 3 1 1 0 1 10 3 1 2 0 1 11 3 1 1 0 1 3 0 0 1 0 2 12 13

Hi, You should look at Pinheiro and Bates (2000) Mixed-effects models in S and S-Plus. It describes how to format the correlation matrix to pass to functions lme and gls. Basically, the correlation matrix has to be one of the corStruct classes, probably corSymm for your example. So in the call to lme (or gls if you really have no random effects), use something like:

Dear R developers: I am a PHD candidate student in the school of public health of Peking University and my major is genetic epidemiology. I am learning the FAM-MDR algorithm, which is used to detect the gene-gene and gene-environment interactions in the data of pedigree. The codes were written by Tom Cattaert of the University of Liege. The algorithms and the sample datasets are available at

Hi, I have a dataframe that contains pedigree information; that is individual, sire and dam identities as separate columns. It also has date of birth. These identifiers are not numeric, or not sequential. Obviously, an identifier can appear in one or two columns, depending on whether it was a parent or not. These should be consistent. Not all identifiers appear in the individual column - it is

Hi All, I'd like to ask for a few clarifications. I am doing some calculations over some biggish datasets. One has ~ 23000 rows, and 6 columns, the other has ~620000 rows and 6 columns. I am using these datasets to perform a simulation of of haplotype coalescence over a pedigree (the datestes themselves are pedigree information). I created a new dataset (same number of rows as the pedigree