I always have trouble with "by" myself. I don't usually get an answer
as quickly as I did here. I usually give up and move to "aggregate" or
"ave".
First I limited the data given to "by" in the first argument to just
the variable under question. Then by the time R sends the information
to your function it may not be named what it was in the global
environment, so I packaged those two functions "runs.test" and
"as.factor" inside a packaging function that converted its name to
"x". Seems to give possibly sensible results. See if this works for you.
> by(exampledata$var, exampledata$interviewer,
function(x) {runs.test(as.factor(x))} )
exampledata$interviewer: a
Runs Test
data: as.factor(x)
Standard Normal = -0.8823, p-value = 0.3776
alternative hypothesis: two.sided
----------------------------------------------------------------------------------------------
exampledata$interviewer: b
Runs Test
data: as.factor(x)
Standard Normal = 0, p-value = 1
alternative hypothesis: two.sided
--
David Winsemius
On Jan 14, 2009, at 3:13 AM, christiaan pauw wrote:
by(exampledata, exampledata$interviewer,
runs.test(as.factor(exampledata$var)))
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
