simr

simr copied to clipboard

Consider a model with random effects `(1 | A) + (1 | B) + (1 | C)`. How to test for significance of effect A? Currently might use `rcompare(~ (1...

pitakakariki

powerCurve not wrapping user supplied tests

7

This works fine: ``` > powerSim(fm, test=function(.) 1, nsim=5) ``` But this doesn't: ``` > powerCurve(fm, test=function(.) 1, nsim=5) Calculating power at 8 sample sizes along x Error in attr(test,...

pitakakariki

'long vectors not supported yet' error

7

I'm trying to use simr with a dataset and model from my experiments and it's generating `long vectors not supported yet: subassign.c:1818` errors for every attempted model fit. Here is...

Shotgunosine

Create powerggplot.R

1

Hello! I'd just like to suggest this ggplot function. If you find it useful, please feel free to edit it to taste.

pablobernabeu

Linked to #228

pablobernabeu

Mislabelled default LR test

2

``` > simdata$g fm doTest(fm, fixed("g")) boundary (singular) fit: see ?isSingular p-value for predictor 'g': 1.434264e-13 -------------------- Test: Kenward Roger (package pbkrtest) > doTest(fm, fixed("g", "lr")) boundary (singular) fit: see...

pitakakariki

``` makeLmer(y ~ x + (1|g), fixef=1:5, VarCorr=1, sigma=1, data=simdata) ``` ``` Linear mixed model fit by REML ['lmerMod'] Formula: y ~ x + (1 | g) Data: simdata REML...

pitakakariki

Taking too much time to run the code

6

Hi, I've extended the model along participant number and used powerCurve with nsim=100. It's been taking me more than 72 hrs, and the code is still running. I've tried the...

sangheek16

Play nice with lmerTest

3

Defaults to "unknown test" in `fixed`.

pitakakariki

`powerCurve` output does not indicate errors, which makes debugging scripts using `simr` difficult - 100% errors are displayed as 0.00 power. ``` r library(simr) ``` (Make some missing data to...

mattansb