similar to: analysis of multivariate normal (missing values)

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

I'm discovering R (very impressive), and didn't find in the docs a simple method for replacing, in a data frame, missing values (NA) with the column's mean (or any other method for reconstructing missing values when needed). Thanks in advance for your help. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hi, I have been trying to dyn.load a library (rq.obj), which will allow me to run a quantile regression function, but so far unsuccessfully. I have tried under windows 98 and R 6.24: 1) dyn.load("d:\\...\\rq.obj") 2) dyn.load("d:\...\rq.obj") 3) 1 and 2 accounting for case sensitivity. 4) dyn.load("d:/.../rq.obj") 5) Place the files in the directory where from I

I compiled R-1.2.2 on my SuSE 7.1, kernel 2.4.2-2 machine with libreadline.4.1 and although I can use vi commands to scroll on the shell (bash) command line, I CANNOT do it in R. Here's the relevant parts of ./configure (and make): anovo at friedman R-1.2.2 > ./configure | grep readline checking for rl_callback_read_char in -lreadline... yes checking for readline/history.h... yes

dear list members, there seems to be a problem with the prelim.norm function (package norm) as number of items in the dataset increases. the output of prelim.norm() is a list with different summary statistics, one of them is the missingness indicator matrix "r". it lists all patterns of missing data and a count of how often each pattern occured in the dataset. as the number of items and

Dear R-users, Is there a procedure to identify neighbors in a regular lattice using either a "rook" or a "queen" criterium? To be more specific, suppose: My lattice: 1 2 3 4 5 6 7 8 9 "Rook" Neighbors 1 has neighbors 2,4 2 has neighbors 1,3,5 ... 5 has neighbors 2,4,6,8 "Queen Neighbors" 1 has neighbors 2,4,5 2 has neighbors 1,3,4,5,6 For each case,

Hi, Before I re-invent the wheel, is there a function in R similar to S+'s arima.sim, i.e., a function that simulates arima processes. ts and tseries packages don't seem to have such function, but I may have overlooked it. Thank you for your time, Alvaro Novo R Version 1.1.1 SuSE 6.4 Linux KDE 2.0 Beta 5 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-

Hi, It seems that there's a bug in summary, in the max. output... but max() alone works fine. > hw04.dframe$area ... [41] 1790 1380 1296 2745 798 2306 438649 1481 1559 2450 ... > summary(hw04.dframe) area Min. : 798 1st Qu.: 1349 Median : 1690 Mean : 6962 3rd Qu.: 2306 Max. :438600 ### should read 438649 or, to the point,

Regardless of which version -- 1.1.1 or 1.2.0 (2000-11-27) -- with a fresh "directory" (i.e. no .RData), I am getting an extremely weird result. R : Copyright 2000, The R Development Core Team Version 1.2.0 Under development (unstable) (2000-11-27) > jj _ seq(0.05,0.95,by=0.002) > sum(jj==0.75) ## WRONG ANSWER [1] 0 > 0.05 + 350*.002 ## Double check that 0.75 is in jj [1]

Hi, I'm struggling to compile R-0.90.1 in a RedHat 6.1 machine. (I have the rpm, but I would like to compile it, as I do under SuSE 6.3). I have also installed quite a few other libraries (not installed by default, in particular, *-devel-* ones), but to no avail thus far. Here's the compilation error after going through ./configure [root at friedman R-0.90.1]# make <snip>

Hi, We have a user reporting the following problem > % R > > R : Copyright 2001, The R Development Core Team > Version 1.2.3 (2001-04-26) > > R is free software and comes with ABSOLUTELY NO WARRANTY. > You are welcome to redistribute it under certain conditions. > Type `license()' or `licence()' for distribution details. > > R is a collaborative project

Hi, I spotted an inconsistency between the output of funfits under R, more specifically "nnreg", and the one obtained under S-plus. Notice the difference in the "Root MSE", "GVC" and "GVC cost=2". As a check mark, I am also including the data summaries, confirming that I was in fact using the same data set. Here's my output from both systems:

Greetings, I'm getting this error message using Joe Shafer's "NORM" package. Error in storage.mode(x) <- "double" : (list) object cannot be coerced to type 'double' I'm not sure what this means. I get this message when running s <- prelim.norm(filename) prelim.norm is used to develop some summary statistics, etc. before implementing the EM

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 I always have following error message when I try to read a big ascii-file: > inzp<-read.data() > Error: heap memory (1953 Kb) exhausted [needed 0 Kb more] read.data() is a small function that reads the ascii-file. When I cut the ascii-file to a small one I don't have this problem. Can I extend this 'heap memory' for reading big data files? How? I'm working with R on a

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

Hello, I'm looking for a more descriptive manual/tutorial/paper for the pan package. The provided manual and example do not contain any useful hints how to specify a model with more than one variable and leaves several questions unanswered. This also applies to the referred paper "Schafer: Imputation of missing covariates under a multivariate linear mixed model." Can anyone

Hi, I forgot to add to my initial message that I am using Windows 98 (if it's relevant). I need to perform survival analysis and could not find the S+ functions: - survfit (fits a Kaplan-Meier survival curve), or - coxph (fits a Cox proportional hazards model) in R. Are there any alternative functions available? Thanks, Alvaro -- -- -- -- -- -- -- -- Alvaro A. Novo

Hi, I would like to confirm w/ the other members if the function persp( ) has limitations such as not being able to label the axis... or, am I missing something? ?persp seems to indicate that's not available, but... Alvaro -- -- -- -- -- -- -- -- Alvaro A. Novo 2116 S. Orchard St., #304

Hi! The postscript device in 0.90 does not produce symbols (pch=) as it should. E.g.: postscript(file="test.ps",width=8,height=6.4,horizontal=FALSE) plot(1:10,1:10,pch=1:5) legend(1,9,letters[1:5],pch=1:5) dev.off() The result looks quite different compared to the X-windows version. The symbol 1 (pch=1) is much to small, while the other symbols seem to big. Is it a bug? How can I