similar to: pb with package mix

Hello, We use the mix package and we have a problem with the DA function. We aren't sure, but it's maybbe a memory problem. We 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

Hi, we try to do a logistic regression with the function glm. But we notice that this function don't give the same results as the SAS proc catmod (differents estimate given). We try to change the contrast on R system with: > options(contrasts=c(unordered="contr.SAS",ordered="contr.poly")) We also try with brlr and logistf functions. Unfortunately, the estimate

Hi, I'm using glm function to do logistic regression and now I want to know if it exists a kind of R-squared with this function in order to check the model. Thank you.

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

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

Dear Professor Schaefer I am experiencing a technical difficulty with your mix package. I would appreciate it if you could help me with this problem. When I run the following code, R 2.5.1 and R 2.6.0 crashes. It's been tested on at least 2 windows machine and it is consistent. Execution code it's self was coped from the help file of imp.mix. Only thing I supplied was a fake dataset.

I have a couple EMC relabled ACBEL/Polytech UPSs which I'm unable to find anything about. What I do know is that ACBEL/Polytech is CoolerMaster, however, they don't make UPSs that I know of. These UPSs were used as part of an EMC Clarion array. Does anyone know anything about these? Thanks, Mark. -------------- next part -------------- An HTML attachment was scrubbed... URL:

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

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

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

Hello list, probably a known issue (but I haven't found a link yet): I'm using 3.0.2a with 'unix charset = ISO8859-15' and get the following then a file with ticks is created: - desired: "???????`'??" - displayed: "???????`'??" ^^^ The euro sign is created correctly in a file name. The filename will be displayed in the filesystem

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? Thanks! Josh

R-experts: I am trying to reproduce some of Paul Allison's results 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

Hi Folks, Good evening, I used the code which i print below, Yes i got the result , but i need the result like this "TruncateHtml is jus....(continued)" some_html = ''<ul><li><a href="http://whatever">TruncateHtml is just like the vanilla truncate rails helper , TruncateHtml is just like the vanilla truncate rails

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

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

Hello all, I'm having problems with setting custom flags on a dovecot imap setup, with mbox mail store. I was running v1.0rc15 on a Debian Etch box, upgraded recently to v1.0.15, but nothing changed. The problem : I read my mail using IMAP. In Thunderbird (various versions, current one being v2.0.0.19), I can set various message flags (Important, Personal, Work, ... mapped to IMAP as

hi. i'm a statistics student and studying bootstrap in R. i'm trying to draw bootstrap samples from a sample, using the following R code: > *boot = numeric(200)* *> {for (i in 1:200)* * boot [i] = (sample(data,replace=T))}* i obtain 200 samples but all of them are the same. i want to obtain different samples. what should i do? can you please help me if possible. thank you.

Hi there, I've seen a strange thing not reported yet AFAIK. We have W2K DCs with SP3 with Samba 3.0.2a everything works fine in regard to winbind, but with Samba 3.0.3 winbind produces schannel len 24 errors and 'wbinfo -t' and 'id DOMAIN\userid' doesn't work. Note that wbinfo -u|g works well and a join was successful as well. I tried to tune my krb5.conf but in the end I

Hello everybody! I need help with a curious problem referring to Word and Samba. I've installed a network with a Samba 3.0.2 PDC and fileserver and about 8 Windows Clients ranging from w2k prof (members of the domain) to XP home. The XP Home users have an account within the domain and so can "mount" the Samba share as a local drive on their computers. All of the users are members