search for: prelim

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

Hey eveyone! I hope someone can help wiht this question. I have a matirux of all zeros and ones and I would like to indentify all unique patterns in the rows andthe number of times the pattern occurs. I changed all zeros to NA tried to use prelim.norm to identify all patterns of missing data in the rows. I got the message Warning message: NAs introduced by coercion Any ideas of how to get this to work? Or are there any way to indentify all the unique patterns in a huge matrix? ( 10000 x 71) Thanks for any suggestions!!...

'sodexquan' is an data frame of numeric variables with some missing values. I would like to approximate these missing values using some algorithm. I thereby chose the mix package. But when trying to run the 'prelim.mix' function (aka the function that will prepare the data for further calculations), R says: prelim.mix(truc,0) Error in as.integer.default(list(c(32, 52, 32, 27, 34, 35, 35, 28, 42, : the object (list) can't be converted in 'integer' I roughly translated the word in english, b...

...in his little green book on missing data (Sage 2002). The dataset for which I am having problems, "usnews", can be found at: http://www.ats.ucla.edu/stat/books/md/default.htm. I am working on a Windows machine with R 2.5 installed, all packages up-to-date. The problem has to do with the prelim.norm() function of the package "norm". Specifically, I need to use this pre-processing function to later use the EM algorithm and DA procedures in the norm package. I am getting an error with the following code. ---------------------- > pre <- prelim.norm(usnews) Error in as.dou...

Have you used Tinn-R and what landmines await the inexperienced? I could not understand why a script that used to work stopped working. Look at these two scenarios I opened an excel spreadsheet and copied several cells to the clipboard Then Scenario 1 Executed from Tinn-R > prelim<-read.delim("clipboard") > str(prelim) 'data.frame': 0 obs. of 1 variable: $ prelim..read.delim.clipboard.: logi Scenario 2 Executed from R editor > prelim<-read.delim("clipboard") > str(prelim) 'data.frame': 18 obs. of 13 variables: -- F...

Hi everyone, I'm trying to use prelim.norm with a "big" matrix (36 columns by 10000 or so rows). I found that prelim.norm has a built-in limit of 30 columns, but I'd still really like to use it for my data. Does anyone know of a different way to do the same thing? Or, would it be easier to try to modify the source code...

hello, We've got a problem with the mix package in order to impute the missing data. After importing the data, the prelim function does not work (only the stlouis data works). We have done : >library(mix) >Manq <- read.table("C:/.../file.txt") >attach(Manq) >save(Manq,file="C:/../R/rw2001/library/mix/data/Manq.rda") >data(Manq) >Manq V1 V2 V3 V4 1 1 1 1 1 2 1 1 3 6 3...

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 algorithm for missing data. Please note that this is a data fusion exercise and there is considerable amounts of missing data by design. Thanks in advance, David -- =====================...

The results, for anyone interested, can be found here http://www.psc.edu/networking/projects/hpn-ssh/results.html Long story short, there doesn't seem to be any notable difference between the two anymore. I still have a few more test combinations to run (<.5ms rtts and against cygwin) so the conclusions might change but, at this point, I'm no longer seeing any sort of performance

Hello all. I am inexperienced with R and am clumsily trying to work through it for specific multiple imputations Id like to run for my thesis.In running the MIX package, I keep getting an error message regarding the use of the prelim.mix command. Error in as.integer.default(list(alcohol = c(1, 1, 1, 1, 1, 1, 1, 1, 1, : (list) object cannot be coerced to integer I cannot find the source of this error. I am assuming its somewhere in the format of my data but cant quite seem to figure it out. I have followed the pro...

...w.anicca-vijja.de/lg/dframe.Rdata")) # load object "dframe" dim(dframe) apply(dframe,2,function(x) sum(is.na(x))) # check how many NAs in the dataset #dframe <- subset(dframe,select=-c(alter,grpzugeh,is1,is4,is6,klassenstufe,mmit,vorai,vorap,voras,vorkf,vorsg,vorvb)) s1 <- prelim.norm(dframe) s1$nmis # re-check of NAs should be identical to above s2 <- prelim.norm(dframe[,1:32])# see below -> still NAs are available - _not_ imputed thetahat1 <- em.norm(s1) theta1 <- da.norm(s1,thetahat1,steps=20,showits=TRUE) # error:...

....2) x3<-rnorm(n,1,2) x4<-floor(rnorm(n)*3) y <-rnorm(n,1*x1+2*x2+3*x3+4*x4,2) w <-rnorm(n,3,1.2) ymis<-y ymis[floor(runif(10,1,n))]<-NA wmis<-w wmis[floor(runif(10,1,n))]<-NA dat<-as.data.frame(cbind(wmis,ymis,x1,x2,x3,x4)) s <- prelim.mix(dat,3) # do preliminary manipulations thetahat <- em.mix(s) # ML estimate for unrestricted model rngseed(1234567) # set random number generator seed newtheta <- da.mix(s,thetahat,steps=100) # data augmentation ximp <- imp.mix(s, newtheta, dat) # impute und...

On Monday 24 February 2003 17:32, Deva Seetharam wrote: > Hello > > using icecast1, it was possible to serve static music files. from my > prelim understanding, it seems like icecast2 can serve only streams. is > this correct? > No. icecast2 fully supports serving of static files. Mike --- >8 ---- List archives: http://www.xiph.org/archives/ icecast project homepage: http://www.icecast.org/ To unsubscribe from this list, send a...

...erent set of individuals (i.e. I need to have a model in place before I receive their information). I expected that all that model information would be stored in a single object. In other words, when I run: imp.mix(s, theta, x) I expected that theta completely specifies the model and s (created by prelim.mix) and x only contain info on incomplete the data set I want to impute. As best as I can tell this is not the case. For instance, I create a 'model object' (a general location parameter list) from the data set trainSet using em.mix as follows: sTrain <- prelim.mix(trainSet,nCategoric...

...have done: > Ent<--read.table("C:/.../File.txt") > attach(Ent) > Ent V1 V2 V3 V4 ... V16 V17 1 1 1 2 6 18 18 2 1 1 1 NA 14 17 3 1 1 2 1 16 14 .... 199 2 1 NA 7 19 18 200 2 1 3 2 14 17 > EntPrelim<-prelim.mix(as.matrix(Ent),9) > EntEM<-em.mix(EntPrelim,maxits=500) > rngseed(1234567) > EntDA<-da.mix(EntPrelim, EntEM, steps=100, showits=TRUE) Steps of data Augmentation: 1... Error in da.mix(EntPrelim, EntEM, steps=100; showits=TRUE): Improper posterior--empty cells If...

Here''s my comments for the preliminary ZFS Crypto review. - Dan Webrev: http://cr.opensolaris.org/~darrenm/zfs-crypto-gate/webrev/ General comments: DEA-1 - SCCS keywords need to be removed DEA-2 - Copyright updated ------------------------------------------------------------------ usr/src/lib/libcryptoutil/common/keyfile.c pk...

Thanks to Fernando Tusell and especially to Brian Ripley for their work on 'mix', leading to an apparently good package mow available on CRAN. Going through the R code for the function prelim.mix, I am wondering why the following method of calculation is used at one point: umd <- as.integer(round(exp(cumsum(log(d))))) (d is a vector containing, in effect, the numbers of levels of the factors in col1, col2, ... of the categorical variables. Therefore umd is a vector containing the...

hello, I want to do a single Imputation method on binary data set. Is it possible to use imp.cat from CAT package? I have a problem defining "theta" when data is binary. Do you know any references on the subject? Thank you, Sigalit. [[alternative HTML version deleted]]

Hi all, When using package MIX, I often get the error: NA/NaN/Inf in foreign function call (arg 1). For example, > s<-prelim.mix(Y,6) Error in prelim.mix(Y, 6) : NA/NaN/Inf in foreign function call (arg 1) or when I ran: MI<-vector("list",5) #<--vector of complete data after MI fit.model.mi<-vector("list",5) rngseed(1234567) #<-- set random number generator seed margins<-c(1,2) #<-...

...g. In my program, I first source("sump6d2.q") to get data from splus, and attach(sump6d2), Whenever I run this part of program, the error pops up: library(mix) Y<-cbind(race2,nodes,htypeed1,htypeed2,ngraded2,blood2,agem40,agem40sq,e r,pr) #<--Y has to be a matrix for mix s<-prelim.mix(Y,6) #----------------------------run multiple imputation------------------------------# MI<-vector("list",10) #<--vector of complete data after MI fit.model.mi<-vector("list",10) rngseed(1234567) #<-- set random number generator seed for (i in 1:10){ cat(&...