similar to: lme with nested factor and random effect

Dear R users, Could you help my with the following problem? I want to repeat a glm analysis with 2 independent variables for all 900 variables (snps) in my data set. So, I want to check whether snp1 has a different effect on my outcome variable in patients and controls(phenotype). And repeat that for snp2 to snp900. Is there an easy way to get a summary of the data, e.g. a list of P values of all

Dear R-experts, My problem is how to handle a 10GB data file containing genotype data. The file is in a particular format (Illumina final report) and needs to be altered and merged with phenotype data for further analysis. PERL seems to be an frequently used solution for this type of work, however I am inclined to think it should be doable with R. How do I open a text-file, line by line,

Hoping someone can offer me some assistance. I'm trying to execute a script and I keep getting this error message about "Error: element 12 is empty". I'm wondering if my syntax is incorrect within legend.list. If anyone has any suggestions to sees something obvious that I am missing, I would greatly appreciate any help. Many Thanks, Patrick > # These are the symbols and

I have a path: path = "/nfs/users/nfs_n/ns9/ Phenotype Analysis/Results/Run_AmplRatio_neg BinaryAll trained without akapn+tnik.csv" I wish to replace the spaces with "\ " so that it can be read by a system call to unix. Using gsub I try: > gsub(" ","\\ ",path) [1] "/nfs/users/nfs_n/ns9/Phenotype Analysis/Results/Run_AmplRatio_neg BinaryAll

Dear all,   I have data from the following experimental design and trying to fit a mixed model with lme function according to following steps but struggling. Any help is deeply appreciated.   1) Experimental design: I have 40 plants each of which has 4 clones. Each clone planted to one of 4 blocks. Phenotypes were collected from each clone for 3 consecutive years. I have genotypes of plants. I

Hey, I have code that can check the quality of a data set we're working with (expression data), and I'm having some trouble writing code that would make the expression data we have tie to other data we want to link it to (called phenotype data). Does anyone have any advice on how I could make a single object that would do this? Other relevant info: I want to use the pdata() function,

Hello R Forum users, I was hoping someone could help me with the following problem. Consider the following "toy" dataset: Accession SNP_CRY2 SNP_FLC Phenotype 1 NA A 0.783143079 2 BQ A 0.881714811 3 BQ A 0.886619488 4 AQ B 0.416893034 5 AQ B 0.621392903 6 AS B 0.031719125 7 AS NA 0.652375037 "Accession"

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

Hi, I'm looking for an efficient code that will enable me to reshuffle data (phenotype) for certain number of individuals and creating a loop that will randomly simulate it for 10000 times *(permutation)*. I also need to find how I keep the information (p value for each SNP) gathered for all the 10000 iterations. My data set looks like this (n=500): Individual # Phenotype SNP1 SNP2

Hello, I have data in a dataframe with 139104 rows which is multiple of 96x1449. i have a phenotype file which contains the phenotype information for the 96 samples. the snp name is repeated 1449X96 samples. I haveto merge the two dataframes based on sid and sen. this is how my two dataframes look like dat<-data.frame(snpname=rep(letters[1:12],12),sid=rep(1:12,each=12),

Dear useRs, First off, sorry about the long post. Figured it's better to give context to get good answers (I hope!). Some time ago I wrote an R function that will get all pairwise interactions of variables in a data frame. This worked fine at the time, but now a colleague would like me to do this with a much larger dataset. They don't know how many variables they are going to have in the

Hi, I am currently writing code to input a few thousand files, run them through the Random Forests package, and then output corresponding results. When I use the code below: zz<-textConnection("ex.lm.out", "w") sink(zz)

Dear All, Does anybody have experience with R package pbatR (http://cran.r-project.org/web/packages/pbatR/index.html)? I am trying to use it to analyze the family-based case-control data, but the package totally doesn?t work on my computer. I contacted the authors of the package, but I haven?t heard anything from them. Following the package manual, I tried the simple example as below:

Hi: I use lm (linear model) to analyze 47 variables , 8 responses So I use loop to finish it . I want the program to show the results that P-value is less than 0.05. How can I cite the P-valus from lm result ? Ping The code: #using LM to model general fati for (j in 48:52) { for (i in 3:46){ gen.fat<-y_x[,j] gen.fat<-as.numeric(gen.fat) snp_marker<-y_x[,i] x<-colnames(y_x)

Dear R-Users, I have a question concerning the determination of breakpoints and comparison of slopes from broken-line regression models. Although this is rather a standard problem in data analysis, all information I gathered so far, did not answer my questions. I added a subsetted example of my data. Basically it is a timeseries of recorded phenotypes in three different groups of plants. You

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

hi there - I find i get a flicker on the penultimate frame of the effect in Firefox (1.04 and 1.06) i don''t get the same behaviour with IE. is this simply a limitation of the browser or am i doing something wrong ? cheers jonah -------------- next part -------------- An HTML attachment was scrubbed... URL:

Hi All, My attempts to build an R package on my Windows 7 computer using R V13.0 ir R V13.1 using R CMD build --binary <filename> have been failing at the penultimate step with the error message "running 'zip' failed" coming after the procedure has completed the MD5 sums step.The same build command on the same package completes successfully on my older Vista

I have a TDM22B (TDM400 PCI + 2xFXO + 2xFXS). One FXO is connected to the PSTN. When Asterisk places a call, it dials using DTMF. If I listen in on the line during the dialing, there is a roughly one second pause between the penultimate and the final digit -- proportionally much longer than the pause between the preceeding digits. This accounts for about 25% of the time taken to dial. Anyone

Colleagues, (Running R 2.7.0) I have a script that I want to delete as it completes execution. The penultimate line of the script (before the quit command) is: file.remove("Scriptname") The script is executed as: R --no-save < Scriptname In OS X and Linux this is successful and returns: > > file.remove("x") > [1] TRUE and the file is deleted In Windows XP,