similar to: Can S-Plus packages be used in R without modification?

Greetings folks, Shafer's S-plus package "norm" for multiple imputation of missing values in multivariate normal data has been most kindly and usefully ported to R by Alvaro A. Novo. Shafer's website http://www.stat.psu.edu/~jls/ lists four S-plus packages in all: NORM - multiple imputation of multivariate continuous data CAT - multiple imputation of multivariate

Dear all! I'd like to do multiple imputation of missing values with s-plus libraries that are provided by Shafer (http://www.stat.psu.edu/~jls/misoftwa.html). I wonder, whether these libraries are compatible or somehow convertible to R (because I don't have S-plus), so that I can use this functions using the R Program. I would be happy if you could tell me, -if it is possible to use

Hi Folks, I'm currently looking at missing data/imputation methods (including multiple imputation). S-Plus has a "missing data library". What similar resources are available within R? Or does one roll one's own? Best wishes to all, Ted. -------------------------------------------------------------------- E-Mail: (Ted Harding) <Ted.Harding at nessie.mcc.ac.uk>

Hi Folks, I'm using the 'norm' package (based on Shafer's NORM) on some data. In outline, (X,Y) are bivariate normal, var(X)=0.29, var(Y)=24.4, cov(X,Y)=-0.277, there are some 900 cases, and some 170 values of Y have been set "missing" (NA). The puzzle is that, repeating the multiple imputation starting from the same random seed, I get different answers from the repeats

Dear R-help, I am trying to impute missing data for the first time using R. The norm package seems to work for me, but the missing values that it returns seem odd at times -- for example it returns negative values for a variable that should only be positive. Does this matter in data analysis, and/or is there a way to limit the imputed values to be within the minimum and maximum of the actual

Hello R experts I am using Hmisc aregImpute and Impute (following example on page 105 of The Hmisc and Design Libraries). *My end goal is to have NAs physically replaced in my dataframe. I have read the help pages and example in above sited pdf file, but to no avail. Here is example of what I did. Ph, my data frame, is attached. > xt <- aregImpute (~ q5 + q22rev02 + q28a, n.impute=10,

I am looking for a book that discusses the theory of multiple imputation (and other methods of dealing with missing data) and, just as importantly, how to implement these methods in R or S-Plus. Ideally, the book would have a structure similar to Faraway (Regression), Pinheiro&Bates (Mixed Effects) and Wood (GAMs) and would be very modern (i.e. published within the last couple of years).

Dear alltogether, I experience very strange behavior of imputation of NA's with the NORM library. I use R 2.2.0, win32. The code is below and the same dataset was also tried with MICE and aregImpute() from HMISC _without_ any problem. The problem is as follows: (1) using the whole dataset results in very strange imputations - values far beyond the maximum of the respective column, >

Hi! I am a relatively new user of R and I use it to prepare my dissertation. I have come to some very usefull and specific libraries written for S-PLUS 4 and would like to use them in R. Is that possible? I just found out that one of these libraries has already been transfered to R, while 3 others have not. For the matter of beeing more exact I''m interested in the dealing with missing

Here is a problem that perhaps someone out here has an idea about. It vaguely reminds me of something I've seen before, but can't place. Can anyone help? For multiple imputation, there are simpler methods available if the patterns of missing data are 'monotone' --- if Vj is missing then all variables Vk, k>j are also missing, vs. more complex methods required when the

Dear David, you asked this question a while ago on the R mailing list and got no answer. I have the same problem and was wondering if you had found a solution Cheers Loic Loic Thibaut, PhD candidate, ARC Centre of Excellence for Coral Reef Studies, School of Marine Biology, James Cook University, Townsville, Qld, 4811, Australia. Tel + 61 747 815 735, Fax: + 61 747 251 570, email:

Hi, I have a set of covariance matrices but not the original data. I want to carry out some exploratory factor analysis. So, I am trying to construct a covariance matrix list as the input for factanal. I can construct a list which includes the cov, the centers, and the n.obs. But it doesn't work. I get an error that says "Error in sqrt(diag(cv)) : Non-numeric argument to mathematical

Hi In a dataset I have length and age for cod. The age, however, is ony given for 40-100% of the fish. What I need to do is to fill inn the NAs in a correct way, so that age has a value for each length. This is to be done for each sample seperately (there are 324 samples), meaning the NAs for sampleno 1 shall be calculated from the known values from sampleno 1. As for example length 55 cm

Hi all, I'm currently working with a dataset that has quite a few missing values and after some investigation I figured that multiple imputation is probably the best solution to handle the missing data in my case. I found several references to functions in S-Plus that perform multiple imputation (NORM, CAT, MIX, PAN). Does R have corresponding functions? I searched the archives but was not

Hello, well, I searched list-archive, cran and the references, but found nothing. Thus: Does anybody around here know anything about Dempster-Shafer Theory, Evidence Theory or Hints in R? Has anybody stumbled about a package that I overlooked or implemented something in this area? I really would like to not implement a hint-model a second time. My apologies if I missed something obvious, but I

Hello all. I am trying to run PAN, multilevel multiple imputation program, in R to impute missing data in a longitudinal dataset. I could successfully run the multiple imputation when I only imputed one variable. However, when I tried to impute a time-varying covariate as well as a response variable, I received an error message, “Error: subscript out of bounds.” Can anyone tell if my commands

Alistair-- I wrote functions to calculate Clinical Significance in Splus for the following article: McGlinchey, J. B., Atkins, D. C., & Jacobson, N. S. (2002). Clinical significance methods: Which one to use and how useful are they? Behavior Therapy, 33, 529-550. I will send the functions to you back-channel. cheers, Dave -- Dave Atkins, PhD Assistant Professor in Clinical Psychology

Hello, I have a general question about combining imputations as well as a question specific to the rms and Hmisc packages. The situation is multiple regression on a data set where multiple imputation has been used to give M imputed data sets. I know how to get the combined estimate of the covariance matrix of the estimated coefficients (average the M covariance matrices from the individual

Hi, I have a big dataset which has many missing values and want to implement Multiple imputation via Monte carlo markov chain by following J Schafer's "Analysis of incomplete multivariate data". I don't know where to begin and is looking for a sample R code that implements multiple imputation with EM, MCMC, etc.... Any help / suggestion will be greatly appreciated. David

In two weeks, Dave Thomas will be in Pasadena doing a three-day Rails workshop. I''m about to sign up for it. I figure three intense days with Dave and a few dozen other Rails developers will produce huge learning leaps for me. Then it occurred to me that it might be useful to ask the list if anyone has attended one of these workshops before and, if so, how valuable you found