search for: nancy1

Hi, I am trying to run PCA on a set of data with dimension 115*300,000. The columns represnt the snps and the row represent the individuals. so this is what i did. #load the data code<-read.table("code.txt", sep='\t', header=F, nrows=300000) # do PCA # pr<-prcomp(code, retx=T, center=T) I am getting the following error message "Error: cannot allocate vector of

Hi, I need to rbind two data frames. Each one has a header . after the rbind I would like to keep the header for each and have the two data frames separated by a line. Is this possible to do in R? For example weight_mean weight_sd.dev > F 14.33333 4.932883 > M 34.66667 10.692677 > > hight_mean hight_sd.dev > F 35.00000 7.071068 > M 34.66667 10.692677 --

Hi, I am conducting an association analysis of genotype and a phenotype such as cholesterol level as an outcome and the genotype as a regressor using multiple linear regression. There are 3 possibilities for the genotype AA, AG, GG. There are 5 people with the AA genotype, 100 with the AG genotype and 900 with the GG genotype. I coded GG genotype as 1, AG as 2 and AA as 3 and the p-value for the

Hi All, this could be a simple question but I am looking into modifying a data frame using a "condition" without the need to loop over that data, would that be possible? I have tried the following > x<-c(4,5,6,6,8) > y<-c("a","b","b","b","c") > data<-data.frame(x,y) > data x y 1 4 a 2 5 b 3 6 b 4 6 b 5 8 c if

Hi, I have used multiple linear regression on a data set and one if the regressor was significant with a p-value =0.01 I need to calculate the power for a multiple linear regression. i.e. do I have enough power to believe the above p-value? -- View this message in context: http://www.nabble.com/power--tp25776305p25776305.html Sent from the R help mailing list archive at Nabble.com.

Hi, I do have a file that has 500000 columns and 40 rows. I want to apply PCA on that data and this is what I did h1<-read.table("Ccode.txt", sep='\t', header=F) # reads the data from the file Ccode.txt h2<-prcomp(na.omit(h1),center=T) but I am getting the following error "Error in svd(x, nu = 0) : 0 extent dimensions" I appreciate if someone can help

Hi All, I have data for two groups, group with 100 points and group B with 15 points. i needed plot these two groups in one scatter plot, each group with a different color. I tried plot3d(data, col = c("red", "blue")[c(rep(1, 100), rep(2, 15))], pch=20) . this works but the graph does not look good. the other option that I tried is scatterplot3d(data, type="p",

Hi All, I have a problem with rbind. I have data that consist of weight height .. etc of 1000 patients. I would like to find the mean and the standard deviation ( for the weight , height etc) for each gender. data<-read.table("data.txt", header=T, sep='\t') fdata=NULL for (i in 1:50){ nn<-names(X)[i] m<-tapply(X[,i],data$gender,mean,na.rm=T) s<-tapply(X[,i],

Hi All, I have data with 400 individuals and the following information Grade: pass or fail coded as 1 for pass and 0 for fail Sex: male or female ( coded as 1 for male and 2 for female ) Age Teaching.method : can be 1,2,3 I want to fit a logistic regression where the outcome if (1=pass or 0 for fail) and the rest of the variables are the regressors. My question is that I am not sure