similar to: package kinship - %*%

Hello, I try to define a global variable. My example: R> a <- "old" R> test <- function () { a <- "new" } R> test() R> a # shoud be "new" This doesn't work. I would like to modify the variable "a" in a procedure. How can I do that. Thank you for helping. Sven Kn侟ppel (Germany-Berlin)

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, my problem is that I don't know how long must be an array of double while calling C from R. R-Code: > array <- c(1,1,1) > save <- .C ( "Ctest" , a = array ) C-Code: void Ctest ( double *array ) { ... array = (double*) realloc ( array , new_number * sizeof(double) ) ; ... } The length of "array" will be compute in C. At the end save$a has a length

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

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,

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

Dear r-help members. I would like to compare species numbers of moths between eight different forests (each sampled for six nights). I would like to do a nested anova to compare species numbers between forests and nights. For more site specific details I wanted to do a Tukey test (TukeyHSD). Unfortunately this test is not working (message: "nicht anwendbare Methode f?r 'TukeyHSD' auf

Hello, I built a package, everything was working very well. Now I included a NAMESPACE to export only some variables. (export(f, g, ...) ) Everything is working fine, except one function. There I get an error: Fehler in UseMethod("open") : keine anwendbare Methode f?r "open" In this function I use the function "open" to generate an output bar. * if (verbose) {

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:

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

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

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:

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 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, 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

Hello, I'm looking for a method to estimate narrow sense heritability of traits in a RIL population. Papers I've checked either use either SAS or SPSS or do not give any details at all. I've found some reference to using variance components in ANOVA, using the kinship or wgaim packages, but I don't have a clue as to how to do any of this. Is there any way fro a very R illiterate

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

Hello, Apologies in advance if this is a stupid question. I am running R on Ubuntu 9. R version 2.9.2 (2009-08-24) I am trying to work with plm. I think the library is installed, as I can do > library(plm) Loading required package: kinship Loading required package: survival Loading required package: splines Loading required package: nlme Loading required package: lattice [1] "kinship