thr3ads.net - R help - [R] contouring x y scatter data [Oct 2011]

If this information is useful, please help other people find it:
Share via:

emorway

2011-Oct-17 18:00 UTC

[R] contouring x y scatter data

Hello,

I'm almost positive R can do the following, I just haven't hit upon the
right package or search terms, however.  Here's what I'm after: 
I've got
concentration output from two different models that I want to qualitatively
compare in a contour plot (or some variant of a contour plot).  The problem
as I see it is that the data is not gridded is the usual regular fashion,
and even if it were, the models are working with two different irregular
computational/numerical grids at which concentrations are calculated. 
Running the code below will produce a plot that will hopefully illuminate
what I'm attempting to describe.  At each location where 'O'
appears, there
is a calculated concentration that I would like to contour.  I'd like to do
the same thing at each 'X' location and then compare, visually, how
close
the contours are (hopefully close).  Both data.frames contain columns/fields
titled "ts181" and "ts1825" which contain the concentrations
at time step
181 days and time step 1825 days.  Is it possible to contour x-y-value data
if it is not laid out on a regular grid?


X<-c(0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000)

Z<-c(50,49.95,49.85,49.75,49.65,49.55,49.45,49.35,49.25,49.15,49.05,48.95,48.85,48.75,48.65,48.55,48.45,48.35,48.25,48.15,48.05,47.95,47.85,47.75,47.65,47.55,47.45,47.35,47.25,47.15,47.05,46.95,46.85,46.75,46.65,46.55,46.45,46.35,46.25,46.15,46.1,49.25,49.2,49.1,49,48.9,48.8,48.7,48.6,48.5,48.4,48.3,48.2,48.1,48,47.9,47.8,47.7,47.6,47.5,47.4,47.3,47.2,47.1,47,46.9,46.8,46.7,46.6,46.5,46.4,46.3,46.2,46.1,46,45.9,45.8,45.7,45.6,45.5,45.4,45.35,48.2,48.15,48.05,47.95,47.85,47.75,47.65,47.55,47.45,47.35,47.25,47.15,47.05,46.95,46.85,46.75,46.65,46.55,46.45,46.35,46.25,46.15,46.05,45.95,45.85,45.75,45.65,45.55,45.45,45.35,45.25,45.15,45.05,44.95,44.85,44.75,44.65,44.55,44.45,44.35,44.3,47,46.95,46.85,46.75,46.65,46.55,46.45,46.35,46.25,46.15,46.05,45.95,45.85,45.75,45.65,45.55,45.45,45.35,45.25,45.15,45.05,44.95,44.85,44.75,44.65,44.55,44.45,44.35,44.25,44.15,44.05,43.95,43.85,43.75,43.65,43.55,43.45,43.35,43.25,43.15,43.1,45.5,45.45,45.35,45.25,45.15,45.05,44.95,44.85,44.75,44.65,44.55,44.45,44.35,44.25,44.15,44.05,43.95,43.85,43.75,43.65,43.55,43.45,43.35,43.25,43.15,43.05,42.95,42.85,42.75,42.65,42.55,42.45,42.35,42.25,42.15,42.05,41.95,41.85,41.75,41.65,41.6,44,43.95,43.85,43.75,43.65,43.55,43.45,43.35,43.25,43.15,43.05,42.95,42.85,42.75,42.65,42.55,42.45,42.35,42.25,42.15,42.05,41.95,41.85,41.75,41.65,41.55,41.45,41.35,41.25,41.15,41.05,40.95,40.85,40.75,40.65,40.55,40.45,40.35,40.25,40.15,40.1,42.5,42.45,42.35,42.25,42.15,42.05,41.95,41.85,41.75,41.65,41.55,41.45,41.35,41.25,41.15,41.05,40.95,40.85,40.75,40.65,40.55,40.45,40.35,40.25,40.15,40.05,39.95,39.85,39.75,39.65,39.55,39.45,39.35,39.25,39.15,39.05,38.95,38.85,38.75,38.65,38.6,41,40.95,40.85,40.75,40.65,40.55,40.45,40.35,40.25,40.15,40.05,39.95,39.85,39.75,39.65,39.55,39.45,39.35,39.25,39.15,39.05,38.95,38.85,38.75,38.65,38.55,38.45,38.35,38.25,38.15,38.05,37.95,37.85,37.75,37.65,37.55,37.45,37.35,37.25,37.15,37.1,39.5,39.45,39.35,39.25,39.15,39.05,38.95,38.85,38.75,38.65,38.55,38.45,38.35,38.25,38.15,38.05,37.95,37.85,37.75,37.65,37.55,37.45,37.35,37.25,37.15,37.05,36.95,36.85,36.75,36.65,36.55,36.45,36.35,36.25,36.15,36.05,35.95,35.85,35.75,35.65,35.6,38,37.95,37.85,37.75,37.65,37.55,37.45,37.35,37.25,37.15,37.05,36.95,36.85,36.75,36.65,36.55,36.45,36.35,36.25,36.15,36.05,35.95,35.85,35.75,35.65,35.55,35.45,35.35,35.25,35.15,35.05,34.95,34.85,34.75,34.65,34.55,34.45,34.35,34.25,34.15,34.1,35,34.95,34.85,34.75,34.65,34.55,34.45,34.35,34.25,34.15,34.05,33.95,33.85,33.75,33.65,33.55,33.45,33.35,33.25,33.15,33.05,32.95,32.85,32.75,32.65,32.55,32.45,32.35,32.25,32.15,32.05,31.95,31.85,31.75,31.65,31.55,31.45,31.35,31.25,31.15,31.1)

