similar to: Inquiry About R Packages for Specific Research Areas

Dear All, I wonder if someone can point me in the right direction here. I'm working with the rms library, R 2.9.2 under Windows XP. I'm trying to arrange two plots side by side for a colleague. mfrow or mfcol do not seem to work, however, so I am obviously missing something important. I know that there have been changes in the graphics from Design to rms, but am just not sure where to

I had added a vignette to the coxme package and all worked well locally, but it failed at CRAN.???? The issue is that the vignette involves using coxme for pedigree data, it doesn't work without the kinship2 package, and I hadn't put in the necessary "if (require(...." logic. The question is, how do I make the entire vignette conditional? If the package isn't available,

Dear LLVM Users, The European LLVM User Group Meeting 2011 schedule is on-line at http://llvm.org/devmtg/2011-09-16/. Thank you every one that registered and submitted papers and posters, There are 5 talks planned, several side-sessions (one of them reserved for lightning talks) and more than 40 people already registered! The current proposed list of side sessions are: - OpenCL (and similar

Hi All: I am working with a dataset on Arsenic toxicity, and I am trying to calculate the 20th, 40th, 60th, 80th, and highest percentiles for a variable, dietary Moisture (variable name dMoist). The inbuilt function quantile(dMoist) would print 0, 25th, 50th, 75th, and 100th percentile. Does there exist a function that can calculate xth percentile (where x = 10th, 20th, ... etc) values?

I would like to fit a hierarchical regression model from Witte et al. (1994; see reference below). It's a logistic regression of a health outcome on quntities of food intake; the linear predictor has the form, X*beta + W*gamma, where X is a matrix of consumption of 82 foods (i.e., the rows of X represent people in the study, the columns represent different foods, and X_ij is the amount of

Dear useRs: Announcing the release of kinship2, version 1.2.1, to CRAN. kinship2 is a branch from the original kinship package with some key updates for the pedigree and kinship functions, with some additional functions that work with the pedigree object. Highlights: * pedigree and pedigreeList objects implemented with S3 methods print and "[," * pedigree plotting routine

Dear useRs: Announcing the release of kinship2, version 1.2.1, to CRAN. kinship2 is a branch from the original kinship package with some key updates for the pedigree and kinship functions, with some additional functions that work with the pedigree object. Highlights: * pedigree and pedigreeList objects implemented with S3 methods print and "[," * pedigree plotting routine

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

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

I have a data frame containing the Id, Mother, Father and Sex from about 10,000 animals in our colony. I am interested in graphing simple family trees for a given subject or small number of subjects. The basic idea is: start with data frame from entire colony and list of index animals. I need to identify all immediate relatives of these index animals and plot the pedigree for them. We're

Hello! I am building a package, which includes also one Fortran subroutine, which works fine if I compile it as a shared library and load it into R via dyn.load(). However, when I launch R CMD check it doesn't stop with checking examples. It's just doing and doing ... I pasted the whole output from R CMD check. Does anyone have any suggestions? I'm still using R 2.0.1.

Hola a todos, A partir de los siguientes datos: d <- list(`1` = structure(list(ped = c(1L, 1L, 1L, 1L, 1L, 1L, 1L), id = 1:7, father = c(2L, 0L, 0L, 2L, 2L, 2L, 2L), mother = c(3L, 0L, 0L, 3L, 3L, 3L, 3L), sex = c(2L, 1L, 2L, 2L, 2L, 1L, 2L), affected = c(1L, 2L, 1L, 1L, 2L, 2L, 2L)), row.names = c("1", "2", "3", "4", "5",

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

Dear R developers: I am a PHD candidate student in the school of public health of Peking University and my major is genetic epidemiology. I am learning the FAM-MDR algorithm, which is used to detect the gene-gene and gene-environment interactions in the data of pedigree. The codes were written by Tom Cattaert of the University of Liege. The algorithms and the sample datasets are available at

Hi, I have a dataframe that contains pedigree information; that is individual, sire and dam identities as separate columns. It also has date of birth. These identifiers are not numeric, or not sequential. Obviously, an identifier can appear in one or two columns, depending on whether it was a parent or not. These should be consistent. Not all identifiers appear in the individual column - it is

Hola a todos, me gustaría realizar una consulta asociada a la generación de valores genéticos del pedigreemm en R. Primero generé el archivo de pedigree incluyendo los parentales para posteriormente estimar la varianza aditiva y los valores genéticos para cada individuo, relacionando los individuos por medio de la matriz de parentesco. Me da todo perfecto, el complemento pedigreemm trabaja muy

Hi, You should look at Pinheiro and Bates (2000) Mixed-effects models in S and S-Plus. It describes how to format the correlation matrix to pass to functions lme and gls. Basically, the correlation matrix has to be one of the corStruct classes, probably corSymm for your example. So in the call to lme (or gls if you really have no random effects), use something like:

I am using dgc.genetics to perform TDT analysis on SNP data from a cohort of trios. I now have a file with about 6008 variables. The first few variables related to the pedigree data such as the pedigree ID the person ID etc. Thereafter each variable is a specific locus or marker. The variables are named by a pattern such as "Genotype.nnnnn" with nnnnn corresponding to a number which

Dear all, I'd like to introduce the paramlink package, offering parametric linkage analysis in R, as well as other likelihood-based pedigree analyses. While a rudimentary version of the package has existed for a while, it was never properly introduced to this forum. The present version is a major update with new and improved functions. Features: * Singlepoint LOD scores for

Dear all, I'd like to introduce the paramlink package, offering parametric linkage analysis in R, as well as other likelihood-based pedigree analyses. While a rudimentary version of the package has existed for a while, it was never properly introduced to this forum. The present version is a major update with new and improved functions. Features: * Singlepoint LOD scores for