search for: xinxi813

Has anyone got any advice about what hardware to buy to run lots of R analysis? Links to studies or other documents would be great as would be personal opinion. We are not currently certain what analysis we shall be running, but our first implementation uses the functions lme and gls from the library nlme. To do one data point currently takes 1.5 seconds on our 3 year old sunfire box, and

Hi everyone, Does anyone here has experience using MICE to impute missing value? I am having problem to pool the imputed dataset for a MANOVA test, could you give me some advice please? Here is my code: > library(mice) >

Hi all, I have a question, is there any R package dealing with latent transition analysis with both categorical and continuous indicators? So far what I found from GOOGLE are only packages dealing with latent class analysis. So what about the longitudinal situation? Any way we could look at the transition from one class to another across time points? Thank you very much. ya [[alternative

How do I remove a "$" character from a string sub() and gsub() with "$" or "\$" as pattern do not work. > sub("$","","ABC$DEF") [1] "ABC$DEF" > sub("\$","","ABC$DEF") Error: '\$' is an unrecognized escape in character string starting "\$" >

I am trying to merge two data frames, but one of the column headings are different in the two frames. How can I rjoin or rbind the tho frames? Johnny # Generate 2 blocks by confounding on abc d1 <- conf.design(c(1,1,1), p=2, block.name="blk", treatment.names = c("A","B","C")) d2 <- conf.design(c(1,1,1), p=2, block.name="blk",

Hi list, I am wondering if there is a way to use EM algorithm to handle missing data and get a completed data set in R? I usually do it in SPSS because EM in SPSS kind of "fill in" the estimated value for the missing data, and then the completed dataset can be saved and used for further analysis. But I have not found a way to get the a completed data set like this in R or SAS. With

Hi, Is it possible to read, say 100 dat files in one time? I want to combine these 100 data set into one single dat file. Any suggestions? Thank you very much. ya [[alternative HTML version deleted]]

Hi everyone, I am trying to get the factor score for each individual case from a principal component analysis, as I understand, both princomp() and prcomp() can not produce this factor score, the principal() in psych package has this option: scores=T, but after running the code, I could not figure out how to show the factor score results. Here is my code, could anyone give me some advice please?

Hello again, I'd like to determine if an Excel file is open or writable. Can anyone help me with that? I write some stats to an .xlsx Excel file using the xlsx package. I can't write to the file unless its closed. How do I determine if the .xlsx file is open or closed so I can write to it? I've looked at file.info and file.access and I couldn't get those to work for me. Any

Hi list, I am wondering if there is a way to use EM algorithm to handle missing data and get a completed data set in R? I usually do it in SPSS because EM in SPSS kind of "fill in" the estimated value for the missing data, and then the completed dataset can be saved and used for further analysis. But I have not found a way to get the a completed data set like this in R or SAS. With

I have data that looks like this: > df1 group id 1 red A 2 red B 3 red C 4 blue D 5 blue E 6 blue F I want a list of the groups containing vectors with the ids. I am avoiding subset(), as it is only recommended for interactive use. Here's what I have so far: df1 <- data.frame(group=c("red", "red", "red", "blue",

Dear list, I have a question about imputing 2 level data in MICE, could you give me some suggestions please? Thank you very much. The data set contains 35634 cases and 1007 variables, 280 of them are categorical variables, and the rest of them are continuous variables. On the second level, there are 198 units. I am trying to impute missing values for 270 categorical variables by using the

Dear list, I am trying to find some place on internet for R cloud computing.I have some questions about it. I followed this webpage: http://www.r-bloggers.com/crdata-org-to-shut-down/ trying to gain some knowledge about where I can find free cloud computing service for R. It looks like there are 2 option: crdata.org and cloudstat.org. Since the first one has been temporarily shut down. I

Dear list, Is there any one use MICE package deal with multilevel missing values here? I have a question about the 2lonly.pmm() and 2lonly.norm(), I get the following error quite often. Here is the code the error, could you give me some advice please? Am I using it in the right way? > ini=mice(bhrm,maxit=0) > pred=ini$pred > pred V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 V11 V12 V13 V14 V15

Dear list, I am trying to impute the two level data, I have a question about a warning. Could you give me some suggestions please? Thank you very much. Here is my code and output of mice package: > ini <- mice(try, maxit=0) > pred=ini$pred > pred FAC1_1 FAC2_1 FAC3_1 FAC4_1 FAC5_1 FAC6_1 FAC7_1 FAC8_1 FAC9_1 FAC10_1 ClassSize_1 ClassSize_2 ClassSize_3 intercept