similar to: R-Help : running MIX package

Hi, I have a dataset with 6 categorical variables. I have used this following code to make the variables u1-u6 ordered factors and this works well. cat1cat2 cat3 cat4 cat5 cat6 ? 0 ? ?? 1 ? ? 1????? 0 ??? 0? ?? 1 ? 1 ? ?? 1 ? ? 0 ? ?? 0 ? ? 0 ? ? 0 ....... .... ############ data<-read,table("example.txt") data <- as.data.frame(lapply(data, ordered)) ############ Now,

Hello, useR:, Suppose I have two plots made by using contour() function, say Cont1 and Cont2 respectively. They have slightly difference because of the two slightly different data I used. I want to see the difference between them so I want to plot Cont2 on Cont1, are there any methods to plot it without filling the frame of Cont1 totally of Cont2. I mean, how I can integreate the two plots

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

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

All, I have been reading Dr. Harrell's excellent "Regression Modeling Strategies" book and trying out the exercises. I understand that contrast( ) is used to obtain contrasts between two variables for given levels of other nuisance variables; is there a way to use contrast( ) to obtain, for example, Scheffe confidence intervals / hypothesis tests for many post hoc contrasts at

On Fri, May 15, 2015 at 5:27 PM, Kaylor, Andrew <andrew.kaylor at intel.com> wrote: > I like the way this sorts out with regard to funclet code generation. > It feels very natural for Windows EH, though obviously not as natural for > non-Windows targets and I think it is likely to block some optimizations > that are currently possible with those targets. > Right, it will

It appears to be legitimate to include multi-level categorical and continuous variables in defining the model for earth (e.g. y ~ cat + cont1 + cont2) but is it also then possible use categoricals in the predict method using the earth result? I tried but it returns an error which is not very informative. Thanks Chris

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.

Dear alltogether, I experience very strange behavior of imputation of NA's with the NORM library. I use R 2.2.0, win32. The code is below and the same dataset was also tried with MICE and aregImpute() from HMISC _without_ any problem. The problem is as follows: (1) using the whole dataset results in very strange imputations - values far beyond the maximum of the respective column, >

I am trying to figure out how to write a conditional statement that will execute only if a specified variable is not empty. Right now my view looks like this: <dt>Building Regulations</dt> <% for link in @cat1 %> <dd><a href="<%= link.url %>"><%= link.title %></a></dd> <% end %> Action: @cat1 = Link.find(:all, :conditions

Hello, I can get the median for each factor, but I'd like another column to go with each factor. The nm column is a long name for the lvls column. So unique work except for the order can get messed up. Example: x =

This is driving me nuts:-) I''ve brought it down to a simplified case where I''m rendering this rjs template: page.insert_html :top, ''cat1'', "<li>Some list item</li>" According to FireBug, this is what I''m actually getting in the browser: try { new Insertion.Top("cat1", "<li>Some list

Hello, I'm looking for a solution for the following problem: 1) I have a folder with several csv files; each contains a set of measurement values 2) The measurements of each file belong to a position in a two dimensional matrix (lets say "B02.csv" belongs to position 2,2 3) The size of the matrix is fix 4) I cannot assure to have a csv file for each position 5) Each position

Hi, I am having some difficulties with matrix operations. It is a little hard to explain it so please bear with me. I have a very large data set, large enough that it needs to be split in parts in order to deal with. I can work things on these "parts" but the problem lies in adding together these parts for the final answer. So that been said, let's say that i split the data in 2

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

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

mydata <- read.table(textConnection(" Id cat1 location item_values p-values sequence a111 1 3002737 100 0.01 1 a112 1 3017821 102 0.05 2 a113 2 3027730 103 0.02 3 a114 2 3036220 104 0.04 4 a115 1 3053984 105 0.03 5 a118 1 3090500 106 0.02 8 a119 1 3103304

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

Hi there, I am confused about fucntion call. After defining a function, I called it within another function. dt<-cars; #a copy of R internal dataset "cars" created; dt$cat1<-ifelse(dt$speed<20,0,1); dt$ind<-ifelse(dt$speed<15,1,2); #group variable; freqtot <- function(data,var){ attach(data,warn.conflicts=FALSE) x1 <- data.frame(table(var))

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