similar to: MiceR

Hi, I have had a quick search to see whether this question has been asked/answered before but haven't found anything directly related to it. Basically, I am wondering if I can run the packages, developed by Shafer for S-Plus, that allow multiple imputation of missing data - NORM, CAT, MIX, and PAN. If not, does anyone know if someone has done the modification that would make these packages

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

Hi, I have been trying to untar the source file for R and having no success. I have used decompression utilities that handle .tar and .tgz file but I keep getting an error that there is a "Header CRC Error". So, I am wondering if this is just me and I am not using the right decompression utility or whether there is an error in the headers for the zip file. Anyone have any advice?

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, 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, Many thanks to those of you who responded to my last post about Schafer's MI packages. I am really pleased to have access to them through R which, I have to say, is an amazing piece of software. I am only sorry that I haven't found it until now. But, to my question. Does anyone know if there is a package developed for evaluating clinical significance using Jacobson and Truax's,

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:

I am attempting to plot frequency of occurrence against (calendar) time. As the output is to be produced as a Powerpoint picture on a dark background I wish to make the plot output yellow > par(col="yellow", col.axis="yellow", col.lab="yellow", col.main="yellow") > plot(t2, t1, ylim=c(0,40), main="Episynlestes run 2",

Hi, how can I extract the t-values after running a pairwise.t.test? The output just list the p-values. Many thanks for your help. Cheers Guido ____________________________________ Guido J. Parra School of Tropical Environment Studies and Geography James Cook University Townsville Queensland 4811 Phone: 61 7 47815824 Fax: 61 7 47814020 Mobile: 0437630843 e-mail:

Hi all, I am wondering if anyone has written (or knows of) a function that will conduct a goodness-of-fit test for a gamma distribution. I am especially interested in test statistics have some asymptotic parametric distribution that is independent of sample size or values of fitted parameters (e.g., a chi-squared distribution with some fixed df), because I want to fit gamma distributions to

I am working with transition matrices of sequences of animal behaviours. What I would like to do is parse the original matrices, adjusting row/column order so that the matrix has its main values in blocks surrounding the diagonal. This would cause behaviours involved in functional groupings (e.g. grooming, resting, foraging etc) to appear as blocks. This can be done manually by applying

Hi I am experiencing a problem with the random number generator. When I call any function that involve RNG such as "runif" or "sample" I get this error: .Random.seed is not an integer vector but of type 'double' I can't coerce the random seed and it's recommended not to alter it anyway. Can anybody help me? Thanks and kind regards, Piero Visconti, PhD

I'm trying to use the biological female and male signs in R2.5.1 under Windows XP. I can access and insert these symbols using word-processors. In general these should be available as \u2640 and \u2642 but I can't make them happen in R using (say) text(5,5, "\u2640") message "invalid \uxxxx sequence" I've replaced Arial with Arial unicode MS in Rdevga and

R2WinBUGS users, I'm getting the error message "monitor could not be set" when I try to monitor (and output to R) a response variable that includes missing data (NAs)...ie.e imputing for y. Setting monitors for other parameters works OK and running the R-generated WinBUGS script file in WinBUGS, also works OK. Below is my R code and the log.txt file output by R2WinBUGS. As I want to

Hello, I recently started using the metafor package (version 1.6-0) in R (2.15.1, 64-bit Windows 7) and noticed that I was getting slightly different values when I manually calculated the standardized mean difference versus what escalc was giving me. Here''s a very simple example: escalc(measure="SMD", m1i=5,m2i=10,n1i=5,n2i=5,sd1i=1,sd2i=2,vtype="LS") The result

Not an R question directly, but has anyone got a method for placing a moderately large number of (near) equi-spaced points on a sphere? I have a nasty feeling platonic solids are needed for exact solutions and I'm thinking of samplings involving around 200 - 1000 regularly-spaced points, Thanks, Richard Rowe Richard Rowe Senior Lecturer Department of Zoology and Tropical Ecology, James

I think this was asked recently but I can't find the item. Is it possible to find the colour at a point on a graphics page? I want to estimate the area of a complex shape by colouring it then random sampling from an enclosing box ... Richard Rowe Richard Rowe Senior Lecturer Department of Zoology and Tropical Ecology, James Cook University Townsville, Queensland 4811, Australia fax (61)7

Dear people, I'm running RedHat 9.0 and R : Version 1.7.1 (2003-06-16) from the help file # Usage: # # filter(x, filter, method = c("convolution", "recursive"), # sides = 2, circular = FALSE, init) # init: for recursive filters only. Specifies the initial values of # the time series just prior to the start value, in reverse # time

I am working with some survival data with missing values. I am using the mice package to do multiple imputation. I have found code in this thread which handles pooling of the MI results: https://stat.ethz.ch/pipermail/r-help/2007-May/132180.html Now I would like to plot a survival curve using the pooled results. Here is a reproducible example: require(survival) require(mice) set.seed(2) dt

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 <-