Question

我有一个带有观察变量的数据框和一个运行数千行的日期戳[Var1，DD，MM，YYYY]。我需要为每年的观察变量拟合分布[指数或伽马]，并获得每年的相关参数。

在Matlab中，它只是

 j=1
 k=1

 for i=1:(no_of_rows-1)

    if  year(i+1) = year(i)
        temp_data_year(j) = Var1(i)
        j=j+1

    else  [a,b]= gamfit(temp)
         param(:,:,k) = [a,b]
         k=k+1
   endif

 end

因此，我将在变量参数的数据中获取每年的参数。

R中有什么可以做到的吗？

谢谢，

Answer 1

喜欢这个。

# creates a sample dataset - you have this already
set.seed(1)             # for reproducible example
df <- data.frame(var1=c(rgamma(365,2,4),rgamma(365,3,5),rgamma(365,1,8)),
                 YYYY=rep(2012:2014,each=365))

# you start here...
library(fitdistrplus)   # for fitdist(...)
aggregate(var1~YYYY,df,function(X)fitdist(X,distr="gamma")$estimate)
#   YYYY var1.shape var1.rate
# 1 2012   1.891706  3.873906
# 2 2013   2.812962  4.778191
# 3 2014   1.031067  7.826776

阅读fitdist(...)包中fitdistrplus的文档。有几种拟合算法可用。

每年循环并申请功能

1 个答案: