[R] Error: attempt to apply non-function

[Duncan Murdoch]

  On 23/09/2010 8:54 AM, 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?

It would help a lot if you posted the actual error message, and the 
results of a call to traceback() just after it.  We can't run your code 
(we don't have x0.trial01 and don't know how you created it).

Duncan Murdoch

> 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.