2016-07-08 5 views
1

나는 18 년 동안 시뮬레이션을하고 관찰 된 데이터를 가지고 있으며이 코드를 사용하여 RMSE를 계산할 수 있습니다. sqrt(mean((df$simulated-df$observed)^2 , na.rm = TRUE))다른 시간 척도로 RMSE를 계산하는 방법

하지만 RMSE를 다른 시간 계열의 전체 기간, 매년, 매월, 매주 및 계절별 시간 척도로 계산해야합니다. 또한, 나는 상관 계수, 분산, 바이어스 및 평균을 수행하고 하나의 파일에 모든 결과를 넣고 싶습니다.

도움이 될 것입니다.

답변

1

hydroTSMhydroGOF 패키지를 체크 아웃해야합니다. 당신은 코드를 공유하고 싶습니다,

# Getting the new numeric goodness-of-fit measures 
    gof(sim = Simulated, obs = Observed) 

    # Plot 'obs' vs 'sim' for the daily, 
    # monthly and annual time series 
    ggof(sim = Simulated, obs = Observed, ylab = "Q (ft3/day)", 
     ftype = "dma", FUN = mean) 

Daily, Monthly and Annual timescale simulation vs observation

ggof(sim = Sim, obs = Obs, ylab = "Q (ft3/day)", 
     ftype = "seasonal", FUN = mean) 

Seasonal simulation vs observation

+0

네, 사용하고 hydroGOF 패키지 : 그들은 당신이

예 찾고있는 모든 것을해야합니까? – irfan

+0

@irfan : 내 편집을 참조하십시오 – Tung

+0

완벽하게 작동합니다. – irfan

0

스케일링이 문제가된다면 미리 데이터를 표준화하고 표준화하는 것이 좋습니다. 이는 날짜를 의미있는 숫자로 인코딩하여 수행 할 수 있습니다.

적합 모델의 요약에서 원하는 대부분의 통계를 추출 할 수 있습니다.이 요약은 str()으로 검색 할 수 있습니다. Metrics 패키지에서 RMSE를 비롯한 많은 유용한 통계를 얻을 수도 있습니다.