[R] regression by group summary error

[arun]

Hi Catalin,

I tried with a subset of the variables.  Infact, there is an option in lmList() to subset

biN <- bi[,c(1,3,22,34)]
str(biN)
'data.frame':    66 obs. of  4 variables:
 $ Exp          : chr  "B" "B" "B" "B" ...
 $ Clona        : Factor w/ 5 levels "A4A","AF2","Max4",..: 3 3 3 3 3 3 3 3 3 3 ...
 $ CR2V.cm.     : num  20.4 19.6 24 22.2 27.6 24.2 19.4 22.2 21.5 28.2 ...
 $ masa.uscat.tr: num  6.2 3.88 7.82 5.65 10.2 ...


 by(biN,biN$Exp, function(x) summary(lmList(masa.uscat.tr~CR2V.cm.|Clona,data=x,na.action=na.omit)))
biN$Exp: B
Call:
  Model: masa.uscat.tr ~ CR2V.cm. | Clona 
   Data: x 

Coefficients:
   (Intercept) 
    Estimate   Std. Error      t value     Pr(>|t|) 
-8.252204377  1.668647582 -4.945444721  0.001127319 
   CR2V.cm. 
    Estimate   Std. Error      t value     Pr(>|t|) 
6.674136e-01 7.218748e-02 9.245559e+00 1.519752e-05 

Residual standard error: 0.6671037 on 8 degrees of freedom

------------------------------------------------------------ 
biN$Exp: C
Call:
  Model: masa.uscat.tr ~ CR2V.cm. | Clona 
   Data: x 

Coefficients:
   (Intercept) 
           Estimate Std. Error   t value   Pr(>|t|)
A4A       -6.941742   6.698275 -1.036348 0.30545738
AF2      -10.368170   4.669365 -2.220467 0.03135312
Max4      -6.766456   2.925135 -2.313211 0.02523457
Monviso    9.493789   3.779740  2.511757 0.01558630
Pannonia  -8.044174   3.617619 -2.223610 0.03112594
   CR2V.cm. 
           Estimate Std. Error   t value     Pr(>|t|)
A4A      0.62752742  0.2513938 2.4961928 1.619954e-02
AF2      0.72802937  0.1609617 4.5229977 4.270833e-05
Max4     0.63517380  0.1377759 4.6101954 3.209599e-05
Monviso  0.04123701  0.1366451 0.3017818 7.641786e-01
Pannonia 0.74293255  0.1187515 6.2561931 1.194664e-07

Residual standard error: 1.473519 on 46 degrees of freedom




A.K.






On Wednesday, November 20, 2013 5:53 AM, catalin roibu <catalinroibu at gmail.com> wrote:

Hello!
I have a problem with R. I want to apply a linear model for data composed by two groups (Exp and Clona). i try this code:
by(bi,bi$Exp, function(x) lmList(masa.uscat.tr~CR2V.cm.|Clona,data=x,na.action=na.omit))

But I want to view to summary for each Exp and i try this code:
by(bi,bi$Exp, function(x) summary(lmList(masa.uscat.tr~CR2V.cm.|Clona,data=x,na.action=na.omit)))

And I have this error:
Error in `[<-`(`*tmp*`, use, use, ii, value = lst[[ii]]) : 
  subscript out of bounds

Please help me to solve this problem!

Thank you very much!

 My data is:

