Re: [R] effect function in the effects package

2011-10-20 Thread Xu Jun
Dear Prof. Fox, I just picked up R not long ago, and apologize that I am not that familiar with some basics. I am trying to replicate what I can do with Stata in R. Thanks for all your help! On Thu, Oct 20, 2011 at 2:47 PM, John Fox wrote: > Dear Xu Jun, > > On Thu, 20 Oct 2011 14:41:30 -0400 >

Re: [R] effect function in the effects package

2011-10-20 Thread John Fox
Dear Xu Jun, On Thu, 20 Oct 2011 14:41:30 -0400 Xu Jun wrote: > Dear Professor Fox, > > Now I got it. It all comes from my unfamiliarity with the effect > function. I forgot the c part in the given.values option, plus it > looks like plot(effect()) does not like factor(warm) in the polr > funct

Re: [R] effect function in the effects package

2011-10-20 Thread Xu Jun
Dear Professor Fox, Now I got it. It all comes from my unfamiliarity with the effect function. I forgot the c part in the given.values option, plus it looks like plot(effect()) does not like factor(warm) in the polr function. So here are the two working lines: ordwarm2$warm2 <- as.factor(ordwarm2

Re: [R] effect function in the effects package

2011-10-20 Thread John Fox
Dear Xu Jun, It's really not possible for me to know the source of the error without a complete, reproducible example, but if I had to guess, I'd guess that there's a scoping problem of some sort, with the value of yr89 coming from somewhere that you don't expect. How about fitting the original

Re: [R] effect function in the effects package

2011-10-20 Thread Xu Jun
Dear Professor Fox, I didn't include codes in between these two commands (polr and plot(effect()). I am still trying to work it out. Maybe it's related to how I coded some of the factor variables as I got some successful runs after I changed the coding for some factor variables. I will keep workin

Re: [R] effect function in the effects package

2011-10-20 Thread Xu Jun
Dear Professor Fox, Thanks a lot! That is an embarrassing error. After I cleaned up and simplified my codes, and ran the following two lines: myologit <- polr(factor(warm) ~ yr89 + male + white + age + ed + prst, + data=ordwarm2, method=c("logistic")) plot(effect("age", myologit, x

Re: [R] effect function in the effects package

2011-10-20 Thread John Fox
Dear Xu Jun, I'm not sure whether this is the source of the error, but it may help to spell the xlevels argument correctly (it is not "xlevles"). I hope this helps, John John Fox Sen. William McMaster Prof. of Social Statistics Department of Soc

[R] effect function in the effects package

2011-10-20 Thread Xu Jun
Dear r-help listers, I am using effects to produce an effect plot after the proportional odds logistic regression model. There is no problem for me to estimate the model, but when it comes to the graphing, I was stuck. see the codes below: #