Displaying 20 results from an estimated 2000 matches similar to: "Random Coefficients using coxme"
2006 Jun 29
2
help with coxme
Hi there,
I have a question on fitting data by coxme. In particular I want to fit a
random intercept and random slope cox model. Using the rats dataset as an
example, I generated another covariate x2 and want to specify a random slope
for x2. Here is my code:
x2=matrix(rep(runif(50), 3), 50, 3)
x2=as.vector(t(x2))
rats2=cbind(rats, x2)
But when I used the coxme function as follows, it gave
2008 Mar 05
1
coxme - fitting random treatment effect nested within centre
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))
>
2008 Dec 28
1
Random coefficients model with a covariate: coxme function
Dear R users:
I'm new to R and am trying to fit a mixed model
Cox regression model with coxme function.
I have one two-level factor (treat) and one
covariate (covar) and 32 different groups
(centers). I'd like to fit a random coefficients model, with treat and covar
as fixed factors and a random intercept, random
treat effect and random covar slope per center.
I haver a couple of
2010 Apr 23
2
Deferred Default Marker
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:
2011 Jul 08
1
coxme for random effects only model
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)
2009 Nov 13
2
error checks
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
2009 Aug 05
1
inheriting C symbols
The package coxme depends heavily on bdsmatrix, to the point of needing
access to some of its C calls. The kinship package (in progress) uses
the R level functions in bdsmatrix, but not the C calls. That is, I
don't want to always export the symbols.
For testing I can use an explicit dyn.load('somepath/bsdmatrix.so',
local=F). How do I incorporate such a dependency in the
2012 Feb 03
1
coxme with frailty--variance of random effect?
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), #
2012 Feb 19
1
coxme: model simplification using LR-test?
Hi
I'm encountering some problems with coxme
My data:
I'm looking at the survival of animals in an experiment with 3 treatments,
which came from 4 different populations, two of which were infected with a
parasite and two of which were not. I'm interested if infected animals
differe from uninfected ones across treatments.
Factor 1: treatment (3 levels)
Factor 2: infection state
2011 Jan 25
1
coxme and random factors
Hi
I would really appreciate some help with my code for coxme...
My data set
I'm interested in survival of animals after an experiment with 4
treatments, which was performed on males and females. I also have two
random factors:
Response variable: survival (death)
Factor 1: treatment (4 levels)
Factor 2: sex (male / female)
Random effects 1: person nested within day (2 people did
2011 Aug 02
2
Extract p value from coxme object
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]]
2012 Jun 30
1
How do I extract coefficient standard errors /CI for a "coxme" model
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
2007 Apr 20
1
Approaches of Frailty estimation: coxme vs coxph(...frailty(id, dist='gauss'))
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
2011 Feb 04
1
GWAF package: lme.batch.imputed(): object 'kmat' not found
Hello, All,
GWAF 1.2
R.Version() is below.
system(lme.batch.imputed(
phenfile = 'phenfile.csv',
genfile = 'CARe_imputed_release.0.fhsR.gz',
pedfile='pedfile.csv',
phen='phen1',
covar=c('covar1','covar2'),
kinmat='imputed_fhs.kinship.RData',
outfile='imputed.FHS.IBC.GWAF.LME.output.0.txt'
))
Gives the error messages:
Error in
2011 Dec 30
0
New version of coxme / lmekin
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
2007 Dec 05
4
coxme frailty model standard errors?
Hello,
I am running R 2.6.1 on windows xp
I am trying to fit a cox proportional hazard model with a shared
Gaussian frailty term using coxme
My model is specified as:
nofit1<-coxme(Surv(Age,cen1new)~ Sex+bo2+bo3,random=~1|isl,data=mydat)
With x1-x3 being dummy variables, and isl being the community level
variable with 4 levels.
Does anyone know if there is a way to get the standard error
2012 Oct 07
1
variances of random effects in coxme
Dear R users,
I'm using the function coxme of the package coxme in order to build Cox
models with complex random effects. Unfortunately, I sometimes get
surprising estimations of the variances of the random effects.
I ran models with different fixed covariates but always with the same 3
random effects defined by the argument
varlist=coxmeMlist(list(mat1,mat2,mat3), rescale = F, pdcheck = F,
2011 Apr 15
1
no solution yet, please help: extract p-value from mixed model in kinship package
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
2007 Apr 20
1
Hiding "Warning messages" in coxme output
Dear list,
I have been trying to use coxme in R 2.3.1.
When I use coxme in the following data sim.fr1, i get
"Warning messages: using 'as.environment(NULL)' is
deprecated"
Why does it occur?
How can I hide such warning message,
especially when coxme is under a loop?
Mohammad Ehsanul Karim (Institute of Statistical
Research and Training, University of Dhaka)
>
2010 Aug 22
2
coxme AIC score and p-value mismatch??
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