similar to: cbind formula definition

Dear R-users, I keep getting an ERROR saying " Error in model.matrix.default(formula, mframe) : model frame and formula mismatch in model.matrix() " when i fit poLCA with more than 63 variables. Below are the details. I am trying to do a Latent Class Analysis using poLCA. My data set contains binary scores of, for instance, 200 students on 100 items. These numbers could even be more

Hi there, I have the following problem Four data frames exist: data1 data2 data3 data4 Now I want to write a loop and temporarily store the data1, data2, data3, data4 in a variable called data. I tried the following... for (i in 1:4) { Data <- paste("data",i,sep="") ... .. } but it doesn't function. I think the problem is the definition of the mode of the pasted

Hi everyone, I'm using the poLCA module for some analysis work. Basically, Im using the command >poLCA(f, data=response,nclass=2) This returns a poLCA object (a list) >From this data, I need to retrieve certain indexes, such as[[5]] as seen below. [[5]] Pr(1) Pr(2) Pr(3) class 1: 0 1.0 0.0 class 2: 0 0.5

Hello, I want to perform a latent class analysis using poLCA package. My formula is: substances <- cbind(subs1, subs2, subs3, subs4, subs5, subs6) ~ gender+age+education+income+occupation+urban+dbehavior+incarceration+treatment+depression+alcriteria I want to include sample weights in the model, I have read that poLCA does not take into account weights, but when I introduce them, it seems

Hi there, Does anyone know a method to calculate the number of different patterns in a given data frame. The variables are of polytomous type and not binary (for the latter i found a package called "countpattern" which unfortunately only functions for binary variables). V1 V2 V3 0 3 1 1 2 0 1 2 0 So, in this case, i would like to get "2" as output. Thanks

Dear all, I'm an Sri Lankan undergraduate student. I'm also a total newbie to R. My aim is to use the poLCA package to do a latent class analysis. I found the documentation very helpful, but need to make a small clarification that has stumped me awhile. In my work, I need to make provision for conditional dependence. I'm told that poLCA lets you do that. Unfortunately, I

Dear Mr. Therneau, Mr. Fox, or to whoever, who has some time... I don't find a solution to use the "survfit" function (package: survival) for a defined pattern of covariates with a Cox-Model including a time dependent interaction term. Somehow the definition of my "newdata" argument seems to be erroneous. I already googled the problem, found many persons having the

The following is evidence of what is surely an undesirable feature. The issue is the handling, in calls to model.frame(), of an explanatory variable that has been derived as an unclassed factor. (Ross Darnell drew this to my attention.) ## Data are slightly modified from p.191 of MASS > worms <- data.frame(sex=gl(2,6), Dose=factor(rep(2^(0:5),2)), +

What's the best approach to running latent class analysis with R? I've downloaded both randomLCA and poLCA packages, but I am interesting in running a standard LCA with individual records (not frequency table) as input data. Wen Gu John Jay College of Criminal Justice445 West 59 StreetNew York, NY 10029 wgu@gc.cuny.edu _________________________________________________________________

Hi, I am working with -poLCA- which uses the package -scatterplot3d- in order to work. However, I work on a restricted server that is not online. The R-version installed there is 2.6. However, running library(poLCA) after installing the most recent -scatterplot3d- on the server leaves me with R telling me that scatterplot3d needs R-version 2.7 or a more recent versions. Since I cannot install the

My problem is I have data with both categorial and numerical data, currently only the categorical number contains missing data, was wondering do I make a new dataframe containing only the categorical columns? How would you use Latent Class Model specifically poLCA to impute the missing data? http://www.sscnet.ucla.edu/polisci/faculty/lewis/pdf/poLCA-JSS-final.pdf The reason why I chose not to

Hi everyone? Does anybody know how can I estimate a Latent class multinomial (or conditional) logit using R? I have tried flexmix, poLCA, and they do not seem to support this model. thanks in advance adan -- View this message in context: http://r.789695.n4.nabble.com/Latent-class-multinomial-or-conditional-logit-using-R-tp4230083p4230083.html Sent from the R help mailing list archive at

Dear R users, I would like to perform Latent Class Logit Models for the analysis of choice experiments in environmental valuation. This kind of analysis is usually performed with NLogit Software (http://www.limdep.com). I attach the results I usually obtain using NLogit and NLogit model specifications. For Random parameter models and Logit Models I usually perform my analysis with the package

Hi, I am running some latent class models with civariates using poLCA in R, and have a question the meaning of some error messages. My model has 8 y-variables (dichotomous) and 2 x-variables (also dichotomous), and I am trialling different numbers of classes. When I run a 2-class model, the model runs fine without any error messages. I used the nrep = 10 option to avoid local maxima. Class %

Dear List, I?m trying to figure out *memory page sharing***using HVM on Xen 4.0.1, but I?m not really able to find some useful information or sources about this issue. I just found two links porviding these sources: http://knol.google.com/k/learning-grant-tables http://blog.chinaunix.net/space.php?uid=20286427&do=blog&id=109114 Both are not working.. I think they are made for an older

Hi, I am using wxRuby 1.9.7 and Ruby 1.8.6. I have built an extension in C++ and used swig to load it into Ruby interpreter, which works great. The extension is an often long executing algorithm, and I''ve noticed that if I fork a new Ruby thread in button clicked event, and in this thread run the algorithm the application freezes for the time of execution. I have looked through the

HI there, I want to add a legend to a plot using the density and angle argument, so patterns with lines in different angles are used in the plot and should be referred to. When I use default settings, the filled boxes are too small. With the cex argument I can enlarge the whole legend, but then the text gets too big. How could I just increase the size of the single boxes and not the text. I

Hello R users,  I am trying to obtain a direct adjusted survival curve. I am sending my whole code (see below). It's basically the larynx cancer data with Stage 1-4. I am using the cox model using coxph option, see the fit3 coxph. When I use the avg.surv option on fit3, I get the following error: "fits<-avg.surv(fit3, var.name="stage.fac", var.values=c(1,2,3,4), data=larynx)

Well, I would suggest using the code already in place in the survival package. Here is my code for your problem. I'm using a copy of the larynx data as found from the web resources for the Klein and Moeschberger book. larynx <- read.table("larynx.dat", skip=12, col.names=c("stage", "time", "age", "year",

Hi everyone, I've been wondering how to correctly build clang/LLVM for the AVR target architecture. Unfortunately documentation is very scarce (or outdated or I didn't find it) and while I've been able to build clang/LLVM for AVR I'm still falling short of compiling an actual binary for the MCU. Here are the steps I've undertaken so far: git clone