[R] multilevel sampling weight

2013-11-18 Thread Pablo Menese Camargo
Anyone is working with the upgrade of any multilevel package for apply sampling weights? Thanks [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting gu

Re: [R] foreign package problem with stata 13

2013-10-02 Thread Pablo Menese Camargo
dataset in stata. > > I don't think foreign library in R supports stata 13 yet. > > Richard > > On Wed, Oct 2, 2013 at 3:11 PM, Pablo Menese Camargo > wrote: > > foreign package does not support dataset saved at stata 13. > > anyone knows any wayt to make

[R] foreign package problem with stata 13

2013-10-02 Thread Pablo Menese Camargo
foreign package does not support dataset saved at stata 13. anyone knows any wayt to make it works? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] SEM polychoric lavaan

2013-08-12 Thread Pablo Menese Camargo
I have this mcfa<- "impulsivity=~imp1r+imp2+imp3+imp4 physical=~phys1+phys2+phys3+phys4 risky=~risk1+risk2+risk3+risk4 selfish=~self1+self2+self3+self4 simple=~simp1+simp2+simp3+simp4r temper=~temp1+temp2+temp3+temp4 control=~impulsivity+physical+risky+selfish+simple+temper" where: the model is a

Re: [R] SEM polychoric

2013-08-12 Thread Pablo Menese Camargo
g a > question. > Regards, > > José > > Prof. José Iparraguirre > Chief Economist > Age UK > > Profesor de Economía > Universidad de Morón > Morón, Buenos Aires, Argentina > > > > -Original Message----- > From: r-help-boun...@r-project.org [

[R] SEM polychoric

2013-08-11 Thread Pablo Menese Camargo
I have this mcfa<- "impulsivity=~imp1r+imp2+imp3+imp4 physical=~phys1+phys2+phys3+phys4 risky=~risk1+risk2+risk3+risk4 selfish=~self1+self2+self3+self4 simple=~simp1+simp2+simp3+simp4r temper=~temp1+temp2+temp3+temp4 control=~impulsivity+physical+risky+selfish+simple+temper" where: the model is a

[R] eRm package, RSM: Error in solve.default(parets$hessian) :

2013-05-21 Thread Pablo Menese Camargo
I want to make an index using IRT, I use eRm package. I have a data frame (24 items, 4 categories, 427 cases) that works if I use PCM (partial credit model), but, when I run RSM (rating scale model) this appear: Error in solve.default(parest$hessian) : Lapack routine dgesv: system is exactly sin

[R] Confirmatory Factor Analysis

2013-03-10 Thread Pablo Menese
Hi, I'm trying to perform a hierachical, second order CFA. That's the thing that I need to leave AMOS. I found some sim.hierarchical and omega stuff, but nothing clear enough. Can anyone help me? I just need a simple and clear manual. Best, Pablo. [[alternative HTML version deleted]] _

Re: [R] issue creating a subset

2013-03-04 Thread Pablo Menese
you create the matrix using > cbind, as in my matrix3 example. > > -Don > > > -- > Don MacQueen > > Lawrence Livermore National Laboratory > 7000 East Ave., L-627 > Livermore, CA 94550 > 925-423-1062 > > > > > > On 2/28/13 5:53 PM, "Pablo M

[R] issue creating a subset

2013-02-28 Thread Pablo Menese
I'm performing item response theory with eRm packages I am excluding the persons that doesn't fit in the infit/outfit persons. for that I created a condition. then I have to create a new subset or matrix but with the condition. So: ORIGINAL matrix<-cbind(item1, item2, item3, item4) IF I PERFORM

Re: [R] Issue attaching a dataset

2013-02-28 Thread Pablo Menese
thanks a lot. works!!! On Tue, Feb 26, 2013 at 7:58 AM, Prof Brian Ripley wrote: > On 26/02/2013 09:55, Milan Bouchet-Valat wrote: > >> Le lundi 25 février 2013 à 20:26 -0300, Pablo Menese a écrit : >> >>> I use to work whit stata dataframe, so, when I use R I type re

[R] Issue attaching a dataset

2013-02-25 Thread Pablo Menese
I use to work whit stata dataframe, so, when I use R I type read.dta Until today I do that without any problem, after type: mydata<-read.dta("C:/dropbox/dataframe.dta") attach(mydata) Everything works great... but today, when I typed: mydata<-read.dta("C:/dropbox/dataframe.dta") attach(mydata)

Re: [R] problem with svyglm

2012-11-28 Thread Pablo Menese
81 .0410122 .4743798 .6359363 So, I thing that the problem is in the survey design... On Tue, Nov 27, 2012 at 11:49 PM, David Winsemius wrote: > > On Nov 27, 2012, at 2:31 PM, Pablo Menese wrote: > > Sorry, it send it alone... >> >> When I use it: >> >

Re: [R] Problems with weight

2012-11-28 Thread Pablo Menese
: > Le mardi 27 novembre 2012 à 18:33 -0300, Pablo Menese a écrit : > > I can't ... I don't know why but I can't > > > > When I use it: > > > > logit <- glm(bach ~ egp4 + programa, weight=wst7, > > family=quasibinomial(link"logit")) >

Re: [R] Problems with weight

2012-11-27 Thread Pablo Menese
> > http://journal.r-project.org/archive/2009-2/RJournal_2009-2_Damico.pdf > > once you're ready to run the regression, use svyglm() instead of glm() and > drop the weights argument (since it will already be part of the survey > design) :) > > > > On Fri, Nov 23, 2012

Re: [R] problem with svyglm

2012-11-27 Thread Pablo Menese
Sorry, it send it alone... When I use it: logit <- glm(bach ~ egp4 + programa, weight=wst7, family=quasibinomial(link"logit")) I reach the same betas that in STATA, but the hypothesis test, the t value, and the std. error is different. I think that the solution can't be so far from this...

Re: [R] problem with svyglm

2012-11-27 Thread Pablo Menese
I colud not, even without attach the dataset. The thing is, when I use this On Fri, Nov 23, 2012 at 5:56 PM, David Winsemius wrote: > > On Nov 23, 2012, at 12:08 PM, Pablo Menese wrote: > > I have this problem. >> >> test <- svydesign(id=~1,weights=~peso) >&g

[R] Problems with weight

2012-11-23 Thread Pablo Menese
Until a weeks ago I used stata for everything. Now I'm learning R and trying to move. But, in this stage I'm testing R trying to do the same things than I used to do in stata whit the same outputs. I have a problem with the logit, applying weights. in stata I have this output . svy: logit bach job

[R] problem with svyglm

2012-11-23 Thread Pablo Menese
I have this problem. test <- svydesign(id=~1,weights=~peso) logit <- svyglm(bach ~ job2 + mujer + egp4 + programa + delay + mdeo + str + evprivate, family=binomial,design=test) then appear: Error in svyglm.survey.design(bach ~ job2 + mujer + egp4 + programa + : all variables must be in desig