«lme4-nlme» 태그된 질문

lme4 및 nlme는 선형, 일반화 선형 및 비선형 혼합 효과 모델을 피팅하는 데 사용되는 R 패키지입니다. 혼합 모델에 대한 일반적인 질문은 [mixed-model] 태그를 사용하십시오.

4
glmer의 랜덤 효과 분산 해석
데이터가 이항 분포 (과일 성숙 또는 성숙하지 않음) 인 수분에 관한 논문을 개정하고 있습니다. 그래서 나는 glmer하나의 무작위 효과 (개별 식물)와 하나의 고정 효과 (치료)를 사용했습니다. 리뷰어는 식물이 과일 세트에 영향을 미쳤는지 알고 싶어하지만 glmer결과를 해석하는 데 문제가 있습니다. 나는 웹을 읽었으며 직접 비교 glm하고 glmer모델에 문제가있을 수 있으므로 그렇게하지 …

2
Nakagawa & Schielzeth (2013) R2glmm 방법을 사용하여 혼합 모형에서
혼합 모델에서 값 계산에 대해 읽었으며 R-sig FAQ,이 포럼의 다른 게시물 (몇몇 링크하지만 평판이 충분하지 않음) 및 R 을 사용하는 것으로 이해하는 몇 가지 다른 참조를 읽었습니다. 혼합 모델의 맥락에서 2 개의 값은 복잡합니다.아르 자형2R2R^2아르 자형2R2R^2 그러나 최근에이 두 논문을 접했습니다. 이 방법들이 유망한 것처럼 보이지만 (나는) 통계학자가 아니므로 다른 …

1
비선형 혼합 모형 (nlme) 예측에 대한 신뢰 구간
비선형 혼합 nlme모형 의 예측에 대한 95 % 신뢰 구간을 얻고 싶습니다 . 아무것도 표준 내에서이 작업을 수행하기 위해 제공되기 때문에 nlme, 나는 "인구 예측 간격"의 방법을 사용하는 것이 올바른 것처럼, 궁금 모델의 맥락에서 벤 Bolker의 책 장에 설명 된 최대 우도에 맞게 의 아이디어를 기반으로, 적합 모형의 분산 공분산 …

1
기본 lme4 옵티마이 저는 고차원 데이터에 대해 많은 반복이 필요합니다
TL은, DR은 : lme4최적화 기본적 모델 파라미터의 수에 선형으로 나타날 것이다 방법 등가보다 느리게 glm그룹 더미 변수 모델. 속도를 높이기 위해 할 수있는 일이 있습니까? 상당히 큰 계층 적 로짓 모델 (~ 50k 행, 100 열, 50 그룹)을 맞추려고합니다. 일반 로짓 모델을 데이터에 그룹화하면 (그룹에 더미 변수가 있음) 제대로 작동하지만 …

2
SAS PROC GLIMMIX가 이항 glmm에 대해 glmer (lme4)와 다른 임의의 기울기를 제공하는 이유는 무엇입니까?
저는 R에 더 익숙한 사용자이며 4 개의 서식지 변수에 대해 5 년 동안 약 35 명의 개인에 대한 임의의 기울기 (선택 계수)를 추정하려고했습니다. 응답 변수는 위치가 "사용 된"(1) 또는 "사용 가능한"(0) 서식지 (아래 "사용") 여부입니다. Windows 64 비트 컴퓨터를 사용하고 있습니다. R 버전 3.1.0에서는 아래 데이터와 표현식을 사용합니다. PS, TH, …