ts181<-c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.113,0.1177,0.1167,0.1192,0.1242,0.1297,0.1317,0.131,0.1296,0.1296,0.1334,0.0014,0,0.0004,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.1074,0.1121,0.1111,0.1137,0.1187,0.1242,0.1262,0.1256,0.1241,0.1241,0.1279,0.0014,0,0.0004,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.0991,0.1039,0.1029,0.1053,0.1103,0.1159,0.1181,0.1175,0.116,0.116,0.1198,0.0015,0,0.0003,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.089,0.0939,0.0928,0.0951,0.1,0.1058,0.1082,0.1077,0.1063,0.1062,0.1101,0.0017,0,0.0003,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.0759,0.0811,0.0798,0.0816,0.0862,0.0923,0.0952,0.095,0.0936,0.0935,0.0973,0.0022,0,0.0003,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.0594,0.0659,0.0651,0.0674,0.0721,0.0785,0.0823,0.0825,0.0811,0.081,0.0846,0.0031,0,0.0002,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.0571,0.0639,0.0628,0.0632,0.0628,0.0616,0.0613,0.0622,0.0629,0.0646,0.0693,0.0046,0,0.0002,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.0561,0.063,0.0618,0.0619,0.061,0.0596,0.0592,0.0602,0.0609,0.0627,0.0676,0.0048,0,0.0001,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.0553,0.0623,0.0611,0.0612,0.0602,0.0587,0.0584,0.0594,0.0602,0.062,0.0669,0.0049,0,0.0001,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.0548,0.0618,0.0607,0.0606,0.0596,0.058,0.0578,0.0588,0.0596,0.0614,0.0664,0.005,0,0.0001,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.0543,0.0614,0.0603,0.0602,0.0591,0.0575,0.0573,0.0583,0.0592,0.061,0.0661,0.0051,0,0.0001,0,0,0,0,0,0,0)

ts1825<-c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0008,0.008,0.0461,0.2037,0.2818,0.3323,0.3652,0.3866,0.4001,0.4087,0.4141,0.4174,0.417,0.408,0.3219,0.2796,0.2366,0.1958,0.1609,0.1309,0.1027,0.0836,0.0689,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0008,0.008,0.0462,0.1898,0.2693,0.3207,0.3542,0.376,0.3898,0.3985,0.404,0.4074,0.4071,0.398,0.3221,0.2797,0.2368,0.196,0.1611,0.131,0.1029,0.0837,0.0691,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0008,0.0081,0.0463,0.1731,0.2541,0.3066,0.3408,0.363,0.3771,0.386,0.3917,0.3953,0.3953,0.3862,0.3228,0.2804,0.2375,0.1966,0.1617,0.1317,0.1035,0.0843,0.0697,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0008,0.0081,0.0459,0.1571,0.2394,0.2932,0.3282,0.3509,0.3653,0.3744,0.3804,0.3842,0.3843,0.3752,0.3241,0.2816,0.2389,0.1979,0.163,0.1329,0.1046,0.0853,0.0707,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0008,0.0081,0.0458,0.1554,0.2381,0.2922,0.3273,0.3502,0.3647,0.3737,0.3794,0.3821,0.3794,0.3652,0.327,0.2843,0.2414,0.2002,0.165,0.1349,0.1064,0.0871,0.0723,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0008,0.008,0.0457,0.1546,0.2375,0.2917,0.327,0.3499,0.3644,0.3734,0.3791,0.3818,0.3791,0.3649,0.3271,0.2844,0.2421,0.2034,0.1684,0.1376,0.1086,0.0892,0.0741,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0008,0.008,0.0457,0.154,0.237,0.2913,0.3267,0.3496,0.3641,0.3732,0.3789,0.3817,0.379,0.3648,0.3271,0.2845,0.2422,0.2035,0.1686,0.1381,0.1123,0.0918,0.0761,0.0525,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0008,0.008,0.0456,0.1535,0.2366,0.291,0.3264,0.3494,0.364,0.3731,0.3788,0.3816,0.3789,0.3647,0.3271,0.2846,0.2422,0.2035,0.1686,0.1382,0.1124,0.0922,0.0791,0.0744,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0007,0.008,0.0456,0.1531,0.2363,0.2908,0.3262,0.3493,0.3638,0.373,0.3787,0.3815,0.3788,0.3646,0.3272,0.2846,0.2423,0.2036,0.1687,0.1382,0.1125,0.0923,0.0792,0.0753,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0007,0.008,0.0455,0.1528,0.2361,0.2906,0.3261,0.3491,0.3637,0.3729,0.3786,0.3814,0.3788,0.3646,0.3272,0.2846,0.2423,0.2036,0.1687,0.1383,0.1125,0.0923,0.0792,0.0754,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0007,0.008,0.0455,0.1526,0.2359,0.2905,0.326,0.3491,0.3636,0.3728,0.3785,0.3813,0.3787,0.3646,0.3272,0.2847,0.2424,0.2036,0.1687,0.1383,0.1125,0.0923,0.0792,0.0754)

FE<-data.frame(x=X,z=Z,ts181=ts181,ts1825=ts1825)


X2<-c(50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950)

