Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- a <- rbinom(100, 1, .5)
- a <- sorted(a)
- holder <- data.frame(1:length(a), a)
- ols <- lm(holder$a ~ holder$X1.length.a.)
- resid <- holder$a - predict(ols)
- plot(predict(ols),resid)
- normErr.Small <- lapply(1:100, function(x) {
- x + rnorm(n = 1, 1, sd = 1)
- })
- dat <- data.frame(x = 1:length(normErr.Small),y = unlist(normErr.small))
- ols.small <- lm(y ~ x, data = dat)
- plot(ols.small, which = 1, main = "normal 1 1")
- #larger variance
- normErr.Large <- lapply(1:100, function(x) {
- x + rnorm(n = 1, 1, sd = 10000)
- })
- dat <- data.frame(x = 1:length(normErr.Large),y = unlist(normErr.Large))
- ols.Large <- lm(y ~ x, data = dat)
- plot(ols.Large, which = 1, main = "normal 1 10000")
- normErr.Off <- lapply(1:100, function(x) {
- x + rnorm(n = 1, 10000, sd = 1)
- })
- dat <- data.frame(x = 1:length(normErr.Off),y = unlist(normErr.Off))
- ols.Off <- lm(y ~ x, data = dat)
- plot(ols.Off, which = 1, main = "normal 10000 1")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement