[R-sig-Geo] Plotting more than one regression line in ggplot

r@i@1290 m@iii@g oii @im@com r@i@1290 m@iii@g oii @im@com
Wed Jun 5 16:21:57 CEST 2019


I am trying to plot, using ggplot, a series of scatter plots with regression lines for several datasets. I started with the following dataset, "onepectCO2MEDIAN". The data for this dataset is as follows:
    onepctCO2MEDIAN
                    x          y
    layer.1   0.000000000  0.0000000
    layer.2   0.006794447  4.9002490
    layer.3   0.014288058  0.1608000
    layer.4   0.022087920  6.6349133
    layer.5   0.030797357 -1.2429506
    layer.6   0.038451072  1.5643374
    layer.7   0.048087904 -2.2659035
    layer.8   0.058677729  2.2070045
    layer.9   0.069261406 -2.3677001
    layer.10  0.080524530 -1.0913506
    layer.11  0.092760246  0.4099940
    layer.12  0.103789609 -0.1259727
    layer.13  0.116953168 -2.4138253
    layer.14  0.129253298  7.0890257
    layer.15  0.141710050 -0.7593539
    layer.16  0.156002052  0.0454416
    layer.17  0.170648172 -1.5349683
    layer.18  0.185318425  6.5524201
    layer.19  0.199463055 -0.8312563
    layer.20  0.213513337 -2.5099183
    layer.21  0.228839271  0.1365968
    layer.22  0.246981293 -1.3719845
    layer.23  0.263012767 -0.8712988
    layer.24  0.278505564  0.6632584
    layer.25  0.293658361  0.7938036
    layer.26  0.310747266  3.4880637
    layer.27  0.325990349 -4.4612208
    layer.28  0.342517540  0.0871734
    layer.29  0.362751633 -1.4171578
    layer.30  0.380199537 -0.9956508
    layer.31  0.394992948  0.3215526
    layer.32  0.414373398  3.1403866
    layer.33  0.430690214 -0.7376099
    layer.34  0.449738145 -2.4860541
    layer.35  0.470167458 -3.4235858
    layer.36  0.489019871  0.4824748
    layer.37  0.507242471 -0.9785386
    layer.38  0.524314284  8.5359684
    layer.39  0.543750525  5.4844742
    layer.40  0.564234197  3.2149367
    layer.41  0.583679616  3.9168916 
    layer.42  0.601459444  4.4907020
    layer.43  0.619924664  6.5410410
    layer.44  0.639932007  4.8068650
    layer.45  0.661347181  8.1510170
    layer.46  0.684117317  0.2697413
    layer.47  0.704829752 -0.1807501
    layer.48  0.725045770  9.7181249
    layer.49  0.745165825  1.5406466
    layer.50  0.765016139 -1.6476041
    layer.51  0.783461511  4.8024603
    layer.52  0.806382924  4.0421516
    layer.53  0.829241335  9.3756512
    layer.54  0.849924415  5.3305050
    layer.55  0.871352434  7.5445803
    layer.56  0.893632233  6.4679547
    layer.57  0.916052133  2.8096065
    layer.58  0.938579470  5.3921661
    layer.59  0.959907651  7.2043689
    layer.60  0.981643587  3.3350806
    layer.61  1.004116774  8.8690707
    layer.62  1.028363466  1.7861299
    layer.63  1.054009140  6.2555038
    layer.64  1.072440803  7.6079236
    layer.65  1.094457805  7.6871483
    layer.66  1.123176277  4.7787764
    layer.67  1.149430871 12.7110502
    layer.68  1.170912921 -0.7156284
    layer.69  1.196743071  1.6490899
    layer.70  1.218625903  3.0363024
    layer.71  1.241868377  4.2974769
    layer.72  1.267941594  1.9543778
    layer.73  1.290708780  3.9986964
    layer.74  1.313222289  4.5179472
    layer.75  1.339045882  0.9337905
    layer.76  1.362803459  3.3050770
    layer.77  1.384450197  3.5422970
    layer.78  1.409720302  5.9973660
    layer.79  1.435851157  0.5081869
    layer.80  1.455592215  7.9661630
    layer.81  1.479495347  9.9460496
    layer.82  1.506051958  3.7908372
    layer.83  1.525728464  2.5735847
    layer.84  1.549362063 10.1404974
    layer.85  1.573440671 13.7408304
    layer.86  1.600278735  0.9335771
    layer.87  1.623879492  9.7588742
    layer.88  1.650029302  1.2769395
    layer.89  1.672362328 13.4970906
    layer.90  1.700221121 10.2087502
    layer.91  1.724793375  1.6811275
    layer.92  1.751070559  6.1178992
    layer.93  1.778022110 -0.1567626
    layer.94  1.803022087  3.8237479
    layer.95  1.830668867  4.4331468
    layer.96  1.855736911  5.9790707
    layer.97  1.882615030 11.3104333
    layer.98  1.909218490  8.2142607
    layer.99  1.938130021 15.3209674
    layer.100 1.963727593  5.8178217
    layer.101 1.993271947  9.6004907
    layer.102 2.022548139  3.4063646
    layer.103 2.050679922  4.7375010
    layer.104 2.078064442  3.0133019
    layer.105 2.104113460  5.5659522
    layer.106 2.133597612 12.0346333
    layer.107 2.164026260 -0.4028320
    layer.108 2.194852829 10.5996780
    layer.109 2.224257946  5.4479584
    layer.110 2.252194643  4.7052374
    layer.111 2.277335048 14.0962019
    layer.112 2.304058313  5.7149016
    layer.113 2.330930233  3.7780072
    layer.114 2.357022762  4.4120620
    layer.115 2.386489272  4.1866085
    layer.116 2.417503953  6.9078802
    layer.117 2.448524356  2.7825739
    layer.118 2.478698969  7.6171786
    layer.119 2.510175705 10.2410603
    layer.120 2.539697886  8.1820711
    layer.121 2.567915559  4.8275494
    layer.122 2.597463250 19.1624883
    layer.123 2.627518773 16.0677109
    layer.124 2.658759236 12.5897081
    layer.125 2.692401528  9.2907988
    layer.126 2.721903205  7.4262502
    layer.127 2.753021359  9.3902518
    layer.128 2.786313415 12.6193550
    layer.129 2.819564104 11.1121040
    layer.130 2.850823164 15.7907100
    layer.131 2.880394101 10.7425287
    layer.132 2.911391258  7.7971430
    layer.133 2.942965150  8.8060858
    layer.134 2.974468350 17.5606266
    layer.135 3.008983612 17.3088605
    layer.136 3.040015221 13.4500543
    layer.137 3.072668672 14.6377884
    layer.138 3.105982423  8.0798552dput(onepctCO2MEDIAN)    dput(onepctCO2MEDIAN)
    structure(list(x = c(0, 0.00679444684647024, 0.014288058038801, 
    0.0220879195258021, 0.0307973567396402,0.0384510718286037,0.0480879042297602, 
    0.0586777292191982, 0.0692614056169987, 0.080524530261755,0.0927602462470531, 
    0.103789608925581, 0.116953168064356, 0.129253298044205, 0.141710050404072, 
    0.156002052128315, 0.170648172497749, 0.185318425297737, 0.199463054537773, 
    0.21351333707571, 0.22883927077055, 0.246981292963028, 0.263012766838074, 
    0.278505563735962, 0.29365836083889, 0.310747265815735, 0.325990349054337, 
    0.342517539858818, 0.362751632928848, 0.380199536681175, 0.39499294757843, 
    0.414373397827148, 0.430690214037895, 0.449738144874573, 0.470167458057404, 
    0.489019870758057, 0.507242470979691, 0.524314284324646, 0.543750524520874, 
    0.56423419713974, 0.583679616451263, 0.601459443569183, 0.619924664497375, 
    0.639932006597519, 0.661347180604935, 0.684117317199707, 0.704829752445221, 
    0.725045770406723, 0.745165824890137, 0.765016138553619, 0.783461511135101, 
    0.806382924318314, 0.829241335391998, 0.84992441534996, 0.871352434158325, 
    0.893632233142853, 0.916052132844925, 0.938579469919205, 0.959907650947571, 
    0.981643587350845, 1.00411677360535, 1.02836346626282, 1.05400913953781, 
    1.07244080305099, 1.09445780515671, 1.12317627668381, 1.14943087100983, 
    1.17091292142868, 1.19674307107925, 1.21862590312958, 1.24186837673187, 
    1.26794159412384, 1.2907087802887, 1.31322228908539, 1.33904588222504, 
    1.36280345916748, 1.38445019721985, 1.40972030162811, 1.43585115671158, 
    1.45559221506119, 1.47949534654617, 1.50605195760727, 1.52572846412659, 
    1.5493620634079, 1.5734406709671, 1.60027873516083, 1.62387949228287, 
    1.65002930164337, 1.67236232757568, 1.70022112131119, 1.72479337453842, 
    1.75107055902481, 1.77802211046219, 1.80302208662033, 1.83066886663437, 
    1.85573691129684, 1.88261502981186, 1.90921849012375, 1.93813002109528, 
    1.96372759342194, 1.99327194690704, 2.02254813909531, 2.05067992210388, 
    2.07806444168091, 2.1041134595871, 2.13359761238098, 2.16402626037598, 
    2.19485282897949, 2.2242579460144, 2.25219464302063, 2.27733504772186, 
    2.30405831336975, 2.33093023300171, 2.35702276229858, 2.38648927211761, 
    2.41750395298004, 2.44852435588837, 2.47869896888733, 2.51017570495605, 
    2.53969788551331, 2.567915558815, 2.59746325016022, 2.62751877307892, 
    2.65875923633575, 2.69240152835846, 2.72190320491791, 2.75302135944366, 
    2.78631341457367, 2.8195641040802, 2.85082316398621, 2.88039410114288, 
    2.91139125823975, 2.94296514987946, 2.97446835041046, 3.00898361206055, 
    3.04001522064209, 3.07266867160797, 3.10598242282867), y = c(0, 
    4.90024901723162, 0.160799993152722, 6.63491326258641, -1.24295055804536, 
    1.56433744259162, -2.26590352245208, 2.20700446463354, -2.36770012911069, 
    -1.09135061899174, 0.409993989292701, -0.125972681525582, -2.41382533818026, 
    7.08902570153028, -0.759353880417294, 0.0454415959640926, -1.53496826259972, 
    6.55242014096194, -0.831256280861552, -2.50991825629084, 0.136596820654013, 
    -1.37198445498419, -0.871298832596736, 0.663258363762466, 0.793803634291308, 
    3.48806373666998, -4.46122081238949, 0.0871733966938564, -1.41715777257774, 
    -0.995650815648318, 0.32155262317503, 3.14038657369241, -0.737609879885404, 
    -2.48605406511292, -3.423585843908, 0.482474753780281, -0.978538630093809, 
    8.53596837794201, 5.48447420320695, 3.21493665820644, 3.91689160157513, 
    4.49070195980797, 6.54104103157039, 4.80686500146557, 8.15101701282067, 
    0.26974132191657, -0.180750068063062, 9.71812491230244, 1.54064657400204, 
    -1.64760408795688, 4.80246028991894, 4.04215159914344, 9.37565121768513, 
    5.33050496938428, 7.54458026088508, 6.46795470819342, 2.80960651433971, 
    5.39216613235986, 7.20436888038562, 3.3350806460997, 8.86907069895943, 
    1.78612988613659, 6.25550382050395, 7.60792364896564, 7.68714830528144, 
    4.77877638957615, 12.7110501777314, -0.715628443181046, 1.64908991824022, 
    3.03630240714679, 4.29747688442346, 1.95437780501881, 3.99869636910933, 
    4.51794724689848, 0.933790484492299, 3.30507700050003, 3.5422970157433, 
    5.99736597322524, 0.508186860060022, 7.96616300581067, 9.94604963036295, 
    3.79083717222623, 2.57358468532258, 10.1404974171776, 13.7408303595752, 
    0.933577123801399, 9.75887417074129, 1.27693947132921, 13.4970905965787, 
    10.2087501765735, 1.68112753028756, 6.1178991508927, -0.156762622680077, 
    3.82374791691426, 4.43314678736265, 5.97907067167507, 11.3104332518482, 
    8.21426074201525, 15.320967360602, 5.81782169471483, 9.6004907412354, 
    3.40636455909704, 4.73750103921864, 3.0133019468806, 5.56595224859066, 
    12.0346332527215, -0.40283199827104, 10.5996779538754, 5.44795836991128, 
    4.70523736412729, 14.096201892183, 5.71490161813391, 3.77800720810782, 
    4.41206200639436, 4.18660847858423, 6.90788020044911, 2.78257393345915, 
    7.61717857379431, 10.2410602647684, 8.18207106836167, 4.82754943871433, 
    19.1624882857155, 16.0677109398509, 12.589708067017, 9.29079879799404, 
    7.42625019725314, 9.39025179806185, 12.6193550331438, 11.1121039747257, 
    15.7907099734986, 10.7425286789233, 7.79714300307344, 8.80608578166101, 
    17.5606266346039, 17.3088604929222, 13.4500543478523, 14.6377884248645, 
    8.07985518296064)), class = "data.frame", row.names = c("layer.1", 
    "layer.2", "layer.3", "layer.4", "layer.5", "layer.6", "layer.7", 
    "layer.8", "layer.9", "layer.10", "layer.11", "layer.12", "layer.13", 
    "layer.14", "layer.15", "layer.16", "layer.17", "layer.18", "layer.19", 
    "layer.20", "layer.21", "layer.22", "layer.23", "layer.24", "layer.25", 
    "layer.26", "layer.27", "layer.28", "layer.29", "layer.30", "layer.31", 
    "layer.32", "layer.33", "layer.34", "layer.35", "layer.36", "layer.37", 
    "layer.38", "layer.39", "layer.40", "layer.41", "layer.42", "layer.43", 
    "layer.44", "layer.45", "layer.46", "layer.47", "layer.48", "layer.49", 
    "layer.50", "layer.51", "layer.52", "layer.53", "layer.54", "layer.55", 
    "layer.56", "layer.57", "layer.58", "layer.59", "layer.60", "layer.61", 
    "layer.62", "layer.63", "layer.64", "layer.65", "layer.66", "layer.67", 
    "layer.68", "layer.69", "layer.70", "layer.71", "layer.72", "layer.73", 
    "layer.74", "layer.75", "layer.76", "layer.77", "layer.78", "layer.79", 
    "layer.80", "layer.81", "layer.82", "layer.83", "layer.84", "layer.85", 
    "layer.86", "layer.87", "layer.88", "layer.89", "layer.90", "layer.91", 
    "layer.92", "layer.93", "layer.94", "layer.95", "layer.96", "layer.97", 
    "layer.98", "layer.99", "layer.100", "layer.101", "layer.102", 
    "layer.103", "layer.104", "layer.105", "layer.106", "layer.107", 
    "layer.108", "layer.109", "layer.110", "layer.111", "layer.112", 
    "layer.113", "layer.114", "layer.115", "layer.116", "layer.117", 
    "layer.118", "layer.119", "layer.120", "layer.121", "layer.122", 
    "layer.123", "layer.124", "layer.125", "layer.126", "layer.127", 
    "layer.128", "layer.129", "layer.130", "layer.131", "layer.132", 
    "layer.133", "layer.134", "layer.135", "layer.136", "layer.137", 
    "layer.138"))
I started with the following to generate the first regression line and scatter plot:    lm<-ggplot(onepctCO2MEDIAN) +  
    geom_jitter(aes(RCP1pctCO2cumulativeMedian[1:138], departurea),   
    colour="blue") + geom_smooth(aes(RCP1pctCO2cumulativeMedian[1:138], 
    departurea), method=lm) 
But I receive this error:    Warning message:
    Computation failed in `stat_smooth()`:
    'what' must be a function or character string
A blue scatter plot is successfully generated, but the problem is that the regression line does not appear, presumably related to the above warning.
Is there a reason for this? I would appreciate any assistance!
	[[alternative HTML version deleted]]



More information about the R-sig-Geo mailing list