On 07/20/2018 01:28 PM, Bert Gunter wrote:> Nothing attached. The mail server strips most attachments for security.Here is my repex, inline: library(ggplot2) conv_df <- structure(list(Adjustment = structure(c(1L, 2L, 3L, 4L, 5L, 1L, ??????????????????????????????????????? 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, ??????????????????????????????????????? 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, ??????????????????????????????????????? 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, ??????????????????????????????????????? 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, ??????????????????????????????????????? 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, ??????????????????????????????????????? 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, ??????????????????????????????????????? 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, ??????????????????????????????????????? 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, ??????????????????????????????????????? 5L), .Label = c("Radar", "MeanFieldBias", "Multiplicative", "Mixed", "ConditionalMerge"), class = "factor"), 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.", "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", "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("Adjustment", "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, 181L, 183L, 184L, 185L, 186L, 191L, 193L, 194L, 195L, 196L, 201L, 203L, 204L, 205L, 206L, 211L, 213L, 214L, 215L, 216L, 221L, 223L, 224L, 225L, 226L, 231L, 233L, 234L, 235L, 236L, 241L, 243L, 244L, 245L, 246L, 251L, 253L, 254L, 255L, 256L, 261L, 263L, 264L, 265L, 266L), class = "data.frame") strat_df <- structure(list(Adjustment = structure(c(1L, 2L, 3L, 4L, 5L, 1L, ??????????????????????????????????????????????????? 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, ??????????????????????????????????????????????????? 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, ??????????????????????????????????????????????????? 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, ??????????????????????????????????????????????????? 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, ??????????????????????????????????????????????????? 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, ??????????????????????????????????????????????????? 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, ??????????????????????????????????????????????????? 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, ??????????????????????????????????????????????????? 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, ??????????????????????????????????????????????????? 5L), .Label = c("Radar", "MeanFieldBias", "Multiplicative", "Mixed", "ConditionalMerge"), class = "factor"), 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.", "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", "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("Adjustment", "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, 182L, 187L, 188L, 189L, 190L, 192L, 197L, 198L, 199L, 200L, 202L, 207L, 208L, 209L, 210L, 212L, 217L, 218L, 219L, 220L, 222L, 227L, 228L, 229L, 230L, 232L, 237L, 238L, 239L, 240L, 242L, 247L, 248L, 249L, 250L, 252L, 257L, 258L, 259L, 260L, 262L, 267L, 268L, 269L, 270L), class = "data.frame") ggplot(strat_df, mapping=aes(x=Intercept, y=Slope)) + ? geom_point(aes(color=Struct, shape=Adjustment), size=6) ggplot(strat_df, mapping=aes(x=Intercept, y=Slope)) + ? geom_point(aes(color=Error, shape=Adjustment), size=6) ggplot(conv_df, mapping=aes(x=Intercept, y=Slope)) + ? geom_point(aes(color=Struct, shape=Adjustment), size=6) ggplot(conv_df, mapping=aes(x=Intercept, y=Slope)) + ? geom_point(aes(color=Error, shape=Adjustment), size=6)> > See the posting guide below and ?dput for how to include data. > + We need your faulty code also, of course. > > Cheers, > Bert > > > > Bert Gunter > > "The trouble with having an open mind is that people keep coming along > and sticking things into it." > -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) > > On Fri, Jul 20, 2018 at 3:15 AM, Micha Silver <tsvibar at gmail.com > <mailto:tsvibar at gmail.com>> wrote: > > Hello: > > I have two data frames (subsetted from a larger one). I am > plotting scatterplots with ggplot2 and coloring by different > variables. When using one of the variables "Error" with one of the > dataframes "conv_df" only a single color is displayed. All the > other variables show colors as expected, and when applying the > same ggplot() functions to the second dataframe, colors are > correctly applied for all variables. > > I'm scratching me head over this for some time. Maybe someone can > see what I'm missing? > > My repex is attached > > Thanks > > -- > Micha Silver > Ben Gurion Univ. > Sde Boker, Remote Sensing Lab > cell: +972-523-665918 > > ______________________________________________ > R-help at r-project.org <mailto:R-help at r-project.org> mailing list -- > To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > <https://stat.ethz.ch/mailman/listinfo/r-help> > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > <http://www.R-project.org/posting-guide.html> > and provide commented, minimal, self-contained, reproducible code. > >-- Micha Silver Ben Gurion Univ. Sde Boker, Remote Sensing Lab cell: +972-523-665918
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
Hi Values in conv_df are almost same for each Error level> aggregate(conv_df$Intercept, list(conv_df$Error), mean)Group.1 x 1 High CML error -3.226313 2 Low CML error -3.226536 3 Med CML error -3.226422> aggregate(conv_df$Slope, list(conv_df$Error), mean)Group.1 x 1 High CML error 8.325558 2 Low CML error 8.324242 3 Med CML error 8.324721> boxplot(split(conv_df$Slope, conv_df$Error)) > boxplot(split(conv_df$Intercept, conv_df$Error)) >so points in ggplot are overplotted and only colour for last plotted level is visible. Cheers Petr> 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 > > ______________________________________________ > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code.Osobn? ?daje: Informace o zpracov?n? a ochran? osobn?ch ?daj? obchodn?ch partner? PRECHEZA a.s. jsou zve?ejn?ny na: https://www.precheza.cz/zasady-ochrany-osobnich-udaju/ | Information about processing and protection of business partner?s personal data are available on website: https://www.precheza.cz/en/personal-data-protection-principles/ D?v?rnost: Tento e-mail a jak?koliv k n?mu p?ipojen? dokumenty jsou d?v?rn? a podl?haj? tomuto pr?vn? z?vazn?mu prohl??en? o vylou?en? odpov?dnosti: https://www.precheza.cz/01-dovetek/ | This email and any documents attached to it may be confidential and are subject to the legally binding disclaimer: https://www.precheza.cz/en/01-disclaimer/