Skip to content

Commit

Permalink
Rebuild under knitr 1.48
Browse files Browse the repository at this point in the history
  • Loading branch information
remlapmot committed Jul 7, 2024
1 parent e310bb4 commit 0e4ea9a
Show file tree
Hide file tree
Showing 45 changed files with 4,638 additions and 5,569 deletions.
Binary file modified data/fig1.dta
Binary file not shown.
Binary file modified data/fig2.dta
Binary file not shown.
Binary file modified data/fig3.dta
Binary file not shown.
Binary file modified data/nhefs-formatted.dta
Binary file not shown.
Binary file modified data/nhefs-highprice.dta
Binary file not shown.
Binary file modified data/nhefs-ps.dta
Binary file not shown.
Binary file modified data/nhefs-wcens.dta
Binary file not shown.
Binary file modified data/nhefs_std.dta
Binary file not shown.
Binary file modified data/nhefs_std1.dta
Binary file not shown.
Binary file modified data/nhefs_std2.dta
Binary file not shown.
Binary file modified data/nhefs_surv.dta
Binary file not shown.
Binary file modified data/observe.mmat
Binary file not shown.
24 changes: 0 additions & 24 deletions docs/11-why-model-r.md
Original file line number Diff line number Diff line change
Expand Up @@ -24,15 +24,9 @@ plot(A, Y)
summary(Y[A == 0])
#> Min. 1st Qu. Median Mean 3rd Qu. Max.
#> 10.0 27.5 60.0 67.5 87.5 170.0
```

``` r
summary(Y[A == 1])
#> Min. 1st Qu. Median Mean 3rd Qu. Max.
#> 50.0 105.0 160.0 146.2 185.0 220.0
```

``` r

A2 <- c(1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4)
Y2 <- c(110, 80, 50, 40, 170, 30, 70, 50, 110, 50, 180,
Expand All @@ -47,21 +41,12 @@ plot(A2, Y2)
summary(Y2[A2 == 1])
#> Min. 1st Qu. Median Mean 3rd Qu. Max.
#> 40.0 47.5 65.0 70.0 87.5 110.0
```

``` r
summary(Y2[A2 == 2])
#> Min. 1st Qu. Median Mean 3rd Qu. Max.
#> 30 45 60 80 95 170
```

``` r
summary(Y2[A2 == 3])
#> Min. 1st Qu. Median Mean 3rd Qu. Max.
#> 50.0 95.0 120.0 117.5 142.5 180.0
```

``` r
summary(Y2[A2 == 4])
#> Min. 1st Qu. Median Mean 3rd Qu. Max.
#> 150.0 187.5 205.0 195.0 212.5 220.0
Expand Down Expand Up @@ -107,15 +92,9 @@ summary(glm(Y3 ~ A3))
#> AIC: 170.43
#>
#> Number of Fisher Scoring iterations: 2
```

``` r
predict(glm(Y3 ~ A3), data.frame(A3 = 90))
#> 1
#> 216.89
```

``` r

summary(glm(Y ~ A))
#>
Expand Down Expand Up @@ -168,9 +147,6 @@ summary(mod3)
#> AIC: 170.39
#>
#> Number of Fisher Scoring iterations: 2
```

``` r
predict(mod3, data.frame(cbind(A3 = 90, Asq = 8100)))
#> 1
#> 197.1269
Expand Down
Loading

0 comments on commit 0e4ea9a

Please sign in to comment.