[R] Error: attempt to apply non-function

[peter dalgaard]

On Sep 23, 2010, at 14:54 , Maas James Dr (MED) wrote:

> This code worked fine for me, then did some cleaning up of formatting using ESS (Emacs) and now I get this error, no idea what is causing it, all the brackets/parentheses seem to be balanced.  What have I done wrong?


Look for ") (" -- missing a "+", presumably. As a general matter, try using tools like traceback() or options(recover=...) to catch such things in the act. They are quite hard to spot by visual inspection.

-pd 

> Thanks
> 
> Jim
> 
> 
> 
> p0.trial01 <- 0.25
> TruOR01 <- 0.80
> num.patients.01 <- 50
> num.trials.01 <- 5
> LOR01.het.in <- 0.00
> num.sims <- 1
> 
> simLOR01 <- vector(length=num.trials.01)
> simLORSE01 <- vector(length=num.trials.01)
> simOR01 <- vector(length=num.trials.01)
> trialnum01 <- vector(length=num.trials.01)
> 
> x0count <- vector(length = num.trials.01)
> x1count <- vector(length = num.trials.01)
> 
> ## Trial 1, comparison of treatment 1 with treatment 0
> 
> for (i in 1:num.trials.01) {
> 
>    het01 <- rnorm(1,0,LOR01.het.in)
>    log.odds.ratio.01 <- log(TruOR01) + het01
>    odds.ratio.01 <- exp (log.odds.ratio.01)
> 
>    p1.trial01 <- (p0.trial01 * odds.ratio.01) / (1 - p0.trial01 +
>                                                  (p0.trial01 * odds.ratio.01))
> 
>    x0.trial01[i] <-  rbinom(1,num.patients.01,p0.trial01)
>    x1.trial01[i] <-  rbinom(1,num.patients.01,p1.trial01)
> 
>    trialnum01[i] <- paste ( c ("trial01-"), i, sep="")
> 
>    simLOR01[i] <- log (x1.trial01 * (num.patients.01 - x0.trial01) /
>                        ((num.patients.01 - x1.trial01) * x0.trial01))
> 
>    simLORSE01[i] <- sqrt ( (1/x0.trial01) + (1/(num.patients.01 -
>                                                 x0.trial01)) (1 / x1.trial01) + (1 / (num.patients.01 -
>                                                                                       x1.trial01)))
> 
>    simOR01[i] <- (x1.trial01 * (num.patients.01 - x0.trial01) /
>                   ((num.patients.01 - x1.trial01) * x0.trial01))
> 
> }
> 
> 
>    ## Output all results to a data.frame called results
> 
>    results <- data.frame (tn1 = trialnum01, SimOR01 = simOR01,
>                           SimLOR01 = simLOR01, SimLORSE01 = simLORSE01,
>                            p1trial = p1.trial01, x0count, x1count )
> 
> set.seed(9321685)
> 
> results
> 
> rm(list=ls())
> 
> 
> ===============================
> Dr. Jim Maas
> University of East Anglia
> 
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help at 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.

-- 
Peter Dalgaard
Center for Statistics, Copenhagen Business School
Solbjerg Plads 3, 2000 Frederiksberg, Denmark
Phone: (+45)38153501
Email: pd.mes at cbs.dk  Priv: PDalgd at gmail.com