1
Lmer 모델이 수렴하지 않습니다
내 데이터는 여기에 설명되어 있습니다 . 반복 측정 ANOVA를 피팅 할 때 "ov (Error () model is singleular error)"가 발생하는 원인은 무엇입니까? lmer내 기본 사례는 다음을 사용하여 상호 작용의 효과를 보려고합니다 . my_null.model <- lmer(value ~ Condition+Scenario+ (1|Player)+(1|Trial), data = my, REML=FALSE) my.model <- lmer(value ~ Condition*Scenario+ (1|Player)+(1|Trial), data = …
12 r  lme4-nlme 

1
ZINB 또는 다른 효과가 혼합 된 카운트 데이터에 적합한 모델을 찾는 데 어려움이 있습니까?
나는 고독한 꿀벌 풍부도에 대해 매우 작은 데이터 세트를 가지고 있는데, 분석하는데 문제가 있습니다. 그것은 계수 데이터이며, 거의 모든 계수는 다른 처리에서 대부분의 0으로 처리됩니다. 6 개의 사이트 중 2 개에 각각 하나씩 두 개의 매우 높은 값이 있으므로 카운트 분포가 매우 긴 꼬리를 갖습니다. R에서 일하고 있습니다. lme4와 glmmADMB의 …

1
lmer와 p- 값에 대한 혼란 : memisc 패키지의 p- 값과 MCMC의 p- 값은 어떻게 비교됩니까?
패키지 의 함수 lmer()가 lme4p- 값을 생성하지 않았다는 인상을 받았습니다 ( lmerp- 값 및 그 밖의 모든 것을 참조하십시오 ). 이 질문에 따라 대신 MCMC에서 생성 된 p 값을 사용 했습니다. lme4혼합 모델의 중요한 효과 및이 질문 : 의 패키지에서 출력에서 p- 값을 찾을 수 없습니다lmer()lm4R . 최근에 memisc 라는 …

2
레벨 당 1 개의 관측치가있는 혼합 모형
glmer비즈니스 데이터에 임의 효과 모델을 적용하고 있습니다. 목표는 지역별 변동을 고려하여 총판 별 판매 실적을 분석하는 것입니다. 다음과 같은 변수가 있습니다. distcode: 약 800 레벨의 총판 ID region: 최상위 지리 ID (북쪽, 남쪽, 동쪽, 서쪽) zone: 중간 수준 지리 내에 region약 30 개 수준이 중첩되어 있습니다. territory: zone약 150 단계 …

1
히든 마르코프 모델에서 "최상의"모델을 선택하기위한 기준
데이터의 잠재 상태 수를 추정하기 위해 HMM (Hidden Markov Model)에 맞추려고하는 시계열 데이터 세트가 있습니다. 이 작업을 수행하는 의사 코드는 다음과 같습니다. for( i in 2 : max_number_of_states ){ ... calculate HMM with i states ... optimal_number_of_states = "model with smallest BIC" ... } 이제 일반적인 회귀 모델에서 BIC는 가장 …

2
R에서 lmer ()를 사용하여 Poisson GLMM에서과 분산을 테스트하는 방법은 무엇입니까?
다음과 같은 모델이 있습니다. > model1<-lmer(aph.remain~sMFS1+sAG1+sSHDI1+sbare+season+crop +(1|landscape),family=poisson) ... 이것은 요약 출력입니다. > summary(model1) Generalized linear mixed model fit by the Laplace approximation Formula: aph.remain ~ sMFS1 + sAG1 + sSHDI1 + sbare + season + crop + (1 | landscape) AIC BIC logLik deviance 4057 4088 -2019 4039 Random effects: …

1
nlmer ()를 사용하여 반복 측정 데이터에 비선형 혼합 효과 모델을 맞추려면 어떻게해야합니까?
반복 측정 데이터를 분석하려고 노력하고 R있습니다. 내 데이터는 본질적으로 다음과 같습니다. 두 개의 치료 그룹이 있습니다. 각 그룹의 모든 과목은 매일 시험을 치르며 점수 (시험에서 정확한 비율)를받습니다. 데이터는 긴 형식입니다. Time Percent Subject Group 1 0 GK11 Ethanol 2 0 GK11 Ethanol 3 0 GK11 Ethanol 4 0 GK11 Ethanol …

2
lmer로부터 자유도 확보
나는 출력을 구성했지만 다음과 같은 lmer 모델을 적합시켰다. Random effects: Groups Name Std.Dev. day:sample (Intercept) 0.09 sample (Intercept) 0.42 Residual 0.023 다음 공식을 사용하여 각 효과에 대한 신뢰 구간을 만들고 싶습니다. ( n - 1 ) s2χ2α / 2 , n - 1, ( n - 1 ) s2χ21 − …

1
전체 성공을 가진 범주 형 변수가있는 이항 glmm
이항 반응 변수와 범주 형 예측 변수로 glmm을 실행하고 있습니다. 무작위 효과는 데이터 수집에 사용 된 중첩 디자인에 의해 제공됩니다. 데이터는 다음과 같습니다. m.gen1$treatment [1] sucrose control protein control no_injection ..... Levels: no_injection control sucrose protein m.gen1$emergence [1] 1 0 0 1 0 1 1 1 1 1 1 0 …

2
3 방향 반복 측정 분산 분석에 해당하는 lme4 :: lmer는 무엇입니까?
내 질문은 이lme4::lmer 모델이 양방향 반복 측정 ANOVA에 해당 하는 모델 을 보여준 이 응답 을 기반으로 합니다 . require(lme4) set.seed(1234) d <- data.frame( y = rnorm(96), subject = factor(rep(1:12, 4)), a = factor(rep(1:2, each=24)), b = factor(rep(rep(1:2, each=12))), c = factor(rep(rep(1:2, each=48)))) # standard two-way repeated measures ANOVA: summary(aov(y~a*b+Error(subject/(a*b)), …

당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.