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@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.
