MGCV의 특정 값 추정

Question

mgcv 패키지의 권장 값인 this thread을 기반으로 주어진 값에서 GAM 모델을 실행하고 정확한 예상치를 얻으려고합니다. 예측 함수를 사용할 때 결과는 GAM의 예상과 일치하지 않습니다. 내가 사용 결과를 얻을 때

mod <- gam(y~s(a, by=b), data = dat) 

:

plot(mod) 

결과는 다릅니다 내가 변수 a에 하나의 플롯 점을 추정 할 때

는 여기가에 관심이 내 모델의 다음 방법을 사용하십시오.

pdat <- with(dat, 
      data.frame(a = 0, 
         b = mean(b,na.rm = TRUE))) 
pred <- predict(mod, pdat, type = "response", se.fit = TRUE) 
pdat <- transform(pdat, fitted = pred$fit) 
pdat <- transform(pdat, upper = fitted + (1.96 * pred$se.fit), 
        lower = fitted - (1.96 * pred$se.fit)) 

누구에게 어떤 제안이 있습니까? 내가 잘못하고있는 것일까?

Answer 1

시행 착오를 통해 내 질문에 대한 답을 찾았습니다. 변수에 의한 my 값 때문에 b 값을 1로 바꾸어야하기 때문에 값의 평균이 아니라 이음새가 있습니다.이 값은 그래프의 정확한 값을 얻습니다.

아마도 다른 사람에게 도움이 될 것입니다.