[R] adding regression curve to xyplot

Tue Nov 30 17:18:41 CET 2004

Inside the panel, you have to only use the parts of x and y for each unique
combination of factor(lonLabels[whichLon100])*factor(latLabels[whichLat100])

So use   "thislm <- lm(x[subscripts]~y[subscripts])" in the panel function.

Matt Austin
Statistician

Amgen 
One Amgen Center Drive
M/S 24-2-C
Thousand Oaks CA 93021
(805) 447 - 7431

> -----Original Message-----
> From: r-help-bounces at stat.math.ethz.ch
> [mailto:r-help-bounces at stat.math.ethz.ch]On Behalf Of Carlisle Thacker
> Sent: Tuesday, November 30, 2004 8:9 AM
> To: r-help at stat.math.ethz.ch
> Subject: [R] adding regression curve to xyplot
> 
> 
> Dear R-listers,
> 
> It seems that predict() behaves differently within panel.xyplot.  Am I
> doing something stupid?
> 
> Thanks,
> 
> Carlisle
> 
> First, without xyplot():
> > lmtest <- lm(t~s,data=subset(P100,whichLon100==1 & whichLat100==1))
> > lmtest
> 
> Call:
> lm(formula = s ~ t, data = subset(P100, whichLon100 == 1 & 
> whichLat100 ==
>    1))
> 
> Coefficients:
> (Intercept)            t
>     33.3307       0.1393
> > range(P100$t)
> [1]  4.050469 24.514543
> >> predict(lmtest,newdata=data.frame(t=range(P100$t)))
>        1        2
> 33.89501 36.74620
> 
> As expected, predict gives two values.  But inside xyplot() 
> predict gives
> 300 values:
> 
> > 
> xyplot(t~s|factor(lonLabels[whichLon100])*factor(latLabels[whi
> chLat100]),
> +   data=P100,pch=".",
> +   panel=function(x,y,...){panel.xyplot(x,y,...)
> +                           thislm <- lm(x~y)
> +                      print(thislm)
> +                           newt <- range(P100$t)
> +                      print(newt)
> +                           news <-
> as.vector(predict(thislm,newdata=data.frame(t=newt)))
> +                      print(news)
> +                           llines(news,newt,col="red")})
> 
> Call:
> lm(formula = x ~ y)
> 
> Coefficients:
> (Intercept)            y
>     33.3307       0.1393
> 
> [1]  4.050469 24.514543
>   [1] 34.16173 35.31284 34.95317 35.62314 34.64448 34.06225 
> 34.19688 34.03391
>   [9] 33.90860 35.28849 35.29911 34.61618 35.19837 35.28966 
> 34.91978 35.32388
>  [17] 35.39351 35.54533 35.58760 34.77039 35.07226 35.35643 
> 35.40594 33.98609
>  [25] 33.96461 33.91158 34.12389 34.33804 34.72235 34.63541 
> 36.20420 36.13635
>  [33] 34.60623 34.32401 35.59768 34.68576 34.64661 34.58804 
> 34.74548 35.09062
>  [41] 34.70290 34.22428 35.12743 34.90726 34.52325 34.90041 
> 35.46984 34.79572
>  [49] 35.47026 35.35481 34.70722 34.08883 33.97586 33.94797 
> 33.96165 33.97715
>  [57] 33.94886 33.99715 34.10887 34.72597 34.20534 34.45016 
> 35.00769 34.56742
>  [65] 34.76119 34.87190 34.60760 34.50113 34.56594 34.65954 
> 34.57076 34.91072
>  [73] 34.97322 34.99036 34.48229 34.58436 34.56767 34.94204 
> 34.95929 34.89107
>  [81] 34.98846 34.51328 34.74252 34.61337 34.62212 34.70953 
> 34.74080 34.71630
>  [89] 34.68599 35.58629 34.77031 34.49366 34.49873 34.48112 
> 35.11327 34.62598
>  [97] 34.50436 34.50705 34.70170 35.01278 35.12273 35.15541 
> 35.06295 35.15366
> [105] 35.15089 35.11312 35.09834 35.05271 34.84838 34.52225 
> 34.34650 34.32637
> [113] 34.42911 34.79282 35.06830 35.11559 35.25736 34.43281 
> 34.45208 34.58407
> [121] 35.28598 34.34246 34.32532 34.47483 34.79309 34.02772 
> 34.01553 35.33372
> [129] 34.19273 34.07622 34.52820 35.27539 35.37703 34.70560 
> 34.77501 34.52084
> [137] 34.68773 35.19217 35.69538 34.21886 34.17226 34.94784 
> 34.03079 33.95669
> [145] 33.92002 33.92972 33.96568 33.92996 33.94236 33.93887 
> 33.92059 33.91746
> [153] 33.98369 34.17190 34.05327 33.96500 34.16136 33.98172 
> 34.02463 34.00802
> [161] 34.05238 34.05088 34.16828 34.14888 33.94711 33.96567 
> 34.12450 33.98988
> [169] 33.98670 33.98489 34.05447 34.03689 34.17137 34.31777 
> 34.05027 33.97821
> [177] 33.99281 34.00209 34.07560 33.98840 33.98930 33.98548 
> 33.99127 35.23652
> [185] 35.17034 35.53546 35.53861 34.02384 34.03402 33.93877 
> 33.93540 33.90636
> [193] 33.92267 34.57678 34.49457 34.46038 34.64318 35.06561 
> 35.06680 34.76466
> [201] 34.00418 35.48020 35.12687 35.28180 35.07281 35.23499 
> 35.14693 34.07914
> [209] 33.99047 34.00648 34.10073 34.00784 33.99185 33.98191 
> 34.03089 33.96755
> [217] 34.17042 34.10803 34.12671 34.09942 34.25003 34.06542 
> 34.07293 34.19048
> [225] 34.06594 34.01929 34.00155 34.11555 33.89754 33.89527 
> 33.89501 33.93433
> [233] 33.89881 33.94577 33.96603 33.95199 34.00284 34.12974 
> 33.98116 34.08882
> [241] 34.00953 34.98977 34.72824 34.68920 35.15207 34.59716 
> 34.64580 34.25087
> [249] 34.96008 34.62607 35.32346 35.16993 34.46759 34.47157 
> 34.46978 34.46092
> [257] 34.46561 34.46489 34.49970 34.49599 34.55441 34.52048 
> 34.58480 34.59267
> [265] 34.50082 34.57311 34.60421 34.57043 34.53131 34.46365 
> 34.41946 34.74566
> [273] 34.64593 34.67921 34.65706 34.71625 34.61487 34.61687 
> 34.58899 34.47081
> [281] 34.58366 35.66823 35.61883 34.71477 34.72518 34.69732 
> 34.66265 34.65186
> [289] 34.81468 34.74725 34.77658 34.77191 34.81481 34.76427 
> 34.70751 34.96763
> [297] 35.11322 35.16852 34.94393 34.89605
> Error in xy.coords(x, y) : x and y lengths differ
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! 
http://www.R-project.org/posting-guide.html