[R] non-numeric argument to binary operator

[Richard M. Heiberger]

It looks like you are attempting to divide a number by a formula.

> 5 / (a ~ b)
Error in 5/(a ~ b) : non-numeric argument to binary operator


On Fri, Jan 20, 2017 at 12:18 PM, Zeki ÇATAV <zcatav at gmail.com> wrote:

> Hi,
> I'm working on following dataset.
>
>
> > dput(uu5)
> structure(list(grup = structure(c(1L, 1L, 2L, 2L), .Label = c("1",
> "2"), class = "factor"), bw_grp = structure(c(1L, 2L, 1L, 2L), .Label =
> c("0-999",
> "1000+"), class = "factor"), deaths = c(6L, 13L, 1L, 2L), pop = c(26L,
> 67L, 41L, 93L), bw = c(500L, 1500L, 500L, 1500L), rates =
> c(23.0769230769231,
> 19.4029850746269, 2.4390243902439, 2.1505376344086), pct =
> c(0.279569892473118,
> 0.720430107526882, 0.305970149253731, 0.694029850746269), logexpo =
> c(3.25809653802148,
> 4.20469261939097, 3.71357206670431, 4.53259949315326), grup1 = c(1,
> 1, 0, 0)), .Names = c("grup", "bw_grp", "deaths", "pop", "bw",
> "rates", "pct", "logexpo", "grup1"), class = c("grouped_df",
> "tbl_df", "tbl", "data.frame"), row.names = c(NA, -4L), vars = list(
>     grup), labels = structure(list(grup = structure(1:2, .Label = c("1",
> "2"), class = "factor")), row.names = c(NA, -2L), class = "data.frame",
> vars = list(
>     grup), drop = TRUE, .Names = "grup"), indices = list(0:1,
>     2:3), drop = TRUE, group_sizes = c(2L, 2L), biggest_group_size = 2L)
>
> When I try;
> uu4 <- mutate(uu4, logexpo = log(pop), grup1 = as.numeric(group == 1 ))
> fit <- glm(deaths ~ agegrp + grup1 + offset(log(pop)), data = uu4, family =
> poisson)
> coef(summary(fit))
> exp(coef(fit)["grup1"])
> data.frame(deviance(fit), df=fit$df.residual)
> deviance(fit)/(fit$null.deviance ~ deviance(fit))
> this code, I get an error related to last line:
>
> Error in deviance(fit)/(fit$null.deviance ~ deviance(fit)) :
>   non-numeric argument to binary operator
>
> What's wrong? How can I correct this problem?
> Thanks in advance.
>
> PS: R version 3.3.2 (2016-10-31) -- "Sincere Pumpkin Patch"
>
> --
> Zeki Çatav
> http://zekicatav.com <http://zekicatav.tyih.gov.tr>
> twitter.com/zcatav
>
>         [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> 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.

	[[alternative HTML version deleted]]