Z2<-c(49.5,49.4,49.3,49.2,49.1,49.0,48.9,48.8,48.7,48.6,48.5,48.4,48.3,48.2,48.1,48.0,47.9,47.8,47.7,47.6,47.5,47.4,47.3,47.2,47.1,47.0,46.9,46.8,46.7,46.6,46.5,46.4,46.3,46.2,46.1,46.0,45.9,45.8,45.7,45.6,48.5,48.4,48.3,48.2,48.1,48.0,47.9,47.8,47.7,47.6,47.5,47.4,47.3,47.2,47.1,47.0,46.9,46.8,46.7,46.6,46.5,46.4,46.3,46.2,46.1,46.0,45.9,45.8,45.7,45.6,45.5,45.4,45.3,45.2,45.1,45.0,44.9,44.8,44.7,44.6,47.5,47.4,47.3,47.2,47.1,47.0,46.9,46.8,46.7,46.6,46.5,46.4,46.3,46.2,46.1,46.0,45.9,45.8,45.7,45.6,45.5,45.4,45.3,45.2,45.1,45.0,44.9,44.8,44.7,44.6,44.5,44.4,44.3,44.2,44.1,44.0,43.9,43.8,43.7,43.6,46.5,46.4,46.3,46.2,46.1,46.0,45.9,45.8,45.7,45.6,45.5,45.4,45.3,45.2,45.1,45.0,44.9,44.8,44.7,44.6,44.5,44.4,44.3,44.2,44.1,44.0,43.9,43.8,43.7,43.6,43.5,43.4,43.3,43.2,43.1,43.0,42.9,42.8,42.7,42.6,45.5,45.4,45.3,45.2,45.1,45.0,44.9,44.8,44.7,44.6,44.5,44.4,44.3,44.2,44.1,44.0,43.9,43.8,43.7,43.6,43.5,43.4,43.3,43.2,43.1,43.0,42.9,42.8,42.7,42.6,42.5,42.4,42.3,42.2,42.1,42.0,41.9,41.8,41.7,41.6,44.5,44.4,44.3,44.2,44.1,44.0,43.9,43.8,43.7,43.6,43.5,43.4,43.3,43.2,43.1,43.0,42.9,42.8,42.7,42.6,42.5,42.4,42.3,42.2,42.1,42.0,41.9,41.8,41.7,41.6,41.5,41.4,41.3,41.2,41.1,41.0,40.9,40.8,40.7,40.6,43.5,43.4,43.3,43.2,43.1,43.0,42.9,42.8,42.7,42.6,42.5,42.4,42.3,42.2,42.1,42.0,41.9,41.8,41.7,41.6,41.5,41.4,41.3,41.2,41.1,41.0,40.9,40.8,40.7,40.6,40.5,40.4,40.3,40.2,40.1,40.0,39.9,39.8,39.7,39.6,42.5,42.4,42.3,42.2,42.1,42.0,41.9,41.8,41.7,41.6,41.5,41.4,41.3,41.2,41.1,41.0,40.9,40.8,40.7,40.6,40.5,40.4,40.3,40.2,40.1,40.0,39.9,39.8,39.7,39.6,39.5,39.4,39.3,39.2,39.1,39.0,38.9,38.8,38.7,38.6,41.5,41.4,41.3,41.2,41.1,41.0,40.9,40.8,40.7,40.6,40.5,40.4,40.3,40.2,40.1,40.0,39.9,39.8,39.7,39.6,39.5,39.4,39.3,39.2,39.1,39.0,38.9,38.8,38.7,38.6,38.5,38.4,38.3,38.2,38.1,38.0,37.9,37.8,37.7,37.6,40.5,40.4,40.3,40.2,40.1,40.0,39.9,39.8,39.7,39.6,39.5,39.4,39.3,39.2,39.1,39.0,38.9,38.8,38.7,38.6,38.5,38.4,38.3,38.2,38.1,38.0,37.9,37.8,37.7,37.6,37.5,37.4,37.3,37.2,37.1,37.0,36.9,36.8,36.7,36.6,39.5,39.4,39.3,39.2,39.1,39.0,38.9,38.8,38.7,38.6,38.5,38.4,38.3,38.2,38.1,38.0,37.9,37.8,37.7,37.6,37.5,37.4,37.3,37.2,37.1,37.0,36.9,36.8,36.7,36.6,36.5,36.4,36.3,36.2,36.1,36.0,35.9,35.8,35.7,35.6,38.5,38.4,38.3,38.2,38.1,38.0,37.9,37.8,37.7,37.6,37.5,37.4,37.3,37.2,37.1,37.0,36.9,36.8,36.7,36.6,36.5,36.4,36.3,36.2,36.1,36.0,35.9,35.8,35.7,35.6,35.5,35.4,35.3,35.2,35.1,35.0,34.9,34.8,34.7,34.6,37.5,37.4,37.3,37.2,37.1,37.0,36.9,36.8,36.7,36.6,36.5,36.4,36.3,36.2,36.1,36.0,35.9,35.8,35.7,35.6,35.5,35.4,35.3,35.2,35.1,35.0,34.9,34.8,34.7,34.6,34.5,34.4,34.3,34.2,34.1,34.0,33.9,33.8,33.7,33.6,36.5,36.4,36.3,36.2,36.1,36.0,35.9,35.8,35.7,35.6,35.5,35.4,35.3,35.2,35.1,35.0,34.9,34.8,34.7,34.6,34.5,34.4,34.3,34.2,34.1,34.0,33.9,33.8,33.7,33.6,33.5,33.4,33.3,33.2,33.1,33.0,32.9,32.8,32.7,32.6,35.5,35.4,35.3,35.2,35.1,35.0,34.9,34.8,34.7,34.6,34.5,34.4,34.3,34.2,34.1,34.0,33.9,33.8,33.7,33.6,33.5,33.4,33.3,33.2,33.1,33.0,32.9,32.8,32.7,32.6,32.5,32.4,32.3,32.2,32.1,32.0,31.9,31.8,31.7,31.6)

