Displaying 1 result from an estimated 1 matches for "unempw".
Did you mean:
unemp
2012 Apr 11
1
Problem with effects package
...3-17
loaded via a namespace (and not attached):
[1] cluster_1.14.2 digest_0.5.2 tools_2.15.0
# copy and paste the following:
library( catdata )
data( unemployment )
unempt <- unemployment
unempt$durbin <- unempt$durbin - 1
library( reshape )
unempt <- melt( table( unempt ) )
unempw <- cast( unempt, age ~ durbin )
names( unempw ) <- c( 'age', 'short', 'long' )
modt <- glm( durbin ~ age, weights = value, family = binomial, data = unempt )
modw <- glm( cbind( short, long ) ~ age, family = binomial, data = unempw )
library( effects )
modt.ef &l...