Ludvig Londos och Kalle Palm (laddy_ludde@msn.com, kalle.palm@gmail.com)
mean(time ~ sex, data = swim)
## F M
## 65.19 54.66
min(time ~ sex, data = swim)
## F M
## 53.52 47.84
Svar: 65.2 och 53.5
early = subset(swim, year < 1920)
mean(early$time)
## [1] 73.84
min(early$time)
## [1] 61.4
Svar: 73.8 och 61.4
slow = subset(swim, time > 60)
mean(slow$time)
## [1] 69.62
min(slow$time)
## [1] 60.2
Svar: 69.6 och 60.2
mod = mm(wage ~ sector, data = CPS85)
mod
##
## Groupwise Model Call:
## wage ~ sector
##
## Coefficients:
## clerical const manag manuf other prof sales
## 7.42 9.50 12.70 8.04 8.50 11.95 7.59
## service
## 6.54
9.5
12.7
Service
sd(fitted(mod))
## [1] 2.197
Svar: 2.20
sd(resid(mod))
## [1] 4.646
Svar: 4.65
mod1 = mm(wage ~ 1, data = CPS85)
mod2 = mm(wage ~ sector, data = CPS85)
var(fitted(mod1))
## [1] 0
var(fitted(mod2))
## [1] 4.825
Svar: mod2
var(resid(mod1))
## [1] 26.41
var(resid(mod2))
## [1] 21.59
Svar: mod1
mean(w$wage)
## [1] 9.024
Svar: 9.02
wageSex = mm(wage ~ sex, data = w)
wageSex
##
## Groupwise Model Call:
## wage ~ sex
##
## Coefficients:
## F M
## 7.88 9.99
Svar: 7.88
gift = subset(w, married == "Married")
mean(gift$wage)
## [1] 9.398
Svar: 9.40
ogiftaKvinnor = subset(w, (married == "Single" & sex == "F"))
mean(ogiftaKvinnor$wage)
## [1] 8.26
Svar: 8.26
sd(g$height)
## [1] 3.583
Svar: 3.58
mod3 = mm(height ~ 1, data = g)
res = resid(mod3)
sd(res)
## [1] 3.583
Svar: 3.58
mod4 = mm(height ~ sex, data = g)
res1 = resid(mod4)
sd(res1)
## [1] 2.508
Svar: 2.51
Svar: Den andra modellen (kallad mod1 i uppgifterna)