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
>
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
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
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
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
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
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
7 matches
Mail list logo
