similar to: Request comments on missing data diagnosis code

Greetings, I'm in the painful process of migrating from SAS to R. In the process I've discovered that there are some basic things that I am getting hung up on. The most basic is the simple recoding of variables. Suppose I create a vector x <- rnorm(10000) and I want to recode all values of x > 1.5 to NA. How would I do that in R? TIA Cheers, Patrick

Greetings again, Seems I've touched a nerve with my previous post and after a few kind folks alerted me to a better distribution method I chose to put the code here: http://egad.psych.arizona.edu/mddiag.final.R.gz Any tips, criticisms, or additions are welcomed. Thanks for your patience. -- Cheers, Patrick _______________________________________________________________ Patrick E.

________________________________ Title: Scientific Programmer ? Simulation Analysis Salary: $50-70k/y Full Description and Application Portal: https://uacareers.com/postings/31668 We seek a scientific programmer to join the agricultural data science initiative at the University of Arizona. Our team is responsible for developing automated pipelines for data collection, simulation modeling,

I'm trying to install the XML package in R 1.6.1 (>install.packages("XML")). The download is OK, but during configuration the file parser.h cannot be found (it looks in libxml/ and gnome-xml/ then gives up). What is missing from my system? -- Thanks, -susan Susan J. Miller Biotechnology Computing Facility Arizona Research Laboratories Bio West 228 University of Arizona

Dear R users, Prof. Ripley just corrected my understanding of the use of weights in glm, which I thought would allow me to correctly use lmer. However I'm still having problems. My data takes the form of # of infected and uninfected individuals that were measured over time under different treatments. I'm using lmer to adjust for the repeated measures over time. In fitting the

Hello, I am a new R user. I am trying to calculate vector correlations for all pairwise comparisons in my data frame without repeats. I am familiar with the expand.grid function, but this includes repeats. Is there a way to use expand.grid and eliminate repeats? Or is there another function that can be used to do this? Thank you. Rebecca -- Rebecca Young Graduate Student Ecology &

platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 7.1 year 2003 month 06 day 16 language R I have a Dell Inspiron 4000 laptop with Windows ME. I recently installed the 1.7 version of R. The following commands are from my

Hello, Can you get eigenvalues in addition to eigevectors using prcomp? If so how? I am unable to use princomp due to small sample sizes. Thank you in advance for your help! Rebecca Young -- Rebecca Young Graduate Student Ecology & Evolutionary Biology, Badyaev Lab University of Arizona 1041 E Lowell Tucson, AZ 85721-0088 Office: 425BSW rlyoung at email.arizona.edu (520) 621-4005

Hi, I am trying to use the Tps procedure to visualize 3 variables in XYZ space, but keep getting the error messages below: Warning messages: 1: In Krig.find.gcvmin(info, lambda.grid, gcv.grid$GCV, Krig.fgcv, : GCV search gives a minimum at the endpoints of the grid search 2: In Krig.find.gcvmin(info, lambda.grid, gcv.grid$GCV.one, Krig.fgcv.one, : GCV search gives a minimum at the endpoints

I am a novice R user, and have been in trouble to get the right mixed model syntax in microarray analyses. There are three factors: Dye(2 levels), Temperature(3 levels) and Array(3 for each Temperature with a total of 9 arrays). I want to treat array as random, and to regard array variation different between Temperatures. So the model I want to seek is: Y = Dye + Temp + Dye + Temp*Dye +

Hello, R-folks: Here is a statement I use to make a data frame: iccdata <- data.frame(i=rep(1:10,rep(2,10)),j=rep(1:2,10), x=c(0.35011,0.11989,0.13081,0.09919,0.16000,0.12000,0.00000,0.00000, 0.44023,0.32977,2.67081,2.63919,0.09050,0.03950,0.44019,0.30981,0.59000, 0.57000,4.03000,3.77000)) Then here are the data: > iccdata i j x 1 1 1 0.35011 2 1 2 0.11989 3 2 1 0.13081 4

We's like to invite all of you to attend the first annual KVM Forum 2007. Please click here http://www.qumranet.com/KVMForum2007.php to be taken to the event site where you can learn more about the KVM Forum 2007 and register online. Hosted by the KVM Developer Community, the purpose of the forum is to bring together developers, testers and other technical individuals from within the community

We's like to invite all of you to attend the first annual KVM Forum 2007. Please click here http://www.qumranet.com/KVMForum2007.php to be taken to the event site where you can learn more about the KVM Forum 2007 and register online. Hosted by the KVM Developer Community, the purpose of the forum is to bring together developers, testers and other technical individuals from within the community

Hi, I want to impute the missing values in my data set multiple times, and then combine the results (like multiple imputation, but manually) to get a mean of the parameter(s) from the multiple imputations. Does anyone know how to do this? I have the following script: y1 <- rnorm(20,0,3) y2 <- rnorm(20,3,3) y3 <- rnorm(20,3,3) y4 <- rnorm(20,6,3) y <- c(y1,y2,y3,y4) x1 <-

Does anyone know R is to be cited in a publication? I've looked everywhere on the website and cannot find this. Thanks, Heather -- Heather Maughan Department of Ecology and Evolutionary Biology Biosciences West 310 University of Arizona Tucson, AZ 85701 Phone: 520-626-5108 Fax: 520-621-9190 hmaughan at u.arizona.edu

Hi. Can anyone point me to any module in R which implements "Little's Chi Square test" for MCAR. The problem is that i have around 60 behavioural variables on a 6 point categorical scale which i need to test for MCAR and MAR. What i can make out from preliminary analysis is that moderate (0.30 to 0.60) correlations may be present in several variable pairs leading me to suspect

L.S., Does anyone know if there is an R library which implements Little's MCAR test for completely at random missing values? It is implemented in SPSS and SAS, and widely mentioned in the literature. Thanks in advance! Sander van Kuijk -- View this message in context: http://r.789695.n4.nabble.com/Little-s-MCAR-test-tp2331137p2331137.html Sent from the R help mailing list archive at

Hey Guys, Ok, I've got a 2.2.9pre3 Samba box compilied with LDAP. I had to use 2.2.9cvs cause 2.2.8a has some bug with iPlanet that doesn't allow it to compilie properly.... Anyway, I've got a 2.2.8 box setup with: security = server password server = mirage (the 2.2.9 Solaris box) so when I try to loggin to the domain that the 2.2.8 box serves I get an error saying the

I'm relatively new to R, and not particularly adept yet, but I was wondering if there was a simply way to simulate missing data that are MAR, MNAR and MCAR. I've got a good work-around for the MCAR data, but it's sort of hard to work with. Josh [[alternative HTML version deleted]]

Hello: I want to use values from the output of linear models done using permuted data to construct a random distribution. The problem I am having is the extraction of a value, say the p-value or the regression coefficient, from the summary of a linear model. When summarizing a linear model I get this: Call: lm(formula = fitness ~ mm) Residuals: Min 1Q Median 3Q Max