Message-ID: <BN7PR08MB5252B87BC91B76072BB9F040FF040@BN7PR08MB5252.namprd08.prod.outlook.com>
Date: 2020-10-13T20:56:42Z
From: Sammie Haskin
Subject: Assessing whether sigma for a random effects parameter is equal to 0
Hello! Given the pistonrings data set from the qcc package in R, I produced the following code to assess whether the standard deviation of random effect of the model was equal to 0. Here is my code:
library(lme4)
library(qcc)
library(RLRsim)
library(nlme)
library(data.table)
fit.pistons <- lmer(formula=diameter ~ sample + (1 | sample), data = pistonrings,REML=T)
fit.pistons0 <- lm(diameter ~ sample, data = pistonrings)
exactLRT(fit.pistons,fit.pistons0)
Here is the output:
LRT = 8.1423, p-value = 0.0014
Is this result implying that the standard deviation for the random effect is significant such that we reject the null hypothesis and that H0: sigma = 0 is false?
[[alternative HTML version deleted]]