similar to: heritability estimation

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,

Hi R help, I wanted to set a fixed heritability for the object "Trait" left of ~ operator. Is this possible using lme, or would I need to use something like gls? This is the model line contained within a function; lme.out <- summary(lme(Trait ~ genotype, random = ~1|sire, weights=~number, na.action="na.omit")) Many thanks for any help.

Dear R users,   I am trying to use "wgaim" package for QTL analysis using mixed model approach. But i am stuck with "asreml" function while using "wgaim" package. Do i need a separate package to activate "asreml" function beside "wgaim" package ?   If so, i tried to download "Asreml-R" package (i guess this is right package

Dear R-list, maybe someone can help me with the following mixed-effects models problem, as I am unable to figure it out with the 'nlme-bible'. I would like to fit (in R, obviously) a so-called animal model (google e. g. "Heritability and genetic constraints of life-history" by Pettay et al.) to estimate the variance component that is due to genetic effects. The covariances of

I'm trying also to understand how to get the between-group variance out of a one-way ANOVA, but I'm beginning to think that in a sense, the variance does not exist. Emma said: *The model is response(i,j)= group(i)+ error(i,j)* Yes, if by group(i) you mean intercept + coefficient[i]. *we assume that group~N(0,P^2) and error~N(0,sigma^2) * Only the error is assumed to be a random

Dear Friends, I am seeking for any help on an error message in lme functions. I use mixed model to analyze a data with compound symmetric correlation structure. But I get an error message: "Error in corMatrix.corCompSymm(object) : NA/NaN/Inf in foreign function call (arg 1)". If I change the correlation structure to corAR1, then no error. I have no clue how to solve this problem.

Dear All, I am in the somewhat unfortunate position of having to reproduce the results previously obtained from (non-metric?) MDS on a "kinship" matrix using Statistica. A kinship matrix measures affinity between groups, and has its maximum values on the diagonal. Apparently, starting with a nxn kinship matrix, all it was needed to do was to feed it to Statistica flagging that the

> 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

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:

Is there a way to restrict connection to my asterisk server to users based on their IP addresses, and not just password. I have some hackers who connect to my server to make illegitimate solicitation calls to people. I had to shutdown the server for now until I find a solution. ANY HELP? Thanks. ond -------------- next part -------------- An HTML attachment was scrubbed... URL:

An attacker is scanning my Asterisk Switch to gain illegitimate access to VoIP call functionality. Using a sip scanning tool, *it* sends REGISTERs with random identities. And when it discovers one identity subscribed in my switch, it tries to authenticate with random passwords using this user name. For the moment, I have replaced this account. And also blocked the IP it has used but each time

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

Hi, I am often re-use commands in R by using the arrows to retype them and modify for the new needs. This way, I don't need to write commands again, and I can modify the command by inserting/deleting/re-writing the command. Today, suddenly, re-writing change the commands as the "insert" function does in Word. for example - if I am standing in the middle of

Hello, I was hoping someone could answer a question for me that may either be statistical or script related. I don't come from a statistics background, so I am not positive if I am using the correct nomenclature or even the correct procedure. Is it possible to model "random coefficients" in a mixed effects cox-regression using coxme from the Kinship package? For example, using

Hi.... I have setup a 4 node cluster on virtual servers on RHEL platform. Not able to get better performance statistics on glusterFS as compared to local file system. Kindly suggest a test run that can be checked to differentiate between them. Regards, Anish Kumar "Confidentiality Warning: This message and any attachments are intended only for the use of the intended recipient(s). are

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

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

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?

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 Im using the kinship package to draw a pedigree. On my data set this works fine but when i add indivudals to the pedigree i keep getting an error i hope someone can help me! This is the code im using: Data<-read.table("Tree.txt", header=T, sep=",") attach(Data) ped<-pedigree(id, dadid, momid, sex, aff) par(xpd=T) plot.pedigree(ped) This is my data looks like