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?

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

