線形回帰関数をループし、結果の勾配を新しい列として追加する

Question

11年間293サイトのデータがあり、各サイトで線形モデルを全年にわたって実行し、取得した勾配値を抽出して追加したいデータに追加します。 、と呼ばれる年

また、私は取得したいのですが、可能な場合

 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 

    1 20.8 16.5 21.2 24.6 23.2 26.4 23.8 24.9 20.3 19.5 NA 
    2 21.9 19.5 19.8 25.9 25.6 26.7 24.4 22.9 21.8 18.6 15.8 
    3 33.7 44.2 44.8 51.9 59.6 61.4 57.7 61.7 46.5 48.3 45.5 
    4 41.5 48.0 47.1 53.9 67.6 58.3 62.0 53.4 50.7 52.2 47.1 
    5 40.7 45.9 39.9 41.3 47.1 47.2 46.0 41.6 35.8 37.7 37.4 
    6 37.5 38.3 37.1 44.3 49.4 52.6 48.0 44.2 35.9 37.4 37.5 

：

は、私はすべてのすべての年のデータでは、only_meansと年のC（2015 2005）を持つベクトルでDFを持っています各モデルのR^2値のホールドとあまりに

編集することを追加します。独立変数は年、従属変数は、測定値は

の予想される出力を提供している点を通るベストフィットのラインの勾配であります

Answer 1

それはあなたのサイトごとのR二乗値を得ることはありませんが、あなたが長い形式でそれらを持っている場合lmList()は、グループによる線形回帰に適合します：

dd <- read.table(header=TRUE,check.names=FALSE,row.names=1,text=" 
    2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 

    1 20.8 16.5 21.2 24.6 23.2 26.4 23.8 24.9 20.3 19.5 NA 
    2 21.9 19.5 19.8 25.9 25.6 26.7 24.4 22.9 21.8 18.6 15.8 
    3 33.7 44.2 44.8 51.9 59.6 61.4 57.7 61.7 46.5 48.3 45.5 
    4 41.5 48.0 47.1 53.9 67.6 58.3 62.0 53.4 50.7 52.2 47.1 
    5 40.7 45.9 39.9 41.3 47.1 47.2 46.0 41.6 35.8 37.7 37.4 
    6 37.5 38.3 37.1 44.3 49.4 52.6 48.0 44.2 35.9 37.4 37.5 
") 

は長い形式に広い変換：

dd$site <- factor(rownames(dd)) 
dd_long <- tidyr::gather(dd,year,value,-site,convert=TRUE) 

フィット：

fit <- lme4::lmList(value~year|site,data=dd_long) 
## Call: lme4::lmList(formula = value ~ year | site, data = dd_long) 
## Coefficients: 
## (Intercept)  year 
## 1 -394.3945 0.2072727 
## 2 667.1091 -0.3209091 
## 3 -1743.9000 0.8927273 
## 4 -842.4727 0.4454545 
## 5 1176.6091 -0.5645455 
## 6 202.8182 -0.0800000