similar to: How to extract p value from the lmekin object obtained by fitting mixed model with function lmekin() in package coxme?

On 09/06/2012 05:00 AM, r-help-request at r-project.org wrote: > Hi, R experts > > I am currently using lmekin() function in coxme package to fit a > mixed effect model for family based genetic data. How can I extract the p > value from a lmekin object? When I print the object in R console, I can > see the p value and Z value are just over there. But I can not extract

Version 2.2 of coxme has been posted to CRAN, Windows versions and mirrors should appear in due course. This is a major update with three features of note: 1. A non-upwardly compatable change: Extractor functions: beta= fixed effects, b=random effects nlme lme4 coxme <2.2 coxme 2.2 lmekin 2.2 ------------------------------------------------------ beta

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

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,

Hello, Why the correlation between the random effects is negative? library(coxme) rats1 <- coxme(Surv(time, status) ~ (1|litter), rats) random.effects(rats1)[[1]] #one value for each of the 50 litters print(rats1) rats2 <- lmekin(time ~ (1|litter), rats) fixed.effects(rats2) random.effects(rats2)[[1]] #one value for each of the 50 litters print(rats2)

I've moved nlme from Depends to Imports in my coxme package. However, a few of the examples for lmekin use one of the data sets from nlme. This is on purpose, to show how the results are the same and how they differ. If I use data(nlme::ergoStool) the data is not found, data(nlme:::ergoStool) does no better. If I add importFrom(nlme, "ergoStool") the error message is that

I'm having a problem rebuilding a package, new to me in R 2.15.0 (Linux) It hits all that contain the line \usepackage[pdftex]{graphics} and leads to the following when running R CMD check on the directory. (I do this often; a final run on the tar.gz file will happen before submission.) Since I float and resize my figures, removing the line is fatal in other ways.

Hi, I had a problem with the lmekin() in kinship package: lmekin() can not be wrapped into another function library(kinship) #creat an example dataset xx<-rnorm(100) yy<-rnorm(100) id<-1:100 test.dat<-as.data.frame(cbind(xx,yy,id)) rm(xx,yy,id) a<-bdsmatrix(rep(10,10),rep(block,10),dimnames=list(c(1:100),c(1:100))) #100x100 block (n=10) diagonal matrix to indicate the

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,

Take a look at lmekin() in the coxme package. The motivating data set for my development of coxme was the Minnesota Family Breast Cancer project: 24050 subjects in 462 families. The random effect is an intercept per subject with sigma^2 K as its variance where K is the kinship matrix (1 for self-self, .5 for parent-child or sib-sib, .25 for uncle-neice, etc). lmekin is a linear models front

The current coxme code has functions that depend on bdsmatrix and others that depend on Matrix, both those pacakges define S4 methods for diag. When loaded, the message appears: replacing previous import ?diag? when loading ?Matrix? Questions: 1. Do I need to worry about this? If so, what can I do about it? I suppose I could add an importFrom directive, but it will be a pain unless there

Dear R experts; I am trying to extract the p values from a coxme object (package coxme). I can see the value in the model output, but I wanted to have the result with a higher number of decimal places. I have searched the mailing list and followed equivalent suggestions for nlme/lme objects, but I wasn't successful. Thanks; Catarina [[alternative HTML version deleted]]

Someone just sent me a data set that causes the lmekin function, part of the kinship package, to fail. In chasing it down I get an error I have never seen before. fit <- lmekin(icam1 ~ factor(center) + age + factor(sex), random= ~1|iid, data=chaidata, varlist=kmat) Error in Y - fitted : non-numeric argument to binary operator Add the recover option, and the offending lines are

Hi, I am new to R and AIC scores but what I get from coxme seems wrong. The AIC score increases as p-values decrease. Since lower AIC scores mean better models and lower p-values mean stronger effects or differences then shouldn't they change in the same direction? I found this happens with the data set rats as well as my own data. Below is the output for two models constructed with the rats

> Message: 24 > Date: Tue, 31 Jan 2006 18:22:52 +0000 > From: "Lutz Ph. Breitling" <lutz.breitling at gmail.com> > Subject: [R] Mixed-effects models / heterogeneous covariances > To: r-help at stat.math.ethz.ch > Message-ID: > <2e38a1c80601311022i2e1be92doa60b80b50b69eb0c at mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > Dear

Dear all, I am using "coxme" function in Kinship library to fit random treatment effect nested within centre. I got 3 treatments (0,1,2) and 3 centres. I used following commands, but got an error. > ugroup=paste(rep(1:3,each=3),rep(0:2,3),sep='/') > mat1=bdsmatrix(rep(c(1,1,1,1,1,1,1,1,1),3),blocksize=rep(3,3),dimnames=list(ugroup,ugroup)) >

Hello, and thanks for your time I'm trying to extract standard errors to produce confidence intervals from a multivariable coxme model object so I can write a function that will print a summary for some reproducible research. As far as I can glean, the SE is produced on-the-fly by the print method. I'll dig into the source code if I have to, but I'd rather not have to. Any help

Dear all, I have encountered the following problem where coxme seems to allow model with only random effect in R 2.11.1 but not in R 2.13.0. Following is the error message using rat example data. Any comment on this is appreciated. In R2.13 > library(coxme) > rat1 <- coxme(Surv(time, status) ~ rx + (1|litter), rats) > rat0 <- coxme(Surv(time, status) ~ (1|litter), rats)

Dear all, This probably stems from my lack of understanding of the model, but I do not understand the variance of the random effect reported in coxme. Consider the following toy example: #------------------------------- BEGINNING OF CODE ------------------------------------------------ library(survival) library(coxme) #--- Generate toy data: d <- data.frame(id = c(1:100), #

Dear List, In documents (Therneau, 2003 : On mixed-effect cox models, ...), as far as I came to know, coxme penalize the partial likelihood (Ripatti, Palmgren, 2000) where as frailtyPenal (in frailtypack package) uses the penalized the full likelihood approach (Rondeau et al, 2003). How, then, coxme and coxph(...frailty(id, dist='gauss')) differs? Just the coding algorithm, or in