Displaying 20 results from an estimated 2000 matches similar to: "lmekin() function in kinship package"
2012 Nov 24
1
Bootstrap lmekin model
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,
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
2011 Apr 14
0
extract p-value from mixed model in kinship package
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,
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      
2010 Apr 06
0
Strange error
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
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))
>
2013 Jun 04
0
Mixed effects model with a phylogenetic tree/ distance, matrix as a random effect
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
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
2012 Sep 06
1
How to extract p value from the lmekin object obtained by fitting mixed model with function lmekin() in package coxme?
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 them
by the coef() function. kinfit$coefficient$fixed (kinfit is the name of the
lmekin object)
2020 Jan 13
0
Error in R CMD check --as-cran ?
I've been fighting a CMD check error for coxme this morning.?? I thought I had it fixed, 
but instead I had forgotton --as-cran on my last test run.? So the version just submitted 
to CRAN has warning messages in the log.
I think it is an issue with CRAN.?? I've sent a message to R-devel asking for help.?? 
Since then, as a guess, I renamed my internal "ismat" function to
2012 Sep 06
0
p value from lmekin()
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
2011 Oct 06
1
multiple defines of diag
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
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
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:
2008 Sep 15
1
'plm'/'kinship' package on Debian Etch?
Hi Folks,
Has anyone got a compiled binary of the package 'plm'
and/or of the package 'kinship' for Debian Etch?
I'm asking because I'd like to install 'plm', but have
failed because it depends on 'kinship', and 'kinship'
refuses to compile on my Debian Etch.
The underlying reason is the Debian refuses to allow me
to install the development
2011 Jun 01
0
Anyone have experience with kinship pedigree plot?
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:
2020 Jan 13
0
as-cran issue
On 13/01/2020 11:02 a.m., Therneau, Terry M., Ph.D. via R-devel wrote:
> Where can I find out (and replicate) what options as-cran turns on?
> 
> The issue: the following lines generate an error in R CMD check --as-cran? for coxme.? But
> there is no error without as-cran nor is there one when I run the code in a terminal window.
> 
> ismat <- function(x)? inherits(x,
2005 Nov 29
1
package kinship - %*%
Hello,
I like to use the package "kinship" (R version 2.2.0). After loading
this package the operator %*% doesn't work.
Example:
R> library(kinship)
R> a <- cbind(1:2,rnorm(2))
R> a%*%a
Error message:
Fehler in a %*% a : keine anwendbare Methode f侟r "%*%"
The Message in English: Error in a %*% a: no applicable message for
"%*%".
What can I do?
2004 Oct 18
1
installing package 'kinship'
Dear All,
I have problem installing the 'kinship' package. I used the
'install.packages' from R command line and had no problem installing other
packages before.
Here is what I get:
* Installing *source* package 'kinship' ...
** libs
gcc -no-cpp-precomp -I/Library/Frameworks/R.framework/Resources/include
-I/usr/local/include   -fno-common  -g -O2 -c agfit6b.c -o agfit6b.o
2005 Nov 06
1
kinship package example data
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,