nlme

    0

    1답변

    nlme 패키지를 사용하여 반복 측정 (MMRM) 모델을 R로 혼합하려고합니다. 데이터 구조는 다음과 같습니다. 각 환자는 세 그룹 (grp) 중 하나에 속하며 치료 그룹 (trt)에 할당됩니다. 환자 결과 (y)는 6 회의 방문 (방문) 동안 측정됩니다. 다른 방문에 대해 이질적 분산이있는 복합 대칭 모델을 사용하고 싶습니다 (SAS의 PROC MIXED

    0

    1답변

    저는 R 및 혼합 모델 분석에 상당히 익숙합니다. 나는 모든 개인에 대한 time에 따라 변수 ln_ahr의 변화에 ​​대해 하나의 추정을 생성 할 수 있도록 싶었다. 나는 이것이 시간의 변화의 기울기로 생각할 수 있다고 생각한다. 여기 내 데이터 (긴 형식)을 구성하는 방법입니다 : v001 주제 식별자 v001 ln_ahr time 13404 283

    0

    1답변

    library(nlme) fm1 <- nlme(height ~ SSasymp(age, Asym, R0, lrc), data = Loblolly, fixed = Asym + R0 + lrc ~ 1, random = Asym ~ 1, start = c(Asym = 103, R0 = -8.5, lrc = -3.

    0

    1답변

    제목에서 lme4 (또는 다른 패키지?) 객체에서 행렬을 추출하려고합니다. 이 행렬은 G라고하고 있습니다 SAS 표기법에서 임의 효과 의 https://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_mixed_sect022.htm 분산 - 공분산 행렬 : 나는

    0

    1답변

    read.affybatch() 명령을 사용하여 열려고하는 Affymetrix 데이터가 affy에 있습니다. 나는 문제없이 함께 biocLite 및 affy 패키지를 다운로드했지만 한 최대한 빨리는 다음과 같은 문제가 나타나면 읽기 명령을 사용하여 내 .CEL 파일을 시도하고 읽을 때 : Library - package hgu133acdf not instal

    2

    1답변

    lme 명령을 사용하여 상호 작용 무작위 효과가 있지만 무작위 효과가없는 LMM에 맞춰야합니다. 즉, 모델을 oats.lmer (아래 참조)에 맞추고 싶지만 nlme 패키지의 lme 함수를 사용하고 싶습니다. 코드는 내가이 oats.lme <- lme(yield~nitro, data = Oats, random = (~1|Block/Variety)) sum

    0

    1답변

    패키지에서 하나의 계승 결과 lme을보고하고 싶습니다. 나는 A에 대한 전반적인 효과를 알고 싶다. 이렇게하려면 내가 널 모델로 모델을 비교하는 것이다 : 나는 다른 주요 효과 모델을 비교하고 있기 때문에 최대 가능성을 사용하고 m1 <- lme(y~A,random=~1|B/C,data=data,weights=varIdent(form = ~1|A),meth

    0

    1답변

    R 패키지 nlme에서 gls 함수를 사용하여 중첩 된 공간 샘플 집합을 분석하고자합니다.이 예제에서는 많은 샘플이 적어도 일부 공간 좌표에서 중첩됩니다. 나는 corStruct 또는 pdMat 객체를 사용하여 응답 변수 (각 공간 샘플에서 측정하는 것)에서 비 독립성을 설명하려고하지만이 작업을 수행하는 방법에 대해 혼란스러워합니다. 공간 샘플 간의 비 독

    0

    1답변

    나는 다중 레벨 모델을 배우기 위해 nlme 패키지를 사용하고, 일어 났을 때 "Discovering Statistics Using R"교과서의 예제를 사용하고있다. Mixed Models Code 데이터 세트 또한 그들의 동반자 웹 사이트에서 다운로드 할 수있는, 허니문 Period.dat입니다. 나는 "타이머"모델을 업데이트하려고 할 때 require(

    0

    1답변

    세로 데이터 분석에서 숙제를 끝내려고했습니다. 질문은 모델 내에서 연령의 단면 및 종단 효과의 차이 (기본 단면 : 기본, 길이 방향 연령 : 연령 변경)를 비교하는 것입니다. 같은 모델 I 코드 : 는 < -lme에 맞게 (logfev1 ~ baseage + agechange + 높이, 임의 = ~ 1 | 방문 ~ = ID, 상관 관계 = corAR1 (