Description Usage Arguments Value Examples
Simulates data from a linear model using the sim.lm() function, each time applying a permutation test to test the null hypothesis that the slope is zero.
1 2 | sim.perm.B(n, nsim, a, b, nperm = 500, sigma.disturb = 1, mu.x = 8,
sigma.x = 2, rdisturb = rnorm, rx = rnorm, het.coef = 0)
|
n |
The number of pairs of observations to simulate. |
nsim |
The number of simulations to conduct. |
a |
The intercept parameter of the linear regression model to be simulated. |
b |
The slope parameter of the linear regression model to be simulated. |
nperm |
The number of permutations to use per simulation to conduct the permutation test. |
sigma.disturb |
The variance of the disturbances in the model y = a + b*x + disturbance. |
mu.x |
The expectation of the x values. |
sigma.x |
The variance of the x values. |
rdisturb |
A function for drawing the random disturbances. rnorm() is the default, which makes the disturbances normally distributed, but you can use any function for random number generation with first argument the sample size, second argument the expectation, and third argument the standard deviation. |
rx |
A function for drawing the random x values. rnorm() is the default, which makes x normally distributed, but you can use any function for random number generation with first argument the sample size, second argument the expectation, and third argument the standard deviation. |
het.coef |
A number introducing some heteroscedasticity (i.e. non-constant variance) to the disturbances. If het.coef = 0 (the default), then the disturbances have constant variance. If it is positive, then the standard deviation of the disturbances increases with x; if negative, then the standard deviation of the disturbances decreases with x. |
A vector of p values of length nsims, one from each permutation test of the hypothesis that the slope is zero.
1 2 | ps <- sim.perm.B(10, 200, a = 3, b = .1, nperm = 200)
mean(ps < .05) #power at .05 significance level
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.