CAFE

댓글

뒤로

[회귀분석]가중최소제곱,신뢰구간,예측구간..

작성자나스키즈| 작성시간12.07.20| 조회수660| 댓글 8

댓글 리스트

  • 작성자 안재형 작성시간12.07.21 얼핏보기에 out=lm(y~x,weights=x) 에서 wegiths=x때문이 아닌가 합니다.
    weights=x를 빼고 다시 해보신후 그래도 직선이 안나오면 알려주세요.
  • 답댓글 작성자 나스키즈 작성자 본인 여부 작성자 작성시간12.07.31 답변 감사드립니다__) 말씀해주신대로 ")"를 지우고
    lim.p=as.data.frame(predict(out,x_axis,level=0.95,interval="predict",weights=x))
    로 수정하여 다시 수행했는데, 저는 결과가 이전과 다름이 없네요... ㅠㅠ 혹시 제공 패키지 나 다른 툴이 필요한 것일까요? 확인부탁드립니다 선생님 ㅠㅠ
  • 답댓글 작성자 나스키즈 작성자 본인 여부 작성자 작성시간12.07.31 해결했습니다. 잘 나오네요!! 감사합니다!!
  • 작성자 나스키즈 작성자 본인 여부 작성자 작성시간12.07.23 네.. weights=x를 빼면 직선으로 나옵니다. 위에 제가 추가적으로 올린 것 처럼, 예측범위를 x축에 따라 곡선의 형태로 표현할 수는 없는 것인가요? 상기 그림은 다른 논문에서 본 것입니다;;
  • 답댓글 작성자 안재형 작성시간12.07.26 이제야 봤습니다. 제 댓글에 답글을 쓰시면 제가 알수 있는데 그냥 여기에 올리면 알수가 없습니다.
    전에는 안돌리다가 이제 돌려봤는데 저는 멀쩡하게 부드러운 그래프가 나오는데요. 아래 명령어에서 맨끝에 ")"가 3개인데 한개 지우세요.
    lim.p=as.data.frame(predict(out,x_axis,level=0.95,interval="predict",weights=)))
    거기서 에러가 나오고 전에 저장되어 있던 lim.p가 대신 작동한게 아닌가 합니다.
  • 작성자 sl100104 작성시간12.07.24 이걸 원하시나요..

    install.packages("HH")
    library(HH)

    lrfd=read.csv("LRFD.csv")
    attach(lrfd)

    out=lm(y~x)
    ci.plot(out)
  • 작성자 나스키즈 작성자 본인 여부 작성자 작성시간12.07.24 네! 맞습니다. 다만, 회귀 분석을 할 때 가중최소제곱을 이용하여 수행하였고, 더불어 신뢰구간 및 예측구간을 표시 할 때 y 값을 분산을 반영하여 그리고 싶습니다. 제가 추가적으로 올린 그래프, "통계학적 분석을 이용한 예측 극한하중의 평가" 에서 95% 예측구간을 x 축을 따라 동일한 분산을 이용해 예측을 한 것이 아니라, 분산을 달리하여 예측한 것이라 판단됩니다. 정확한 설명인 것인지는 모르겠으나, 제가 이해하기로는 그렇습니다... 그렇기에 x축의 값이 0.2 이하인 구간에서는 y값의 분산이 크기 때문에 예측값의 범위도 넓어 진 것이라 판단됩니다. 어떠한 방식으로 저러한 예측범위를 도출할 수 있을까요?
  • 답댓글 작성자 안재형 작성시간12.07.26 y 예측값의 분산은 왠만한 회귀분석책을 보면 나옵니다. 예측값의 분산은 추정치의 분산보다 큽니다. 그래서 신뢰구간도 넓어집니다.
맨위로

카페 검색

카페 검색어 입력폼
카카오 이모티콘
사용할 수 있는 카페앱에서
댓글을 작성하시겠습니까?
이동시 작성중인 내용은 유지되지 않습니다.