Message-ID: <71BFAC6FE64A9C44A4DDAF2F786B79A80947233D@DB3PRD0410MB380.eurprd04.prod.outlook.com>
Date: 2013-01-22T17:16:58Z
From: Jason Musil
Subject: plot.mob() fails with cut() error "'breaks' are not unique"
DeaR all,
I am using mob() for model based partitioning, with a dichotomous variable (participant's correct/incorrect response to a test item) regressed onto a continuous predictor related to a given property of the test item. Although this variable is continuous, the value of this variable for many items in this particular analysis is 0. The partitioning criterion is self-reported ability in a related area.
> mob1 <- mob(
correct ~ circular.mean | srp.dimension,
control=mob_control(alpha=.001),
model=glinearModel,
family=binomial()
)
> plot(mob1)
Error in cut.default(x, breaks = breaks, include.lowest = TRUE) :
'breaks' are not unique
The same persists if I specify either a desired number of breaks, or explicit breakpoints (e.g. breaks=3 or breaks=c(-0.1,0.1,0.5)). I guess this is to do with the funny distribution of the predictor variable, but I'm not sure what to do about it.
Many thanks and apologies if this doesn't fit the mailing list---it is my first posting!
Jason Musil