similar to: cox model for haplotypes

Hi all, I've been looking in R for an EM algorithm adjusted for multiple-locus haplotypes frequencies, but failed in 100%. Has anyone heard of anything of this kind in R? Thanks in advance, Marcin

A non-text attachment was scrubbed... Name: not available Type: text/enriched Size: 1819 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/20030214/7f3adccd/attachment.bin

Hi everyone, I am doing 5 years mortality predictive index score with survival analysis using a Cox proportional hazard model where I have a continous predictive variable and a right censored response which is the mortality, and the individuals were followed a maximum of 7 years. I'd like to asses the discrimination ability of survival analysis Cox model by computing a ROC curve and area

Hi, Does anybody know how to obtain the same information as the SPSS command "DISPLAY DICTIONARY". I know the read.spss function in the package foreign, but it doesn't give either the missing values defined by the user under SPSS or the format variables (for example if they are dates or not). Thanks in advance. Isaac Subirana (isubirana@imim.es) [[alternative HTML

Hi R-users, When I import a data base from SPSS with some variable with a name of more than 8 characters (which is possible in version 12 or later of SPSS software) I find that the name of these variables have been truncated to 8 characters. I used the 'read.spss' function from the 'foreign' package. Is there a function in R that solves this 'length variable names'

Dear R mailing list, I have a dataset with genotypes from trios and I would like to infer haplotypes for each mother, father and child. The package that I could find that can do this is tdthap. But when the mother is homozygous (e.g., 2/2) the haplotype is called as not possible to infer (0); I would prefer for it to call the genotype (2). From what I understand it is doing what I would like

Dear all, I'm using the package "Genetics", and I'm interested in the computation of D' statistics for Linkage Disequilibrium, for which the LD() command has been realised. Unfortunately I don't find any reference on "how" the D' is computed by the LD() function. In the package documentation it is generally referred as "MLE" estimation, but

I have a question about how to compare a GLM with a GAM model using anova function. A GLM is performed for example: model1 <-glm(formula = exitus ~ age+gender+diabetes, family = "binomial", na.action = na.exclude) A second nested model could be: model2 <-glm(formula = exitus ~ age+gender, family = "binomial", na.action = na.exclude) To compare these two GLM

I see something new and unexpected here. > update.packages() trying URL `http://cran.r-project.org/src/contrib/PACKAGES' Content type `text/plain; charset=iso-8859-1' length 163467 bytes opened URL .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... ......... downloaded

Hi, My query is regarding permutation test and reshuffling of genotype/phenotype data I have been using the haplo.stats package of R. for haplotype analysis and I would like to perform an analysis which I'm requesting your advice. I have a data set of individuals genotyped for 12 SNP and a dichotomous phenotype. At first, I have tested each of those SNP independently in order to bypass

Dear list, I am trying to write a package for simulating meioses in R. We defined a class 'haplotype' which contains the basic units of our simulation: setClass("haplotype",representation(snp = "numeric",qtl = "list", hID = "numeric",phID0 = "numeric",phID1 = "numeric"),

Hi All, is there any chance of vectorising the two ifelse() statements in the following code: for(i in gp){ new[i,1] = ifelse(srow[i]>0, new[srow[i],zippo[i]], sample(1:100, 1, prob =Y1, rep = T)) new[i,2] = ifelse(drow[i]>0, new[drow[i]>0,zappo[i]], sample(1:100, 1, prob =Y1, rep = T)) } Where I am forced to check if the value of drow and srow are >0 for each line... in

Dear R-help listers, I am now asking for helps on how to draw stacked ellipses to illustrate the shared and specific of multiple objects using R. My problem comes from my population genetics study. Now, I genotyped three species, and I get known about the amount of shared and specific haplotypes in each of the species and their combinations. I want to illustrate this result in three stack

I've just uploaded an R version of Daniel Schaid's 'haplo.score' library (from http://www.mayo.edu/statgen/) to CRAN. >From the DESCRIPTION file: [haplo.score is a] suite of routines that can be used to compute score statistics to test associations between haplotypes and a wide variety of traits, including binary, ordinal, quantitative, and Poisson. These methods assume

Dear friends - a couple of papers in PNAS (lastly:framework for making better predictions by directly estimating variables' predictivity, Lo et al PNAS 2016; 113:14277-14282) have focused interest on mapping complex traits to multiple loci spread all over the genome. I have been around on the relevant taskview(s) I hope but fail to see that the backward haplotype transmission association

We have R 2.2.1 installed on a Linux cluster that seems to have problems loading either of our shared object libraries for packages. This seems to be happening on both local and global versions of packages that we install. However, we have only noticed this problem in the past 3 months on this R installation, whereas some users had success before then. It could be that something on our system

Re-sending this from 10/7. Please help! I'm really clueless how to fix this: System: Solaris 7, Software: R_1.7.1 for unix > R CMD INSTALL -l /Rdir/library haplo.stats_1.1.0.tar.gz -works like a charm-- no syntax errors, etc But problems are when I'm loading the library within R. > library(haplo.stats, lib.loc="/people/biostat3/sinnwell/Rdir/library") > Error

Hi All, I'd like to ask for a few clarifications. I am doing some calculations over some biggish datasets. One has ~ 23000 rows, and 6 columns, the other has ~620000 rows and 6 columns. I am using these datasets to perform a simulation of of haplotype coalescence over a pedigree (the datestes themselves are pedigree information). I created a new dataset (same number of rows as the pedigree

Hi all, Hope someone can help me. I start out with a matrix called haplo with 600 rows and two columns. To start with, the elements of the matrix are character strings of 20 numbers. I then want to create separate objects, called ha1-ha600, that are the concatenated 20 + 20 numbers, as numeric. I am having no trouble doing anything except getting the index i to paste to the object ha(i) in

Hi, there, I cannot get accurate value for calculation. for example: ld<-sqrt(1*0.05*0.95*0.05*0.95) 0.05*0.95-ld=-6.938894e-18 0.05*0.95-ld==0 is False. I met this problem in my program, how can I handle it. Thanks. xj.