[R] scatter plot coloring problem
Micha Silver
t@v|b@r @end|ng |rom gm@||@com
Fri Jul 20 13:27:36 CEST 2018
On 07/20/2018 02:01 PM, Micha Silver wrote:
>
>
> On 07/20/2018 01:28 PM, Bert Gunter wrote:
>> Nothing attached. The mail server strips most attachments for security.
Here is a (slightly shorter) repex.
The fourth plot shows the problem. The first three are as expected.
library(ggplot2)
strat_df <- structure(list(Slope = c(1.15639473681994, 0.972278300619073,
1.14313365332712,
1.14399372216612, 1.22529134790727, 2.14326711679831,
0.54859156211142,
1.17046713623601, 1.17453878322687, 3.31323770780669,
0.966673750291528,
0.746319725592914, 1.02480114419885, 1.0229752524756,
1.00132723720128,
1.15639473681994, 0.972278300619073, 1.14313365332712,
1.14399372216612,
1.22529134790727, 2.14326711679831, 0.54859156211142,
1.17046713623601,
1.17453878322687, 3.31323770780669, 0.966673750291528,
0.746319725592914,
1.02480114419885, 1.0229752524756, 1.00132723720128,
1.15639473681994,
0.972278300619073, 1.14313365332712, 1.14399372216612,
1.22529134790727,
2.14326711679831, 0.54859156211142, 1.17046713623601,
1.17453878322687,
3.31323770780669, 0.966673750291528, 0.746319725592914,
1.02480114419885, 1.0229752524756,
1.00132723720128, 1.00564012971235,
0.781822304249351,
0.793232303666035, 0.798538338102623, 1.20657643535607,
2.16318241002351,
1.42268702026443, 0.812412401599228, 0.817089629984966,
3.27413789688355,
0.840735716836696, 0.629583577106973, 0.710708013029288,
0.714002133616219,
0.985501434987176, 1.03494070727933, 0.956962513783496,
0.793549829222589,
0.798851646709461, 1.20660885706037, 2.18674952848921,
1.5634217520085, 0.81273050967827,
0.817407079251797, 3.27417031858785,
0.877344267863596,
0.848664941338144, 0.710978765510023, 0.714268901361411,
0.985533856691474,
1.03983096388345, 0.985084623809766, 0.794138608594684,
0.79943259180892,
1.20665856218731, 2.21336885559571, 1.72189603514386,
0.813327118056704,
0.818002445250057, 3.27422002371479, 0.877473770321168,
0.848389484199237,
0.711483989718558, 0.714766443910276, 0.985583561818414
), Intercept = c(-1.0892464119326, 0.246565490954958, -1.1639008288321,
-1.17146025612969, -1.93566050996642,
-3.96706179209372, 3.9202695728101,
-1.04771345550341, -1.07809518302298,
-19.7820283549024, 0.36822137463493,
2.207658818211, -0.016876438597312,
-0.00842619567899302, 0.0306144651329568,
-1.0892464119326, 0.246565490954958, -1.1639008288321,
-1.17146025612969,
-1.93566050996642, -3.96706179209372, 3.9202695728101,
-1.04771345550341,
-1.07809518302298, -19.7820283549024,
0.36822137463493, 2.207658818211,
-0.016876438597312, -0.00842619567899302,
0.0306144651329568,
-1.0892464119326, 0.246565490954958, -1.1639008288321,
-1.17146025612969,
-1.93566050996642, -3.96706179209372, 3.9202695728101,
-1.04771345550341,
-1.07809518302298, -19.7820283549024,
0.36822137463493, 2.207658818211,
-0.016876438597312, -0.00842619567899302,
0.0306144651329568,
-1.42588814587741, -0.0938807577163176,
-0.760427987439108, -0.802681471353857,
-4.38731629113287, -4.72258266663199,
1.08443985346667, -0.679055757307681,
-0.713928732007598, -22.2121480319997,
-0.187477151457156, 1.28512191810652,
0.0386762170827232, 0.00928219285158222,
-2.45430083552569, -1.5052623193495,
-0.571094020563348, -0.760898951384971,
-0.80312079180565, -4.38543778472072,
-4.74043340113669, 0.976337655364814,
-0.679498378418568, -0.714366826235178,
-22.2102695255875, -0.264069816907669,
0.823783304969757, 0.0386132206604378,
0.00924930329540931, -2.45242232911355,
-1.5204762452006, -0.664489263100851,
-0.761704010560802, -0.803866601426402,
-4.38177522849324, -4.78366916314869,
0.714779631860822, -0.68026224065537,
-0.715122068793879, -22.20660696936,
-0.262931522643713, 0.827039881514761,
0.0385397097353969, 0.00923344703169212,
-2.44875977288606), Struct = c("Bias only", "Bias
only", "Bias only",
"Bias only", "Bias
only", "Lognormal error dist.", "Lognormal error dist.",
"Lognormal error
dist.", "Lognormal error dist.", "Lognormal error dist.",
"FFT-based struct.",
"FFT-based struct.", "FFT-based struct.",
"FFT-based struct.",
"FFT-based struct.", "Bias only", "Bias only",
"Bias only", "Bias
only", "Bias only", "Lognormal error dist.",
"Lognormal error
dist.", "Lognormal error dist.", "Lognormal error dist.",
"Lognormal error
dist.", "FFT-based struct.", "FFT-based struct.",
"FFT-based struct.",
"FFT-based struct.", "FFT-based struct.",
"Bias only", "Bias
only", "Bias only", "Bias only", "Bias only",
"Lognormal error
dist.", "Lognormal error dist.", "Lognormal error dist.",
"Lognormal error
dist.", "Lognormal error dist.", "FFT-based struct.",
"FFT-based struct.",
"FFT-based struct.", "FFT-based struct.",
"FFT-based struct.",
"Bias only", "Bias only", "Bias only", "Bias only",
"Bias only", "Lognormal
error dist.", "Lognormal error dist.",
"Lognormal error
dist.", "Lognormal error dist.", "Lognormal error dist.",
"FFT-based struct.",
"FFT-based struct.", "FFT-based struct.",
"FFT-based struct.",
"FFT-based struct.", "Bias only", "Bias only",
"Bias only", "Bias
only", "Bias only", "Lognormal error dist.",
"Lognormal error
dist.", "Lognormal error dist.", "Lognormal error dist.",
"Lognormal error
dist.", "FFT-based struct.", "FFT-based struct.",
"FFT-based struct.",
"FFT-based struct.", "FFT-based struct.",
"Bias only", "Bias
only", "Bias only", "Bias only", "Bias only",
"Lognormal error
dist.", "Lognormal error dist.", "Lognormal error dist.",
"Lognormal error
dist.", "Lognormal error dist.", "FFT-based struct.",
"FFT-based struct.",
"FFT-based struct.", "FFT-based struct.",
"FFT-based struct."),
Error = c("Low CML error", "Low CML error",
"Low CML error", "Low CML error", "Low CML error", "Low CML error",
"Low CML error", "Low CML error", "Low CML error", "Low CML error",
"Low CML error", "Low CML error", "Low CML error", "Low CML error",
"Low CML error", "Med CML error", "Med CML error", "Med CML error",
"Med CML error", "Med CML error", "Med CML error", "Med CML error",
"Med CML error", "Med CML error", "Med CML error", "Med CML error",
"Med CML error", "Med CML error", "Med CML error", "Med CML error",
"High CML error", "High CML error", "High CML error", "High CML error",
"High CML error", "High CML error", "High CML error", "High CML error",
"High CML error", "High CML error", "High CML error", "High CML error",
"High CML error", "High CML error", "High CML error", "Low CML error",
"Low CML error", "Low CML error", "Low CML error", "Low CML error",
"Low CML error", "Low CML error", "Low CML error", "Low CML error",
"Low CML error", "Low CML error", "Low CML error", "Low CML error",
"Low CML error", "Low CML error", "Med CML error", "Med CML error",
"Med CML error", "Med CML error", "Med CML error", "Med CML error",
"Med CML error", "Med CML error", "Med CML error", "Med CML error",
"Med CML error", "Med CML error", "Med CML error", "Med CML error",
"Med CML error", "High CML error", "High CML error", "High CML error",
"High CML error", "High CML error", "High CML error", "High CML error",
"High CML error", "High CML error", "High CML error", "High CML error",
"High CML error", "High CML error", "High CML error", "High CML error"
)), .Names = c("Slope",
"Intercept", "Struct", "Error"), row.names = c(2L,
7L, 8L, 9L, 10L, 12L, 17L, 18L, 19L, 20L, 22L, 27L, 28L, 29L,
30L, 32L, 37L, 38L, 39L, 40L, 42L, 47L, 48L, 49L, 50L, 52L, 57L,
58L, 59L, 60L, 62L, 67L, 68L, 69L, 70L, 72L, 77L, 78L, 79L, 80L,
82L, 87L, 88L, 89L, 90L, 92L, 97L, 98L, 99L, 100L, 102L, 107L,
108L, 109L, 110L, 112L, 117L, 118L, 119L, 120L, 122L, 127L, 128L,
129L, 130L, 132L, 137L, 138L, 139L, 140L, 142L, 147L, 148L, 149L,
150L, 152L, 157L, 158L, 159L, 160L, 162L, 167L, 168L, 169L, 170L,
172L, 177L, 178L, 179L, 180L), class = "data.frame")
conv_df <- structure(list(Slope = c(1.04807118675076, 0.998948448847309,
0.999431048648743,
1.06271129629311, 1.07572617384751, 2.22760090861169,
0.932482182264868,
1.02936746325801, 2.52988042582817, 2.95633562179131,
18.8227507300988,
0.974588681842621, 108.07393074932, 0.870832606501646,
1.00567956715787, 1.04807118675076,
0.998948448847309, 0.999431048648743,
1.06271129629311, 1.07572617384751,
2.22760090861169, 0.932482182264868,
1.02936746325801, 2.52988042582817,
2.95633562179131, 18.8227507300988,
0.974588681842621, 108.07393074932,
0.870832606501646, 1.00567956715787,
1.04807118675076,
0.998948448847309, 0.999431048648743, 1.06271129629311,
1.07572617384751, 2.22760090861169,
0.932482182264868, 1.02936746325801,
2.52988042582817, 2.95633562179131,
18.8227507300988, 0.974588681842621,
108.07393074932, 0.870832606501646,
1.00567956715787, 0.938306957617447,
0.696202687574185,
0.698287792550519, 1.01157976748338, 1.07216134522869,
2.11775044458303,
0.658564756192898, 0.718936009019651, 2.46328731519796,
2.94817482855984, 13.4010203977042,
0.758075017209242, 75.8087725102359,
0.825432059460383,
1.00237002354864, 0.93840246618556, 0.696484580425202,
0.698549805163105,
1.01161728893021, 1.07215296972693, 2.11785160868381,
0.658846969311393,
0.719217656820784, 2.46333881438479, 2.94816645305808,
13.4028764261544,
0.758261665458076, 75.8196908351967, 0.82547163245355,
1.00236164804689,
0.938596807939808, 0.697064148960971, 0.699062697365857,
1.01169996923602, 1.07214387920809,
2.11805571501552, 0.659427992302351,
0.719759915668507,
2.46344926105519, 2.94815736253924, 13.4060642247861,
0.758655662542348,
75.8383487496718, 0.825555603203494, 1.00235255752804
), Intercept = c(-0.017716276271862, 0.000961091319158938,
-0.000703510322277703,
-0.00845866525088005, -0.0966678762281092,
-0.389386368120609,
0.068730059320282, 0.0367676804603267,
-0.0134454494180847, -2.61898714802878,
-7.30136425060206, -0.21274280699377,
-43.1291289782452, -0.125688969456484,
-0.307472718728007, -0.017716276271862,
0.000961091319158938,
-0.000703510322277703, -0.00845866525088005,
-0.0966678762281092,
-0.389386368120609, 0.068730059320282,
0.0367676804603267, -0.0134454494180847,
-2.61898714802878, -7.30136425060206,
-0.21274280699377, -43.1291289782452,
-0.125688969456484, -0.307472718728007,
-0.017716276271862, 0.000961091319158938,
-0.000703510322277703, -0.00845866525088005,
-0.0966678762281092,
-0.389386368120609, 0.068730059320282,
0.0367676804603267, -0.0134454494180847,
-2.61898714802878, -7.30136425060206,
-0.21274280699377, -43.1291289782452,
-0.125688969456484, -0.307472718728007,
-0.0859314660447012,
0.0135310033957625, -3.43619764283222e-05,
-0.0523622551444108,
-0.475294485394066, -0.460272581935994,
0.0625183970346825, 0.0262518388075903,
-0.0839509467307793, -2.95707142967006,
-5.52875677718304, -0.152998423411686,
-32.1368671911073, -0.164063261306885,
-0.685459757083477, -0.0858163634696177,
0.01352935620949, 9.46328483051025e-05,
-0.052203322505647, -0.47489015022079,
-0.460135002763724, 0.0626429365104076,
0.02640499766902, -0.0838075052075929,
-2.95666709449679, -5.52852300483326,
-0.153159000797531, -32.135862357057,
-0.163909219046001, -0.685055421910201,
-0.0856060516008676,
0.0135479894875152, 0.000302270025863519,
-0.0519314642867551,
-0.474126407682765, -0.45988984866136,
0.0628609678242435, 0.026651390571036,
-0.0835647473472818, -2.95590335195876,
-5.52851009407545, -0.153425312176928,
-32.1365449540903, -0.163646588624356,
-0.684291679372176), Struct = c("Bias only",
"Bias only", "Bias only", "Bias only", "Bias only", "Lognormal error dist.",
"Lognormal error dist.", "Lognormal error dist.", "Lognormal error dist.",
"Lognormal error dist.", "FFT-based struct.", "FFT-based struct.",
"FFT-based struct.", "FFT-based struct.", "FFT-based struct.",
"Bias only", "Bias only", "Bias only", "Bias only", "Bias only",
"Lognormal error dist.", "Lognormal error dist.", "Lognormal error dist.",
"Lognormal error dist.", "Lognormal error dist.", "FFT-based struct.",
"FFT-based struct.", "FFT-based struct.", "FFT-based struct.",
"FFT-based struct.", "Bias only", "Bias only", "Bias only", "Bias only",
"Bias only", "Lognormal error dist.", "Lognormal error dist.",
"Lognormal error dist.", "Lognormal error dist.", "Lognormal error dist.",
"FFT-based struct.", "FFT-based struct.", "FFT-based struct.",
"FFT-based struct.", "FFT-based struct.", "Bias only", "Bias only",
"Bias only", "Bias only", "Bias only", "Lognormal error dist.",
"Lognormal error dist.", "Lognormal error dist.", "Lognormal error dist.",
"Lognormal error dist.", "FFT-based struct.", "FFT-based struct.",
"FFT-based struct.", "FFT-based struct.", "FFT-based struct.",
"Bias only", "Bias only", "Bias only", "Bias only", "Bias only",
"Lognormal error dist.", "Lognormal error dist.", "Lognormal error dist.",
"Lognormal error dist.", "Lognormal error dist.", "FFT-based struct.",
"FFT-based struct.", "FFT-based struct.", "FFT-based struct.",
"FFT-based struct.", "Bias only", "Bias only", "Bias only", "Bias only",
"Bias only", "Lognormal error dist.", "Lognormal error dist.",
"Lognormal error dist.", "Lognormal error dist.", "Lognormal error dist.",
"FFT-based struct.", "FFT-based struct.", "FFT-based struct.",
"FFT-based struct.", "FFT-based struct."), Error = c("Low CML error",
"Low CML error", "Low CML error", "Low CML error", "Low CML error",
"Low CML error", "Low CML error", "Low CML error", "Low CML error",
"Low CML error", "Low CML error", "Low CML error", "Low CML error",
"Low CML error", "Low CML error", "Med CML error", "Med CML error",
"Med CML error", "Med CML error", "Med CML error", "Med CML error",
"Med CML error", "Med CML error", "Med CML error", "Med CML error",
"Med CML error", "Med CML error", "Med CML error", "Med CML error",
"Med CML error", "High CML error", "High CML error", "High CML error",
"High CML error", "High CML error", "High CML error", "High CML error",
"High CML error", "High CML error", "High CML error", "High CML error",
"High CML error", "High CML error", "High CML error", "High CML error",
"Low CML error", "Low CML error", "Low CML error", "Low CML error",
"Low CML error", "Low CML error", "Low CML error", "Low CML error",
"Low CML error", "Low CML error", "Low CML error", "Low CML error",
"Low CML error", "Low CML error", "Low CML error", "Med CML error",
"Med CML error", "Med CML error", "Med CML error", "Med CML error",
"Med CML error", "Med CML error", "Med CML error", "Med CML error",
"Med CML error", "Med CML error", "Med CML error", "Med CML error",
"Med CML error", "Med CML error", "High CML error", "High CML error",
"High CML error", "High CML error", "High CML error", "High CML error",
"High CML error", "High CML error", "High CML error", "High CML error",
"High CML error", "High CML error", "High CML error", "High CML error",
"High CML error")), .Names = c("Slope", "Intercept", "Struct",
"Error"), row.names = c(1L, 3L, 4L, 5L, 6L, 11L, 13L, 14L, 15L,
16L, 21L, 23L, 24L, 25L, 26L, 31L, 33L, 34L, 35L, 36L, 41L, 43L,
44L, 45L, 46L, 51L, 53L, 54L, 55L, 56L, 61L, 63L, 64L, 65L, 66L,
71L, 73L, 74L, 75L, 76L, 81L, 83L, 84L, 85L, 86L, 91L, 93L, 94L,
95L, 96L, 101L, 103L, 104L, 105L, 106L, 111L, 113L, 114L, 115L,
116L, 121L, 123L, 124L, 125L, 126L, 131L, 133L, 134L, 135L, 136L,
141L, 143L, 144L, 145L, 146L, 151L, 153L, 154L, 155L, 156L, 161L,
163L, 164L, 165L, 166L, 171L, 173L, 174L, 175L, 176L), class = "data.frame")
pl.strat.1 <- ggplot(strat_df, mapping=aes(x=Intercept, y=Slope)) +
geom_point(aes(color=Struct),size=4) +
scale_y_continuous(limits=c(0.3, 3.0)) +
scale_x_continuous(limits=c(-7, 7))
pl.strat.2 <- ggplot(strat_df, mapping=aes(x=Intercept, y=Slope)) +
geom_point(aes(color=Error), size=4) +
scale_y_continuous(limits=c(0.3, 3.0)) +
scale_x_continuous(limits=c(-7, 7))
pl.conv.1 <- ggplot(conv_df, mapping=aes(x=Intercept, y=Slope)) +
geom_point(aes(color=Struct), size=4) +
scale_y_continuous(limits=c(0.3, 3.0)) +
scale_x_continuous(limits=c(-7, 7))
pl.conv.2 <- ggplot(conv_df, mapping=aes(x=Intercept, y=Slope)) +
geom_point(aes(color=Error), size=4) +
scale_y_continuous(limits=c(0.3, 3.0)) +
scale_x_continuous(limits=c(-7, 7))
print(pl.strat.1)
print(pl.strat.2)
print(pl.conv.1)
print(pl.conv.2)
>
>
--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918
More information about the R-help
mailing list