ts181.2<-c(0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.144,0.148,0.154,0.156,0.156,0.157,0.157,0.158,0.159,0.160,0.165,0.001,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.126,0.131,0.137,0.140,0.140,0.140,0.140,0.141,0.142,0.143,0.148,0.001,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.109,0.114,0.119,0.122,0.123,0.123,0.123,0.124,0.125,0.126,0.131,0.001,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.092,0.096,0.102,0.105,0.105,0.106,0.106,0.107,0.108,0.109,0.114,0.002,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.075,0.079,0.084,0.087,0.088,0.088,0.089,0.090,0.091,0.092,0.096,0.002,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.058,0.062,0.066,0.069,0.071,0.071,0.072,0.072,0.073,0.075,0.078,0.002,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.042,0.047,0.049,0.053,0.056,0.056,0.056,0.057,0.058,0.059,0.062,0.003,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.038,0.041,0.041,0.041,0.041,0.042,0.042,0.043,0.044,0.045,0.047,0.003,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.037,0.040,0.040,0.039,0.039,0.039,0.040,0.040,0.040,0.040,0.039,0.004,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.036,0.039,0.039,0.038,0.038,0.039,0.039,0.039,0.039,0.039,0.037,0.004,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.035,0.039,0.038,0.037,0.037,0.038,0.038,0.038,0.039,0.039,0.037,0.004,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.035,0.038,0.038,0.037,0.037,0.037,0.038,0.038,0.038,0.038,0.036,0.004,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.034,0.038,0.037,0.036,0.036,0.037,0.037,0.037,0.038,0.038,0.036,0.004,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.034,0.037,0.037,0.036,0.036,0.037,0.037,0.037,0.037,0.037,0.036,0.004,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.034,0.037,0.037,0.036,0.036,0.036,0.037,0.037,0.037,0.037,0.035,0.004,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000)

ts1825.2<-c(0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.008,0.037,0.191,0.264,0.313,0.346,0.368,0.384,0.401,0.407,0.411,0.411,0.403,0.279,0.230,0.189,0.156,0.126,0.102,0.079,0.051,0.054,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.038,0.175,0.249,0.299,0.333,0.356,0.373,0.389,0.395,0.400,0.400,0.391,0.285,0.236,0.193,0.160,0.129,0.104,0.081,0.053,0.055,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.038,0.161,0.236,0.287,0.322,0.345,0.362,0.379,0.385,0.389,0.389,0.381,0.290,0.240,0.197,0.163,0.132,0.107,0.083,0.054,0.055,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.149,0.225,0.276,0.311,0.335,0.352,0.369,0.375,0.380,0.380,0.371,0.293,0.243,0.200,0.165,0.134,0.109,0.084,0.055,0.056,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.148,0.224,0.275,0.310,0.333,0.348,0.358,0.365,0.370,0.370,0.361,0.297,0.247,0.203,0.169,0.137,0.112,0.086,0.057,0.056,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.147,0.223,0.275,0.310,0.333,0.348,0.358,0.365,0.369,0.369,0.354,0.302,0.252,0.208,0.173,0.141,0.115,0.089,0.059,0.057,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.147,0.223,0.275,0.309,0.332,0.348,0.358,0.365,0.369,0.368,0.354,0.302,0.253,0.211,0.176,0.144,0.118,0.092,0.060,0.058,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.147,0.223,0.274,0.309,0.332,0.348,0.357,0.365,0.369,0.368,0.353,0.302,0.253,0.211,0.176,0.146,0.120,0.094,0.062,0.058,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.146,0.223,0.274,0.309,0.332,0.347,0.357,0.364,0.369,0.368,0.353,0.302,0.253,0.211,0.176,0.147,0.121,0.097,0.064,0.055,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.146,0.222,0.274,0.309,0.332,0.347,0.357,0.364,0.369,0.368,0.353,0.302,0.253,0.211,0.176,0.147,0.121,0.097,0.066,0.046,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.146,0.222,0.274,0.308,0.332,0.347,0.357,0.364,0.369,0.368,0.353,0.302,0.253,0.211,0.176,0.147,0.121,0.097,0.066,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.146,0.222,0.274,0.308,0.332,0.347,0.357,0.364,0.368,0.368,0.353,0.302,0.253,0.211,0.176,0.147,0.121,0.097,0.066,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.146,0.222,0.273,0.308,0.332,0.347,0.357,0.364,0.368,0.368,0.353,0.302,0.253,0.211,0.176,0.147,0.121,0.097,0.066,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.146,0.222,0.273,0.308,0.331,0.347,0.357,0.364,0.368,0.368,0.353,0.302,0.253,0.211,0.176,0.147,0.121,0.097,0.066,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.145,0.222,0.273,0.308,0.331,0.347,0.357,0.364,0.368,0.368,0.353,0.302,0.253,0.211,0.176,0.147,0.121,0.097,0.066,0.000)
FD<-data.frame(x=X2,z=Z2,ts181=ts181.2,ts1825=ts1825.2)

