Skip to content
This repository has been archived by the owner on Jun 27, 2024. It is now read-only.

Commit

Permalink
Merge branch 'main' into pa/partial_within
Browse files Browse the repository at this point in the history
  • Loading branch information
palday authored Sep 12, 2023
2 parents 47b4e4b + 7ebe238 commit 9062918
Show file tree
Hide file tree
Showing 4 changed files with 9 additions and 12 deletions.
5 changes: 1 addition & 4 deletions kb07.qmd
Original file line number Diff line number Diff line change
Expand Up @@ -10,7 +10,6 @@ Begin by loading the packages to be used.
#| output: false
using AlgebraOfGraphics
using CairoMakie
using DataFrameMacros
using DataFrames
using MixedModels
using ProgressMeter
Expand All @@ -19,7 +18,6 @@ using SMLP2023: dataset
CairoMakie.activate!(; type="svg")
import ProgressMeter
ProgressMeter.ijulia_behavior(:clear)
```

Expand Down Expand Up @@ -108,8 +106,7 @@ The pseudo-random number generator is initialized to a fixed value for reproduci

```{julia}
Random.seed!(1234321)
hide_progress = true
kbm02samp = parametricbootstrap(2000, kbm02; hide_progress)
kbm02samp = parametricbootstrap(2000, kbm02)
kbm02tbl = kbm02samp.tbl
```

Expand Down
10 changes: 5 additions & 5 deletions kkl15.qmd
Original file line number Diff line number Diff line change
Expand Up @@ -38,21 +38,19 @@ Here we provide some of the corresponding analyses with _MixedModels.jl_ and a m
#| output: false
using Arrow
using AlgebraOfGraphics
using AlgebraOfGraphics: density
using CairoMakie
using CategoricalArrays
using Chain
using DataFrameMacros
using DataFrames
using MixedModels
using MixedModelsMakie
using MixedModelsMakie: qqnorm
using ProgressMeter
using Random
using SMLP2023: dataset
using StatsBase
using AlgebraOfGraphics: density
import ProgressMeter
ProgressMeter.ijulia_behavior(:clear)
CairoMakie.activate!(; type="svg")
```
Expand Down Expand Up @@ -432,7 +430,7 @@ Interesting.
Here we

- generate a bootstrap sample
- compute shortest covergage intervals for the LMM parameters
- compute shortest coverage intervals for the LMM parameters
- plot densities of bootstrapped parameter estimates for residual, fixed effects, variance components, and correlation parameters

## Generate a bootstrap sample
Expand Down Expand Up @@ -487,6 +485,7 @@ The shortest coverage interval for the `GM` ranges from 376 to 404 ms and the as
#| code-fold: true
#| label: fig-betadensitym1
#| fig-cap: ' Kernel density estimate from bootstrap samples of the fixed effects for model m1L'
#=
rn = renamer([
"(Intercept)" => "GM",
"CTR: sod" => "spatial effect",
Expand All @@ -499,6 +498,7 @@ rn = renamer([
"CTR: sod, CTR: dod" => "spatial, attraction",
"CTR: dos, CTR: dod" => "object, attraction",
])
=#
draw(
data(tbl) *
mapping(
Expand Down
4 changes: 2 additions & 2 deletions kwdyz11.qmd
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
---
title: "RePsychLing Kliegl et al. (2010)"
title: "RePsychLing Kliegl et al. (2011)"
jupyter: julia-1.9
author: "Reinhold Kliegl"
---
Expand Down Expand Up @@ -284,7 +284,7 @@ We generate 2500 samples for the 15 model parameters (4 fixed effect, 4 VCs, 6 C
```{julia}
#| code-fold: true
Random.seed!(1234321)
samp = parametricbootstrap(2500, m1; hide_progress=true)
samp = parametricbootstrap(2500, m1)
tbl = samp.tbl
```

Expand Down
2 changes: 1 addition & 1 deletion sleepstudy.qmd
Original file line number Diff line number Diff line change
Expand Up @@ -216,7 +216,7 @@ Simulating and fitting a substantial number of model fits, 5000 in this case, ta

```{julia}
rng = Random.seed!(42) # initialize a random number generator
m1bstp = parametricbootstrap(rng, 5000, m1; hide_progress=true)
m1bstp = parametricbootstrap(rng, 5000, m1)
tbl = m1bstp.tbl
```

Expand Down

0 comments on commit 9062918

Please sign in to comment.