http://www.commres.net/ 2026-04-17T08:24:43+00:00 http://www.commres.net/ http://www.commres.net/_media/wiki/logo.png text/html 2023-10-26T06:28:06+00:00 Anonymous (anonymous@undisclosed.example.com) http://www.commres.net/regression/note?rev=1698301686&do=diff #setwd("user location") #Working directory set.seed(123) #Standardizes the numbers generated by rnorm; see Chapter 5 n <- 100 #Number of participants; graduate students x <- rnorm(n, 175, 7) #IV; hours since dawn m <- 0.7*x + rnorm(N, 0, 5) #Suspected mediator; coffee consumption y <- 0.4*m + rnorm(N, 0, 5) #DV; wakefulness md <- data.frame(x, m, y) attach(md) cor(md) m01 <- lm(y~x, data=md) summary(m01) m02 <- lm(m~x, data=md) summary(m02) m03 <- lm(y~x+m, data=md) summary(m03) install.packa… text/html 2024-09-30T01:41:34+00:00 Anonymous (anonymous@undisclosed.example.com) http://www.commres.net/regression/note2?rev=1727660494&do=diff rm(list = ls()) dat <- read.csv("http://commres.net/wiki/_media/regression01-bankaccount.csv") head(dat) attach(dat) mod <- lm(bankaccount ~ income, data = dat) summary(mod) df.tot <- length(bankaccount) - 1 ss.tot <- var(bankaccount)*df.tot var.tot <- var(bankaccount) df.x <- length(income)-1 ss.x <- var(income)*df.x m.y <- mean(bankaccount) m.x <- mean(income) sp.xy <- sum((income-m.x)*(bankaccount-m.y)) cov.xy <- sp.xy / df.tot cov.xy cov(income, bankaccount) sd.x <- sd(income) sd.y <…