similar to: an issue about missing data

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

Dear list I'm using the reshape command and want to reshape a wide data set to a long one e.g. I have the variables y1,y2,y3,age1,age2,age3,sex,ethnic I want my new long data set to consist of the variables y (which has been created from y1,y2,y3), age (which has been created from age1,age2,age3), sex and ethnic I have tried to use the command:

Hi, I am an R newbie and one thing I am having trouble with binding variables that I have created within one data frame into a new data frame when using na.omit(). To illustrate this problem I will give the example I am working on and the approah I have been using:- data.frame1<-filepath.... attach(data.frame1) #create a new variable using a function new.variable<-rep(1,length(weight3))

Hi all, I am running regressions with many covariates, most of which remain the same each time (control variables). Instead of writing 30 demographic variables every regression, is there a way I could call them all at once using a variable called, perhaps "demog"? I have tried: > demog <- list(age1, age2, age3) but I get an error when I try to call a list in a regression. I also

Hi all, I am hoping someone can help me out with this: If I have dataframe of years and ages and the first column and first row are filled with leading values: Df<- age1 age2 age3 Yr1 1 0.4 0.16 Yr2 1.5 0 0 Yr3 0.9 0 0 Yr4 1 0 0 Yr5 1.2 0 0 Yr6 1.4 0 0 Yr7 0.8 0 0 Yr8 0.6 0 0 Yr9 1.1 0 0 Now the rest of the cells need to be filled according to the previous year and age

Hi, I'm using kinship2 to calculate heritabilty, but I would like calculate in GenABEL too. I trying the code: > require(kinship2) > require(GenABEL) > pedig = with(Dados, pedigree(id=IID, dadid=PAT, momid=MAT, sex=SEX, famid=FID, missid=0)) > kmat = kinship(pedig) > (mod1 = polygenic(altura ~ SEX + idade, data=Dados, kin=kmat)) Erro em intI(i, n = d[1],

Dear Sir/Madam, Since a few day now I try to use the command "polygenic" from the GenAbel package. However, I keep bumping up against an error message: "Error in polygenic(Testo, kin = kinship, data = data1) : dimension of outcome and kinship.matrix do not match". My data exists of 1240 individuals with 74 markers. It mainly consists of small families (2 or more brothers,

Hi all, The last e-mails about beginners gave me the courage to post a question; from a beginner's perspective, there are a lot of questions that I'm tempted to ask. But I'm trying to find the answers either in the documentation, either in the about 15 free books I have, either in the help archives (I often found many similar questions posted in the past). Being an (still actual)

Hi all! I am trying to install package - GenABEL in a Linux (RedHat) environment using R 2.6.2. The CRAN Mirror starts up, I select the site, but I end up with the following error: ****************************************************************** Warning: unable to access index for repository http://cran.stat.ucla.edu/src/contrib Warning messages: 1: In open.connection(con, "r") :

Dear all, since yesterday, I have been experiencing problems with the package GenABEL. When I try to load the package (library(GenABEL)) I get the following error message: Loading required package: MASS Error : .onLoad failed in loadNamespace() for 'GenABEL', details: call: stringSplit[[1]] error: subscript out of bounds Error: package/namespace load failed for ?GenABEL? The funny

Hi All; I have just successfully installed Microsoft: Age of Empires II, Age of Empires II-Conquerors expansion, Age of Empires III on Wine 1.3.7. No dll overides, no pissing about nothing, just dropped in the CD, mapped it to using winecfg and installed. Empires II installed directly from disk and would not run after the install until I applied the no_Cd_patch and not are fine. AoEIII has 3

Hi, I had originally posted regarding an error when trying to install package - GenABEL - it has now become clear that R is not able to connect to the net. Below are a couple of things I've tried with the resulting errors I am running R 2.6.2 on Linux e15 X86 How can I fix this? Thank you! 1) download.file("http://www.nytimes.com", destfile<-tempfile()) Error: trying URL

Does anyone use the R library GenABEL? I am using it to calculate SNP interactions. I have a list of 100 SNPs, I need to look at the interaction between each of two SNPs among the list. my question is how to perform this in GenABEL. I want to use the "lm" function, but don't know how to use the SNP information. for example: result <- (lm(y~SNP1+SNP2+SNP1*SNP2)) the problem here

Hello all, First I have done extensive searches on this forum and others and nothing seems to work. So I decided to post thinking someone could point me to the write post or give me some help. I have drawn a 100 samples from a fictitious population (N=1000), and then randomly selected 25% of the 100 samples. I would like to now merge the data.frame from the 100 samples with the data.frame for

On Thu, Jan 31, 2013 at 10:29 AM, Weijia Wang <wwang.nyu@gmail.com> wrote: > Hi, > > > > I have a new question about subsetting in R. > > > > Say we have this data frame: > > > > PT_ID Blood_Pressure OBS_TYPE > > 92 1900 90.0 DBP > > 94 1900 90.0 DBP > > 174 2900 140.0 SBP > > 176 2900

Hi all! I am using GenABEL on R for GWAS analysis. I am having a couple of issues: First, I am having a problem reading files (.map, & .ped, size 900Mb, using windows 32-bit) onto R in the "convert.snp.ped" statement. I am thinking this problem is likely due to the large size of the files & my version of R is not able to handle them, since I can read in smaller files.

Greeting R helpers J I am not familiar with R but I have to use it to analyze data set that I have (30,000 20,000) I want to change the structure of the dataset and I am wondering how that might be possible in R A main data looks like this: some entities are empty Age No. Age No. Age No. Center1 5 2 8 7

Dear R help, I have the following data frame: structure(list(prochi = c("ind_1", "ind_1", "ind_1", "ind_1", "ind_1", "ind_1", "ind_1", "ind_1", "ind_1", "ind_1"), date_1st_event = structure(c(14784, 14784, 14784, 14784, 14784, 14784, 14784, 14784, 14784, 14784 ), class = "Date"),

Hi all! I am desperately trying to figure out the solution to this error, but nothing as of yet is working. As noted in an earlier post I am using GenABEL. In an attempt to read in the phenotype file, in the format .dat, R keeps giving me the error "more columns than column names" I have tried to read in the data without the headers; I have also tried to trim the data to remove any

Dear colleagues, I had a question wrt the car package. How do I evaluate whether a simpler multivariate regression model is adequate? For instance, I do the following: ami <- read.table(file = "http://www.public.iastate.edu/~maitra/stat501/datasets/amitriptyline.dat", col.names=c("TCAD", "drug", "gender", "antidepressant","PR",