놀랍게도 Google을 사용하여 다음 질문에 대한 답변을 찾을 수 없었습니다.
나는 시간에 대략 시그 모이 드 성장 행동을 보여주는 여러 개인의 생물학적 데이터를 가지고 있습니다. 따라서 표준 물류 성장을 사용하여 모델링하고 싶습니다.
P(t) = k*p0*exp(r*t) / (k+p0*(exp(r*t)-1))
p = 0은 t = 0에서의 시작 값이고, k는 t-> 무한에서의 점근 한 한계이고 r은 성장 속도이다. 내가 볼 수있는 한, 나는 nls를 사용하여 이것을 쉽게 모델링 할 수 있습니다 (내 이해가 부족합니다 : 왜 시간과 데이터를 스케일링하여 표준 로짓 회귀를 사용하여 비슷한 것을 모델링 할 수 없습니까? 비율, 그러나 드물게 http://www.stata-journal.com/article.html?article=st0147 이 탄젠트에 대한 다음 질문은 모델이 이상치> 1)를 처리 할 수 있는지 여부입니다.
이제 세 가지 매개 변수 k, p0 및 r에 고정 (주로 범주 형) 및 임의의 (개인 ID 및 연구 ID) 효과를 허용하고 싶습니다. nlme가 이것을 수행하는 가장 좋은 방법입니까? SSlogis 모델이 내가하려고하는 일에 합리적 인 것 같습니다. 맞습니까? 다음 중 적합한 모델이 있습니까? 시작 값을 올바르게 얻을 수 없으며 update ()는 고정 효과가 아닌 임의의 효과에 대해서만 작동하는 것 같습니다-힌트?
nlme(y ~ k*p0*exp(r*t) / (k+p0*(exp(r*t)-1)), ## not working at all (bad numerical properties?)
data = data,
fixed = k + p0 + r ~ var1 + var2,
random = k + p0 + r ~ 1|UID,
start = c(p0=1, k=100, r=1))
nlme(y ~ SSlogis(t, Asym, xmid, scal), ## not working, as start= is inappropriate
data = data,
fixed = Asym + xmid + scal ~ var1 + var2, ## works fine with ~ 1
random = Asym + xmid + scal ~ 1|UID,
start = getInitial(y ~ SSlogis(Dauer, Asym, xmid, scal), data = data))
비선형 혼합 모델, 특히 비선형 모델에 익숙하지 않기 때문에 초보자 질문이있는 자습서 / FAQ 링크 또는 권장 사항을 읽어보십시오.