similar to: problem with da.mix

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

I am desperately trying to impute missing data using 'imp.mix' but always run into this yucky error message to which I cannot find the solution. It's the first time I am using mix and I'm trying really hard to understand, but there's just this one step I don't get...perhaps someone knows the answer? Thanks! Jens My code runs:

Dear R-users, I haven?t found rules for logical operators. I need to select data according the following rule: Condition A & (Condition B | Condition C) How should I write it?? Is Condition A & Condition B | Condition C correct or will R execute (Condition A & Condition B) | Condition C ? Thanks for your help. Delphine Fontaine Delphine Fontaine Statistician Statistics

Does anyone know if for clinical studies the FDA would accept statistical analyses performed with R ? Delphine Fontaine

Dear R-users, I would like to create a new data frame composed of 2 columns of another data frame. But it does not give me what I want... > casesCNST[1:10,] case X1 X2 X3 X4 expected 1 A1 0 0 0 0 E 2 A2 0 0 0 1 C 3 A3 0 0 0 2 C 4 A4 0 0 0 3 C 5 A5 0 0 0 4 C 6 A6 0 0 1 0 C 7 A7 0 0 1 1 C 8

Hello all I'm running da.norm function in R for climate data rngseed(1234567) theta1=da.norm(mydata, thetahat, steps=1000,showits=T) param1=getparam.norm(mydata,theta1) As I understand the 1000 steps represent the markov chain values. Is there a way to plot them? Something like plot(1:1000, param1$mu[]). I just can't find a way to extract them out of my theta1. Thank you, Andrey.

I made an ANOVA with repeated mesures (aov(Mesure~Distance*Genre*Correct+Error(Sujet/(Distance*Genre*Correct)), data)) and I would like to extract the pvalues. The output is: ----------------------------------------------------------- Error: Sujet Df Sum Sq Mean Sq F value Pr(>F) Residuals 21 97.082 4.623 Error: Sujet:Distance Df Sum Sq Mean Sq F value Pr(>F) Distance

Dear R-users, I would like to write Greek letters followed by subscripts in a graph (on the X-axis and in a legend). I would appreciate any help. thank you, delphine courvoisier ************************************ Delphine Courvoisier Clinical Epidemiology Division University of Geneva Hospital +4122 37 29029

I am conducting a simulation study generating multivariate normal data, deleting observations to create a data set with missing values and then using multiple imputation via da.norm in Schafer's norm package. >From da.norm, I get the following error message: "Error: NA/NaN/Inf in foreign function call (arg 2)" The frequency of the error message seems to depend on the ratio of n

Dear Sweave-users, I want to print listing using sweave. Because my tables are very big, I use the longtable option. But, is it possible to recall the first line of the table (e.g the colnames line) on each new page ? Thanks for your help. Delphine

Hi, I would like to simulate data with a binary outcome and a set of predictors that are correlated. I want to be able to fix the number of event (Y=1) vs. non-event (Y=0). Thus, I fix this and then simulate the predictors. I have 2 questions: 1. When the predictors are continuous, I can use mvrnorm(). However, if I have continuous, ordinal and binary predictors, I'm not sure how to simulate

Hello all, When I used commands "ecm.mix and dabipf.mix" to do a simulation (sample size is small 100), I got an error : Steps of ECM, missing value where True/False needed. I've checked the menu, and the option "prior" of ecm.mix said that if structural zeros appear in the table, hyperparameters for those cells should be set to NA. However, it didn't say how to do

Hello, I am trying to translate a matlab program into R, and I have some problems to find out how to translate the PPVAL function. I know that ppval(pp,xx) returns the value at the points xx of the piecewise polynomial contained in pp( constructed by spline), but I don't know how to do the same thing in R. Thanks for you help Delphine

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

Hello all, When I ran my R program in R version 2.1.1 in windows XP and 2000, and I downloaded a package call "MIX", I got an application error: RGui:Rgui.exe - Application Error: The instruction at "0x1009d8a1" referenced memory at "0xbde48f58". The memory could not be "read". Click on ok to terminate the program click on cancel to debug the program

Hi All, I am trying to creat a simple example byusing ipf function in R, but i could not get it succefully...I am very new to R, does anyone could help, to instruct me about this ipf fucntion? Actually, this is what I mean 50 | 50 ---------------------- 33.4| 28.57 | 14.29 33.3| 23.81 | 4.762 33.3| 9.523 | 19.05

Hi! I've a problem with a game, which don't want to start... I received this error Code: daniele at daniele-laptop:~$ wine .wine/drive_c/Programmi/Delphine\ Software/Darkstone/Darkstone.exe fixme:d3d:IWineD3DImpl_FillGLCaps OpenGL implementation supports 32 vertex samplers and 32 total samplers fixme:d3d:IWineD3DImpl_FillGLCaps Expected vertex samplers + MAX_TEXTURES(=8) >

Dear guRus, I would like to loop over a medium amount of Sweave code, including both R and LaTeX chunks. Is there any way to do so? As an illustration, can I create a .tex file like this using a loop within a .Rnw file, where the "1,2,3" comes from some iteration variable in R? ################################################ \documentclass{article} \usepackage{Sweave} \begin{document}

I've got a data frame having numerical data by zip code: ZIP DATA 94111 12135.545 93105 321354.65654 94111 545.555 94706 558858.66 ... ... I'm using this function to group records by ZIP and calculate the median of DATA: aggregate(d$DATA, list(Zip = d$ZIP), FUN=median, na.rm=T) but what I really want to do is to calculate several statistics (median,

Hi all, I created a predicted variable based on a model, but somehow not all observations have a predicted value. When I tired to add the predicated value to the main data set (data$pr <- pr) , it said: replacement has 34333 rows, data has 34347. Do you know how to solve that? Thanks, Manqing [[alternative HTML version deleted]]