plot(FE$x,FE$z,xlim=range(FE$x),ylim=c(0.8*range(FE$z)[1],1.1*range(FE$z)[2]))
points(FD$x,FD$z,pch=4)
legend("bottomleft",c("Finite Element Mesh
Nodes","Finite Difference Grid
Cell Nodes"),pch=c(1,4))
 


--
View this message in context:
http://r.789695.n4.nabble.com/contouring-x-y-scatter-data-tp3912828p3912828.html
Sent from the R help mailing list archive at Nabble.com.

Carlos Ortega

2011-Oct-17 21:37 UTC

head link

[R] contouring x y scatter data

Hi,

Please check if this representation helps you.

library(lattice)
library(latticeExtra)
p1.gr <- cloud( ts181 ~ x * z, data=FD, panel.3d.cloud = panel.3dbars,
scales = list(arrows = FALSE, just = "right"))
p2.gr <- cloud( ts1825 ~ x * z, data=FD, panel.3d.cloud = panel.3dbars,
scales = list(arrows = FALSE, just = "right"))
print(p1.gr, position = c(0, 0, 0.5, 1), more = TRUE)
print(p2.gr, position = c(0.5, 0, 1, 1))

Regards,
Carlos Ortega
www.qualityexcellence.es

2011/10/17 emorway <emorway@usgs.gov>
> Hello,
>
> I'm almost positive R can do the following, I just haven't hit upon
the
> right package or search terms, however.  Here's what I'm after: 
I've got
> concentration output from two different models that I want to qualitatively
> compare in a contour plot (or some variant of a contour plot).  The problem
> as I see it is that the data is not gridded is the usual regular fashion,
> and even if it were, the models are working with two different irregular
> computational/numerical grids at which concentrations are calculated.
> Running the code below will produce a plot that will hopefully illuminate
> what I'm attempting to describe.  At each location where 'O'
appears, there
> is a calculated concentration that I would like to contour.  I'd like
to do
> the same thing at each 'X' location and then compare, visually, how
close
> the contours are (hopefully close).  Both data.frames contain
> columns/fields
> titled "ts181" and "ts1825" which contain the
concentrations at time step
> 181 days and time step 1825 days.  Is it possible to contour x-y-value data
> if it is not laid out on a regular grid?
>
>
>
>
X<-c(0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,5!
>
> 
00,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,0,100,200,300,400,500,600,700,800,900,1000!
>  ,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200,2300,2400
>
>
,2500,2600,2700,2800,2900,3000,3100,3200,3300,3400,3500,3600,3700,3800,3900,4000)
>
>
>
Z<-c(50,49.95,49.85,49.75,49.65,49.55,49.45,49.35,49.25,49.15,49.05,48.95,48.85,48.75,48.65,48.55,48.45,48.35,48.25,48.15,48.05,47.95,47.85,47.75,47.65,47.55,47.45,47.35,47.25,47.15,47.05,46.95,46.85,46.75,46.65,46.55,46.45,46.35,46.25,46.15,46.1,49.25,49.2,49.1,49,48.9,48.8,48.7,48.6,48.5,48.4,48.3,48.2,48.1,48,47.9,47.8,47.7,47.6,47.5,47.4,47.3,47.2,47.1,47,46.9,46.8,46.7,46.6,46.5,46.4,46.3,46.2,46.1,46,45.9,45.8,45.7,45.6,45.5,45.4,45.35,48.2,48.15,48.05,47.95,47.85,47.75,47.65,47.55,47.45,47.35,47.25,47.15,47.05,46.95,46.85,46.75,46.65,46.55,46.45,46.35,46.25,46.15,46.05,45.95,45.85,45.75,45.65,45.55,45.45,45.35,45.25,45.15,45.05,44.95,44.85,44.75,44.65,44.55,44.45,44.35,44.3,47,46.95,46.85,46.75,46.65,46.55,46.45,46.35,46.25,46.15,46.05,45.95,45.85,45.75,45.65,45.55,45.45,45.35,45.25,45.15,45.05,44.95,44.85,44.75,44.65,44.55,44.45,44.35,44.25,44.15,44.05,43.95,43.85,43.75,43.65,43.55,43.45,43.35,43.25,43.15,43.1,45.5,45.45,45.35,45.25,45.15,45.05,44.95,44.85,44.75,44.6!
>
> 
5,44.55,44.45,44.35,44.25,44.15,44.05,43.95,43.85,43.75,43.65,43.55,43.45,43.35,43.25,43.15,43.05,42.95,42.85,42.75,42.65,42.55,42.45,42.35,42.25,42.15,42.05,41.95,41.85,41.75,41.65,41.6,44,43.95,43.85,43.75,43.65,43.55,43.45,43.35,43.25,43.15,43.05,42.95,42.85,42.75,42.65,42.55,42.45,42.35,42.25,42.15,42.05,41.95,41.85,41.75,41.65,41.55,41.45,41.35,41.25,41.15,41.05,40.95,40.85,40.75,40.65,40.55,40.45,40.35,40.25,40.15,40.1,42.5,42.45,42.35,42.25,42.15,42.05,41.95,41.85,41.75,41.65,41.55,41.45,41.35,41.25,41.15,41.05,40.95,40.85,40.75,40.65,40.55,40.45,40.35,40.25,40.15,40.05,39.95,39.85,39.75,39.65,39.55,39.45,39.35,39.25,39.15,39.05,38.95,38.85,38.75,38.65,38.6,41,40.95,40.85,40.75,40.65,40.55,40.45,40.35,40.25,40.15,40.05,39.95,39.85,39.75,39.65,39.55,39.45,39.35,39.25,39.15,39.05,38.95,38.85,38.75,38.65,38.55,38.45,38.35,38.25,38.15,38.05,37.95,37.85,37.75,37.65,37.55,37.45,37.35,37.25,37.15,37.1,39.5,39.45,39.35,39.25,39.15,39.05,38.95,38.85,38.75,38.65,38.55,38.45,38!
>  .35,38.25,38.15,38.05,37.95,37.85,37.75,37.65,37.55,37.45,37.35,37.25,
>
>
37.15,37.05,36.95,36.85,36.75,36.65,36.55,36.45,36.35,36.25,36.15,36.05,35.95,35.85,35.75,35.65,35.6,38,37.95,37.85,37.75,37.65,37.55,37.45,37.35,37.25,37.15,37.05,36.95,36.85,36.75,36.65,36.55,36.45,36.35,36.25,36.15,36.05,35.95,35.85,35.75,35.65,35.55,35.45,35.35,35.25,35.15,35.05,34.95,34.85,34.75,34.65,34.55,34.45,34.35,34.25,34.15,34.1,35,34.95,34.85,34.75,34.65,34.55,34.45,34.35,34.25,34.15,34.05,33.95,33.85,33.75,33.65,33.55,33.45,33.35,33.25,33.15,33.05,32.95,32.85,32.75,32.65,32.55,32.45,32.35,32.25,32.15,32.05,31.95,31.85,31.75,31.65,31.55,31.45,31.35,31.25,31.15,31.1)
>
>
>
ts181<-c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.113,0.1177,0.1167,0.1192,0.1242,0.1297,0.1317,0.131,0.1296,0.1296,0.1334,0.0014,0,0.0004,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.1074,0.1121,0.1111,0.1137,0.1187,0.1242,0.1262,0.1256,0.1241,0.1241,0.1279,0.0014,0,0.0004,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.0991,0.1039,0.1029,0.1053,0.1103,0.1159,0.1181,0.1175,0.116,0.116,0.1198,0.0015,0,0.0003,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.089,0.0939,0.0928,0.0951,0.1,0.1058,0.1082,0.1077,0.1063,0.1062,0.1101,0.0017,0,0.0003,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.0759,0.0811,0.0798,0.0816,0.0862,0.0923,0.0952,0.095,0.0936,0.0935,0.0973,0.0022,0,0.0003,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.0594,0.0659,0.0651,0.0674,0.0721,0.0785,0.0823,0.0825,0.0811,0.081,0.0846,0.0031,0,0.0002,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.0571,0.0639,0.0628,0.0632,0.062!
>
> 
8,0.0616,0.0613,0.0622,0.0629,0.0646,0.0693,0.0046,0,0.0002,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.0561,0.063,0.0618,0.0619,0.061,0.0596,0.0592,0.0602,0.0609,0.0627,0.0676,0.0048,0,0.0001,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.0553,0.0623,0.0611,0.0612,0.0602,0.0587,0.0584,0.0594,0.0602,0.062,0.0669,0.0049,0,0.0001,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.0548,0.0618,0.0607,0.0606,0.0596,0.058,0.0578,0.0588,0.0596,0.0614,0.0664,0.005,0,0.0001,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0002,0,0,0.0543,0.0614,0.0603,0.0602,0.0591,0.0575,0.0573,0.0583,0.0592,0.061,0.0661,0.0051,0,0.0001,0,0,0,0,0,0,0)
>
>
>
ts1825<-c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0008,0.008,0.0461,0.2037,0.2818,0.3323,0.3652,0.3866,0.4001,0.4087,0.4141,0.4174,0.417,0.408,0.3219,0.2796,0.2366,0.1958,0.1609,0.1309,0.1027,0.0836,0.0689,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0008,0.008,0.0462,0.1898,0.2693,0.3207,0.3542,0.376,0.3898,0.3985,0.404,0.4074,0.4071,0.398,0.3221,0.2797,0.2368,0.196,0.1611,0.131,0.1029,0.0837,0.0691,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0008,0.0081,0.0463,0.1731,0.2541,0.3066,0.3408,0.363,0.3771,0.386,0.3917,0.3953,0.3953,0.3862,0.3228,0.2804,0.2375,0.1966,0.1617,0.1317,0.1035,0.0843,0.0697,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0008,0.0081,0.0459,0.1571,0.2394,0.2932,0.3282,0.3509,0.3653,0.3744,0.3804,0.3842,0.3843,0.3752,0.3241,0.2816,0.2389,0.1979,0.163,0.1329,0.1046,0.0853,0.0707,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0008,0.0081,0.0458,0.1554,0.2381,0.2922,0.3273,0.3502,0.3647,0.3737,0.3794,0.3821,0.3794,0.3652,0.327,0.2843,0.2414,0.2002,0.165,0.1349,0.1064,0.0871,0.0723,0,0,0,0,0,!
>
> 
0,0,0,0,0,0,0,0,0,0,0,0,0,0.0008,0.008,0.0457,0.1546,0.2375,0.2917,0.327,0.3499,0.3644,0.3734,0.3791,0.3818,0.3791,0.3649,0.3271,0.2844,0.2421,0.2034,0.1684,0.1376,0.1086,0.0892,0.0741,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0008,0.008,0.0457,0.154,0.237,0.2913,0.3267,0.3496,0.3641,0.3732,0.3789,0.3817,0.379,0.3648,0.3271,0.2845,0.2422,0.2035,0.1686,0.1381,0.1123,0.0918,0.0761,0.0525,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0008,0.008,0.0456,0.1535,0.2366,0.291,0.3264,0.3494,0.364,0.3731,0.3788,0.3816,0.3789,0.3647,0.3271,0.2846,0.2422,0.2035,0.1686,0.1382,0.1124,0.0922,0.0791,0.0744,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0007,0.008,0.0456,0.1531,0.2363,0.2908,0.3262,0.3493,0.3638,0.373,0.3787,0.3815,0.3788,0.3646,0.3272,0.2846,0.2423,0.2036,0.1687,0.1382,0.1125,0.0923,0.0792,0.0753,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0007,0.008,0.0455,0.1528,0.2361,0.2906,0.3261,0.3491,0.3637,0.3729,0.3786,0.3814,0.3788,0.3646,0.3272,0.2846,0.2423,0.2036,0.1687,0.1383,0.1125,0.0923,0.0792,0.0754,0,0!
>  ,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.0007,0.008,0.0455,0.1526,0.2359,0.290
>
>
5,0.326,0.3491,0.3636,0.3728,0.3785,0.3813,0.3787,0.3646,0.3272,0.2847,0.2424,0.2036,0.1687,0.1383,0.1125,0.0923,0.0792,0.0754)
>
> FE<-data.frame(x=X,z=Z,ts181=ts181,ts1825=ts1825)
>
>
>
>
X2<-c(50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950!
>
> 
,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,18!
>  50,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,32
>
>
50,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950,50,150,250,350,450,550,650,750,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1850,1950,2050,2150,2250,2350,2450,2550,2650,2750,2850,2950,3050,3150,3250,3350,3450,3550,3650,3750,3850,3950)
>
>
>
Z2<-c(49.5,49.4,49.3,49.2,49.1,49.0,48.9,48.8,48.7,48.6,48.5,48.4,48.3,48.2,48.1,48.0,47.9,47.8,47.7,47.6,47.5,47.4,47.3,47.2,47.1,47.0,46.9,46.8,46.7,46.6,46.5,46.4,46.3,46.2,46.1,46.0,45.9,45.8,45.7,45.6,48.5,48.4,48.3,48.2,48.1,48.0,47.9,47.8,47.7,47.6,47.5,47.4,47.3,47.2,47.1,47.0,46.9,46.8,46.7,46.6,46.5,46.4,46.3,46.2,46.1,46.0,45.9,45.8,45.7,45.6,45.5,45.4,45.3,45.2,45.1,45.0,44.9,44.8,44.7,44.6,47.5,47.4,47.3,47.2,47.1,47.0,46.9,46.8,46.7,46.6,46.5,46.4,46.3,46.2,46.1,46.0,45.9,45.8,45.7,45.6,45.5,45.4,45.3,45.2,45.1,45.0,44.9,44.8,44.7,44.6,44.5,44.4,44.3,44.2,44.1,44.0,43.9,43.8,43.7,43.6,46.5,46.4,46.3,46.2,46.1,46.0,45.9,45.8,45.7,45.6,45.5,45.4,45.3,45.2,45.1,45.0,44.9,44.8,44.7,44.6,44.5,44.4,44.3,44.2,44.1,44.0,43.9,43.8,43.7,43.6,43.5,43.4,43.3,43.2,43.1,43.0,42.9,42.8,42.7,42.6,45.5,45.4,45.3,45.2,45.1,45.0,44.9,44.8,44.7,44.6,44.5,44.4,44.3,44.2,44.1,44.0,43.9,43.8,43.7,43.6,43.5,43.4,43.3,43.2,43.1,43.0,42.9,42.8,42.7,42.6,42.5,42.4,42.3,42.2,42.1,42.0,41.!
>
> 
9,41.8,41.7,41.6,44.5,44.4,44.3,44.2,44.1,44.0,43.9,43.8,43.7,43.6,43.5,43.4,43.3,43.2,43.1,43.0,42.9,42.8,42.7,42.6,42.5,42.4,42.3,42.2,42.1,42.0,41.9,41.8,41.7,41.6,41.5,41.4,41.3,41.2,41.1,41.0,40.9,40.8,40.7,40.6,43.5,43.4,43.3,43.2,43.1,43.0,42.9,42.8,42.7,42.6,42.5,42.4,42.3,42.2,42.1,42.0,41.9,41.8,41.7,41.6,41.5,41.4,41.3,41.2,41.1,41.0,40.9,40.8,40.7,40.6,40.5,40.4,40.3,40.2,40.1,40.0,39.9,39.8,39.7,39.6,42.5,42.4,42.3,42.2,42.1,42.0,41.9,41.8,41.7,41.6,41.5,41.4,41.3,41.2,41.1,41.0,40.9,40.8,40.7,40.6,40.5,40.4,40.3,40.2,40.1,40.0,39.9,39.8,39.7,39.6,39.5,39.4,39.3,39.2,39.1,39.0,38.9,38.8,38.7,38.6,41.5,41.4,41.3,41.2,41.1,41.0,40.9,40.8,40.7,40.6,40.5,40.4,40.3,40.2,40.1,40.0,39.9,39.8,39.7,39.6,39.5,39.4,39.3,39.2,39.1,39.0,38.9,38.8,38.7,38.6,38.5,38.4,38.3,38.2,38.1,38.0,37.9,37.8,37.7,37.6,40.5,40.4,40.3,40.2,40.1,40.0,39.9,39.8,39.7,39.6,39.5,39.4,39.3,39.2,39.1,39.0,38.9,38.8,38.7,38.6,38.5,38.4,38.3,38.2,38.1,38.0,37.9,37.8,37.7,37.6,37.5,37.4,37.3,37.2,3!
>  7.1,37.0,36.9,36.8,36.7,36.6,39.5,39.4,39.3,39.2,39.1,39.0,38.9,38.8,3
>
>
8.7,38.6,38.5,38.4,38.3,38.2,38.1,38.0,37.9,37.8,37.7,37.6,37.5,37.4,37.3,37.2,37.1,37.0,36.9,36.8,36.7,36.6,36.5,36.4,36.3,36.2,36.1,36.0,35.9,35.8,35.7,35.6,38.5,38.4,38.3,38.2,38.1,38.0,37.9,37.8,37.7,37.6,37.5,37.4,37.3,37.2,37.1,37.0,36.9,36.8,36.7,36.6,36.5,36.4,36.3,36.2,36.1,36.0,35.9,35.8,35.7,35.6,35.5,35.4,35.3,35.2,35.1,35.0,34.9,34.8,34.7,34.6,37.5,37.4,37.3,37.2,37.1,37.0,36.9,36.8,36.7,36.6,36.5,36.4,36.3,36.2,36.1,36.0,35.9,35.8,35.7,35.6,35.5,35.4,35.3,35.2,35.1,35.0,34.9,34.8,34.7,34.6,34.5,34.4,34.3,34.2,34.1,34.0,33.9,33.8,33.7,33.6,36.5,36.4,36.3,36.2,36.1,36.0,35.9,35.8,35.7,35.6,35.5,35.4,35.3,35.2,35.1,35.0,34.9,34.8,34.7,34.6,34.5,34.4,34.3,34.2,34.1,34.0,33.9,33.8,33.7,33.6,33.5,33.4,33.3,33.2,33.1,33.0,32.9,32.8,32.7,32.6,35.5,35.4,35.3,35.2,35.1,35.0,34.9,34.8,34.7,34.6,34.5,34.4,34.3,34.2,34.1,34.0,33.9,33.8,33.7,33.6,33.5,33.4,33.3,33.2,33.1,33.0,32.9,32.8,32.7,32.6,32.5,32.4,32.3,32.2,32.1,32.0,31.9,31.8,31.7,31.6)
>
>
>
ts181.2<-c(0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.144,0.148,0.154,0.156,0.156,0.157,0.157,0.158,0.159,0.160,0.165,0.001,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.126,0.131,0.137,0.140,0.140,0.140,0.140,0.141,0.142,0.143,0.148,0.001,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.109,0.114,0.119,0.122,0.123,0.123,0.123,0.124,0.125,0.126,0.131,0.001,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.092,0.096,0.102,0.105,0.105,0.106,0.106,0.107,0.108,0.109,0.114,0.002,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,!
>
> 
0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.075,0.079,0.084,0.087,0.088,0.088,0.089,0.090,0.091,0.092,0.096,0.002,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.058,0.062,0.066,0.069,0.071,0.071,0.072,0.072,0.073,0.075,0.078,0.002,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.042,0.047,0.049,0.053,0.056,0.056,0.056,0.057,0.058,0.059,0.062,0.003,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.038,0.041,0.041,0.041,0.041,0.042,0.042,0.043,0.044,0.045,0.047,0.003,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.00!
>  0,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.
>
>
001,0.037,0.040,0.040,0.039,0.039,0.039,0.040,0.040,0.040,0.040,0.039,0.004,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.036,0.039,0.039,0.038,0.038,0.039,0.039,0.039,0.039,0.039,0.037,0.004,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.035,0.039,0.038,0.037,0.037,0.038,0.038,0.038,0.039,0.039,0.037,0.004,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.035,0.038,0.038,0.037,0.037,0.037,0.038,0.038,0.038,0.038,0.036,0.004,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.034,0.038,0.037,0.036,0!
>
> 
.036,0.037,0.037,0.037,0.038,0.038,0.036,0.004,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.034,0.037,0.037,0.036,0.036,0.037,0.037,0.037,0.037,0.037,0.036,0.004,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.001,0.034,0.037,0.037,0.036,0.036,0.036,0.037,0.037,0.037,0.037,0.035,0.004,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000)
>
>
>
ts1825.2<-c(0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.008,0.037,0.191,0.264,0.313,0.346,0.368,0.384,0.401,0.407,0.411,0.411,0.403,0.279,0.230,0.189,0.156,0.126,0.102,0.079,0.051,0.054,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.038,0.175,0.249,0.299,0.333,0.356,0.373,0.389,0.395,0.400,0.400,0.391,0.285,0.236,0.193,0.160,0.129,0.104,0.081,0.053,0.055,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.038,0.161,0.236,0.287,0.322,0.345,0.362,0.379,0.385,0.389,0.389,0.381,0.290,0.240,0.197,0.163,0.132,0.107,0.083,0.054,0.055,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.149,0.225,0.276,0.311,0.335,0.352,0.369,0.375,0.380,0.380,0.371,0.293,0.243,0.200,0.165,0.134,0.109,0.084,0.055,0.056,0.000,0.000,0.000!
>
> 
,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.148,0.224,0.275,0.310,0.333,0.348,0.358,0.365,0.370,0.370,0.361,0.297,0.247,0.203,0.169,0.137,0.112,0.086,0.057,0.056,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.147,0.223,0.275,0.310,0.333,0.348,0.358,0.365,0.369,0.369,0.354,0.302,0.252,0.208,0.173,0.141,0.115,0.089,0.059,0.057,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.147,0.223,0.275,0.309,0.332,0.348,0.358,0.365,0.369,0.368,0.354,0.302,0.253,0.211,0.176,0.144,0.118,0.092,0.060,0.058,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.147,0.223,0.274,0.309,0.332,0.348,0.357,0.365,0.369,0.368,0.353,0.302,0.253,0.211,0.176,0.146,0.120,0.094,0.062,0.058,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.0!
>  00,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0
>
>
.039,0.146,0.223,0.274,0.309,0.332,0.347,0.357,0.364,0.369,0.368,0.353,0.302,0.253,0.211,0.176,0.147,0.121,0.097,0.064,0.055,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.146,0.222,0.274,0.309,0.332,0.347,0.357,0.364,0.369,0.368,0.353,0.302,0.253,0.211,0.176,0.147,0.121,0.097,0.066,0.046,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.146,0.222,0.274,0.308,0.332,0.347,0.357,0.364,0.369,0.368,0.353,0.302,0.253,0.211,0.176,0.147,0.121,0.097,0.066,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.146,0.222,0.274,0.308,0.332,0.347,0.357,0.364,0.368,0.368,0.353,0.302,0.253,0.211,0.176,0.147,0.121,0.097,0.066,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.146,0.222,0.273,0.308,!
>
> 
0.332,0.347,0.357,0.364,0.368,0.368,0.353,0.302,0.253,0.211,0.176,0.147,0.121,0.097,0.066,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.146,0.222,0.273,0.308,0.331,0.347,0.357,0.364,0.368,0.368,0.353,0.302,0.253,0.211,0.176,0.147,0.121,0.097,0.066,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.002,0.009,0.039,0.145,0.222,0.273,0.308,0.331,0.347,0.357,0.364,0.368,0.368,0.353,0.302,0.253,0.211,0.176,0.147,0.121,0.097,0.066,0.000)
> FD<-data.frame(x=X2,z=Z2,ts181=ts181.2,ts1825=ts1825.2)
>
>
>
plot(FE$x,FE$z,xlim=range(FE$x),ylim=c(0.8*range(FE$z)[1],1.1*range(FE$z)[2]))
> points(FD$x,FD$z,pch=4)
> legend("bottomleft",c("Finite Element Mesh
Nodes","Finite Difference Grid
> Cell Nodes"),pch=c(1,4))
>
>
>
> --
> View this message in context:
>
http://r.789695.n4.nabble.com/contouring-x-y-scatter-data-tp3912828p3912828.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help@r-project.org mailing list
> 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.
>
	[[alternative HTML version deleted]]

Possibly Parallel Threads

Search for more reasonably related threads

R help - Oct 2011 - contouring x y scatter data

[R] contouring x y scatter data

[R] contouring x y scatter data

Possibly Parallel Threads