[R] effect.clm

[David Barron]

I'm having a problem using the xlevels option with effect (or Effect) with
objects of class clm, which seems to have no impact on the output (whereas
it does for models of class polr).  Can anyone suggest how to fix this?

Thanks,
David

Here's an example:

mod.wvs <- MASS::polr(poverty ~ gender + religion + degree +
country*poly(age,3), data=WVS)

 Effect('age', mod.wvs, xlevels = list(age = 1:10))

Re-fitting to get Hessian


age effect (probability) for Too Little
age
        1         2         3         4         5         6         7
  8         9        10
0.8403661 0.8270171 0.8133710 0.7995001 0.7854780 0.7713778 0.7572714
0.7432277 0.7293120 0.7155850

age effect (probability) for About Right
age
        1         2         3         4         5         6         7
  8         9        10
0.1295679 0.1399772 0.1505377 0.1611853 0.1718564 0.1824888 0.1930231
0.2034039 0.2135803 0.2235068

age effect (probability) for Too Much
age
         1          2          3          4          5          6
 7          8          9         10
0.03006598 0.03300572 0.03609134 0.03931456 0.04266563 0.04613341
0.04970548 0.05336837 0.05710766 0.06090822

But with clm:

+                 data=WVS)
> Effect('age', mod.wvs, xlevels = list(age = 1:10))

age effect (probability) for Too Little
age
       20        30        40        50        60        70        80
 90
0.5980510 0.5268304 0.4936603 0.4803375 0.4674829 0.4359217 0.3684641
0.2589983

age effect (probability) for About Right
age
       20        30        40        50        60        70        80
 90
0.3031130 0.3453414 0.3629620 0.3696094 0.3757686 0.3897338 0.4129734
0.4227230

age effect (probability) for Too Much
age
        20         30         40         50         60         70
80         90
0.09883598 0.12782815 0.14337772 0.15005316 0.15674848 0.17434449
0.21856247 0.31827868

	[[alternative HTML version deleted]]