ExpPlotClonaProvmasa uscat trmasa usc. Rammasa usc total
B36Max4 PPuieti6.1998484852.6393258438.839174328
B36Max4 PPuieti3.878751.47985.35855
B36Max4 PPuieti7.8227027033.3285207111.15122341
B36Max4 PPuieti5.6453846151.9952380957.640622711
B36Max4 PPuieti10.23.51486486513.71486486
B36Max4 PPuieti8.8155454553.3562790712.17182452
B36Max4 PPuieti5.0333333331.6071428576.64047619
B36Max4 PPuieti6.6934883722.6302083339.323696705
B36Max4 PPuieti6.0210126581.602935787.623948438
B36Max4 PPuieti10.205825243.76831460713.97413985
B33Max4Butasi5.8994174761.7453947377.644812213
B33Max4Butasi3.2614285711.3357352944.597163866
B33Max4Butasi3.3595081971.4566412214.816149418
B33Max4Butasi5.0363636362.0977931037.13415674
B33Max4Butasi3.1221621621.6120125794.734174741
B33Max4Butasi5.0424742273.2469164278.289390653
B33Max4Butasi5.0582558142.7242990657.782554879
B33Max4Butasi4.9778181821.7135042746.691322455
B33Max4Butasi3.1952941181.2434117654.438705882
B33Max4Butasi1.8318181821.0090909092.840909091
B30AF2Butasi6.981951221.7648.74595122
B30AF2Butasi5.8583333331.6866233777.54495671
B30AF2Butasi10.6253.0413.665
C39AF6Sade10.271252.28319327712.55444328
C39AF6Sade9.4734883721.90941441411.38290279
C39AF6Sade10.4825212.4825
C39AF6Sade11.615795452.43113636414.04693182
C39AF6Sade8.1850746271.8093333339.99440796
C39AF6Sade11.045106382.2367295613.28183594
C39AF6Sade9.0666666672.47378571411.54045238
C39AF6Sade10.127876113.09763157913.22550769
C39AF6Sade9.1712903231.82122641510.99251674
C39AF6Sade12.18468752.26259036114.44727786
C42PannoniaSade9.2752.48217391311.75717391
C42PannoniaSade7.211.778.98
C42PannoniaSade11.369393943.11178082214.48117476
C42PannoniaSade7.852968751.9434751779.796443927
C42PannoniaSade8.252.5404761910.79047619
C42PannoniaSade8.6692771082.18707142910.85634854
C42PannoniaSade8.5108860762.0534410.56432608
C42PannoniaSade9.3622222225.52553191514.88775414
C42PannoniaSade11.084819282.57319327713.65801255
C42PannoniaSade10.174626873.00322580613.17785267
C45MonvisoSade12.996938783.21608391616.21302269
C45MonvisoSade11.114565221.88571428613.0002795
C45MonvisoSade8.1293333331.5326666679.662
C45MonvisoSade9.9430434782.3830088512.32605233
C45MonvisoSade11.98058143.08092391315.06150531
C45MonvisoSade10.313766232.21052631612.52429255
C45MonvisoSade9.9475862071.86083333311.80841954
C45MonvisoSade12.242615382.16685714314.40947253
C45MonvisoSade13.566506022.41437125715.98087728
C45MonvisoSade10.875742572.92234042613.798083
C48AF2Sade9.3345454552.20182291711.53636837
C48AF2Sade9.7476404492.81178082212.55942127
C48AF2Sade14.295412844.50688524618.80229809
C48AF2Sade10.264516133.01432258113.27883871
C48AF2Sade13.309242423.96066176517.26990419
C48AF2Sade13.232284265.0054687518.23775301
C48AF2Sade14.82770274.45044776119.27815046
C48AF2Sade15.236695287.55932203422.79601731
C48AF2Sade13.271985823.75825242717.03023824
C48AF2Sade13.713644445.69044444419.40408889
C51Max4Sade8.8608849565.61358695714.47447191
C51Max4Sade13.291532855.65306122418.94459407
C51Max4Sade9.6098507463.38571428612.99556503


-- 

---
Catalin-Constantin ROIBU
Lecturer PhD, Forestry engineer
Forestry Faculty of Suceava
Str. Universitatii no. 13, Suceava, 720229, Romania
office phone       +4 0230 52 29 78, ext. 531
mobile phone   +4 0745 53 18 01
                       +4 0766 71 76 58
FAX:                          +4 0230 52 16 64
silvic.usv.ro