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 <[email protected]> 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]] > > ______________________________________________ > [email protected] 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]] ______________________________________________ [email protected] 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.
