set.seed(0xBEEF)
n <- 300
err <- 2
x <- rnorm(n, 0, 5)
y <- rnorm(n, 0, 5)
df <- data.frame(x, y)
df <- within(df, {
diff <- x-y
absdiff <- abs(x-y)
sqdiff <- (x-y)^2
z.sq <- sqdiff + rnorm(n, 0, err)
})
r1 <- RSA(z.sq~x*y, df)
b1 <- bootRSA(r1, model="SSD", R=5000, parallel="multicore", ncpus=2)
r1.boot.CI <- CI.boot(b1) # compute percentile confidence intervals and percentile p-value of the bootstrapped values
Run the code above in your browser using DataLab