[R] a similar question

arun smartpink111 at yahoo.com
Wed Mar 27 00:02:01 CET 2013


HI Elisa,

You need to review the formulas.  Some of them (especially between 84 & 85 or 4 & 84 is not making sense, possibly typos).  I changed it according to the pattern of the formulas.
source("ElisaNew.txt")
#mat1 # dataset

res<-do.call(cbind,lapply(seq_len(nrow(mat1)),function(i) do.call(rbind,lapply(split(rbind(mat1[i,],mat1[-i,]),1:nrow(rbind(mat1[i,],mat1[-i,]))), function(x) {x1<- rbind(mat1[i,],x);x1<-as.data.frame(x1);
if({indx<-colSums(is.na(x1));indx[1]==0 & indx[2]==0 & indx[3]==1 & indx[4]==2}) #3 peaks 2 peaks comparison
 {x2<- x1[rev(order(x1$Peak3,x1$Npeak3)),];
with(x2,{abs(Peak1[1]-Peak1[2])+abs(Peak2[1]-Peak2[2])+abs(Peak1[1]-Peak3[2]) + abs(Npeak1[1]-Npeak1[2])+abs(Npeak2[1]-Npeak2[2])+abs(Npeak1[1]-Npeak3[2])})
}

else if({indx[1]==0 & indx[2]==1 & indx[3]==2 & indx[4]==2}) #2 peaks 1 peak comparison
{x3<- x1[rev(order(x1$Peak2,x1$Npeak2)),];
with(x3,{abs(Peak1[1]-Peak1[2])+ abs(Peak1[1]-Peak2[2])+ abs(Npeak1[1]-Npeak1[2])+ abs(Npeak1[1]-Npeak2[2])})
}
else if({indx[1]==0 & indx[2]==0 & indx[3]==1 & indx[4]==1}) #4 peaks 2 peaks comparison
{x4<- x1[rev(order(x1$Peak3,x1$Peak4,x1$Npeak3,x1$Npeak4)),];
with(x4,{abs(Peak1[1]-Peak1[2])+abs(Peak2[1]-Peak2[2])+ abs(Peak1[1]-Peak3[2])+ abs(Peak2[1]-Peak4[2])+ abs(Npeak1[1]-Npeak1[2])+abs(Npeak2[1]-Npeak2[2])+ abs(Npeak1[1]-Npeak3[2])+ abs(Npeak2[1]-Npeak4[2])})
}

else if({indx[1]==0 & indx[2]==0 & indx[3]==0 & indx[4]==1}) #4 peaks 3 peaks comparison
{x5<- x1[rev(order(x1$Peak4,x1$Npeak4)),];
with(x5,{abs(Peak1[1]-Peak1[2])+abs(Peak2[1]-Peak2[2])+ abs(Peak3[1]-Peak3[2])+ abs(Peak1[1]-Peak4[2])+ abs(Npeak1[1]-Npeak1[2])+abs(Npeak2[1]-Npeak2[2])+ abs(Npeak3[1]-Npeak3[2])+ abs(Npeak1[1]-Npeak4[2])})
}
else if({indx[1]==0 & indx[2]==1 & indx[3]==1 & indx[4]==1}) #4 peak 1 peak comparison
{x6<- x1[rev(order(x1$Peak2,x1$Peak3,x1$Peak4,x1$Npeak2,x1$Npeak3,x1$Npeak4)),];
with(x6,{abs(Peak1[1]-Peak1[2])+ abs(Peak1[1]-Peak2[2])+abs(Peak1[1]-Peak3[2])+ abs(Peak1[1]-Peak4[2])+ abs(Npeak1[1]-Npeak1[2])+ abs(Npeak1[1]-Npeak2[2])+abs(Npeak1[1]-Npeak3[2])+ abs(Npeak1[1]-Npeak4[2])})
}
else if({indx[1]==0 & indx[2]==1 & indx[3]==1 & indx[4]==2}) # 3 peaks 1 peak comparison
{x7<- x1[rev(order(x1$Peak2,x1$Peak3,x1$Npeak2,x1$Npeak3)),];
with(x7,{abs(Peak1[1]-Peak1[2])+abs(Peak1[1]-Peak2[2])+abs(Peak1[1]-Peak3[2])+ abs(Npeak1[1]-Npeak1[2])+abs(Npeak1[1]-Npeak2[2])+abs(Npeak1[1]-Npeak3[2])})
}
else {x1[is.na(x1)]<-0; x1; #cases where peaks are similar
with(x1,{abs(Peak1[1]-Peak1[2])+ abs(Peak2[1]-Peak2[2])+ abs(Peak3[1]-Peak3[2])+abs(Peak4[1]-Peak4[2]) + abs(Npeak1[1]-Npeak1[2])+ abs(Npeak2[1]-Npeak2[2])+ abs(Npeak3[1]-Npeak3[2])+abs(Npeak4[1]-Npeak4[2])})
}

}))))

res2<-do.call(cbind,lapply(seq_len(ncol(res)),function(i) c(c(tail(res[seq(1,i,1),i],-1),0),res[-c(1:i),i])))
row.names(res2)<-1:nrow(res2)
dim(res2)
#[1] 124 124
res2[1:5,1:5]
#       [,1]      [,2]      [,3]     [,4]      [,5]
#1   0.00000  65.59415  86.62556 407.9987 104.78294
#2  65.59415   0.00000  42.14256 307.3830  39.18879
#3  86.62556  42.14256   0.00000 314.4331  33.88839
#4 407.99871 307.38297 314.43309   0.0000 266.78887
#5 104.78294  39.18879  33.88839 266.7889   0.00000
   


A.K.



________________________________
From: eliza botto <eliza_botto at hotmail.com>
To: "smartpink111 at yahoo.com" <smartpink111 at yahoo.com> 
Sent: Tuesday, March 26, 2013 1:53 PM
Subject: RE: a similar question



Dear Arun,
Few day ago I asked a question about Subtracted the data in row in a special format. 
         Peak1     Peak2    Peak3    Peak4    Npeak1    Npeak2   Npeak3   Npeak4
  [1,] 193.08000 158.09978       NA       NA 115.16080  60.80880       NA       NA
  [2,] 162.84205 138.72360       NA       NA 102.91291  57.07668       NA       NA
  [3,] 171.14483 140.97639       NA       NA  88.44784  39.95476       NA       NA
  [4,]  85.67824  85.03097 91.48871       NA  66.20270  79.96853 57.34178       NA
  [5,] 151.95125 134.14709       NA       NA  90.26042  46.00769       NA       NA
[84,]  80.66790  78.34769 81.12842 80.83636  64.33888  74.85277 70.20992 73.83368
[85,] 163.67442        NA       NA       NA  34.83554        NA       NA       NA

The substraction should take place in the following way...
between station 1 and 2
=abs(1,peak1-2,peak1)+abs(1,peak2-2,peak2)+abs(1,Npeak1-2,Npeak2)+abs(1,peak2-2,Npeak2)
between station 1 and 4
=abs(1,peak1-4,peak1)+abs(1,peak2-4,peak2)+abs(1,peak1-4,peak3)+abs(1,Npeak1-2,Npeak2)+abs(1,peak2-2,Npeak2)+abs(1,Npeak1-4,Npeak3)
between 1 and 84
=abs(1,peak1-84,peak1)+abs(1,peak2-84,peak2)+abs(1,peak1-84,peak3)+abs(1,peak2-84,peak4)+abs(1,Npeak1-84,Npeak2)+abs(1,peak2-84,Npeak2)+abs(1,Npeak1-84,Npeak3)+abs(1,Npeak2-84,Npeak4)
between 4 and 84
=abs(4,peak1-84,peak1)+abs(4,peak2-84,peak2)+abs(4,peak3-84,peak3)+abs(4,peak1-84,peak4)+abs(4,Npeak1-84,Npeak2)+abs(4,peak2-84,Npeak2)+abs(4,Npeak3-84,Npeak3)+abs(4,Npeak1-84,Npeak4)+abs(4,Npeak2-84,Npeak4)
between 84 and 85
=abs(84,peak1-85,peak1)+abs(84,peak2-85,peak1)+abs(84,peak3-85,peak1)+abs(84,peak1-85,peak1)+abs(84,Npeak1-85,Npeak1)+abs(84,peak2-85,Npeak1)+abs(84,Npeak3-85,Npeak1)+abs(84,Npeak1-85,Npeak1)+abs(84,Npeak2-85,Npeak1)

Below the dput of the actual data which you can directly use in R. I have also attached the previous email in which you solved a similar problem.


> dput(t)
structure(c(193.08, 162.8420477, 171.1448253, 85.67824387, 151.951247, 
97.46414674, 98.41596866, 127.6501823, 104.8280078, 193.4970353, 
199.0847667, 176.8056056, 149.0130919, 203.7191973, 158.7761534, 
151.1471448, 173.3468917, 188.7535598, 163.6243114, 104.3552439, 
144.2004793, 252.1751352, 188.3153014, 243.4516946, 120.458769, 
133.7964864, 92.27249234, 88.29083763, 133.2562227, 234.0726695, 
146.3646113, 149.369954, 147.2320252, 121.7396289, 77.9141118, 
87.61102932, 98.70319262, 93.7890117, 175.5433822, 93.7512927, 
94.45820072, 94.37275588, 103.1373724, 150.1042785, 138.7761189, 
188.7021462, 181.7828322, 206.3456388, 91.71107659, 144.3045039, 
182.093652, 178.012344, 158.4503902, 193.3935179, 141.9339532, 
122.2903509, 124.924136, 165.6671209, 199.2754651, 257.9251498, 
143.5997757, 173.1754596, 199.841469, 192.026085, 253.5559367, 
171.3493161, 169.2122863, 179.3942947, 152.9254538, 144.2730434, 
150.1081282, 132.6742234, 122.7814038, 162.845, 130.1472071, 
124.6752462, 128.0815213, 134.8649078, 129.7350348, 130.4832835, 
133.9640353, 111.7047607, 107.1153957, 80.6679022, 163.6744174, 
176.751121, 113.6244811, 328.6653384, 178.7597623, 213.5203283, 
147.8495361, 129.0191864, 120.1288664, 122.9193881, 226.82416, 
176.4431185, 185.3799489, 256.6602768, 177.5516373, 173.920279, 
173.027047, 208.0799694, 287.5454583, 234.4848, 181.837518, 132.7681561, 
125.6683659, 135.3671454, 178.9265525, 122.1269353, 114.6421333, 
176.4445375, 165.7989654, 176.5906123, 133.8557001, 182.2715586, 
169.9412251, 102.8431087, 118.1335804, 119.1328332, 104.7254147, 
190.1609391, 167.3936217, 220.7415613, 158.0997778, 138.7235962, 
140.9763937, 85.0309682, 134.147085, 88.43381717, 86.43197359, 
NA, 117.0416937, NA, 205.3162583, NA, NA, NA, NA, NA, NA, NA, 
141.9852952, 86.38999489, 180.192557, NA, NA, NA, 96.04636759, 
121.4449107, 101.1689267, 100.9065793, 121.4365135, NA, 162.7953145, 
163.3919406, 165.0108464, NA, 82.75190372, 87.15608139, 92.19705073, 
100.0923214, 175.1546501, 91.80899546, 96.81049963, 96.21629331, 
98.27490488, 154.7203677, 160.9197376, NA, NA, NA, NA, 136.6377027, 
166.613283, NA, NA, NA, 123.814623, 111.5650536, 110.0346231, 
138.1385062, NA, 219.2386009, 144.7317462, NA, NA, NA, NA, 138.0290855, 
136.4628199, 148.8132978, 131.3641005, 140.0425449, 170.5778213, 
107.3884982, 103.3489554, 131.591, 106.1624695, 111.7056902, 
102.3693357, 109.8419872, 103.3258681, 109.5187098, 156.0512068, 
114.6485948, 112.0649662, 78.34769086, NA, 144.952144, 106.4101364, 
337.7452996, NA, NA, 133.5147875, 122.3758765, 129.618582, 134.3367399, 
184.1394062, 142.4989924, 150.8782566, 261.879342, 144.6833765, 
NA, NA, 189.5035968, 236.8890365, NA, NA, 136.9814538, 128.7529761, 
153.3453913, NA, 124.7495929, 132.0709614, NA, NA, NA, 123.7208288, 
167.2707171, 161.433516, 123.2784396, 98.89750185, 100.9949358, 
103.1478139, 198.1966297, 188.9813008, NA, NA, NA, NA, 91.48870805, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 81.12842099, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 80.83635787, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, 115.1608029, 102.9129105, 88.44783724, 
66.20269535, 90.26042297, 55.93820907, 53.02374046, 38.46943507, 
41.93951961, 54.0228969, 131.1574231, 43.77704671, 42.73510486, 
46.55967105, 39.90845944, 39.34973343, 39.561716, 40.44960504, 
67.83789663, 39.77542181, 58.76160903, 55.55681811, 55.19231796, 
52.00997323, 56.23502902, 51.53584282, 42.716094, 40.96713125, 
51.18528191, 56.52295992, 85.97792004, 86.10068421, 78.73743715, 
39.65016889, 55.65975623, 59.81778569, 65.35525897, 44.5505851, 
113.3965602, 46.91458173, 44.93247935, 44.75638566, 48.76291865, 
70.81354773, 70.92721813, 46.77599113, 39.93167964, 38.79052359, 
40.03252748, 50.83829694, 54.88977092, 43.30494123, 58.38927899, 
122.4550136, 64.94447616, 48.11149672, 48.18377043, 96.80220343, 
44.79825334, 155.0829943, 44.06588188, 57.11966214, 44.20848045, 
43.13966225, 53.16048482, 86.35957119, 90.76594376, 56.84816651, 
54.54970929, 63.04449199, 74.49043212, 50.64394155, 61.07315829, 
64.4603, 54.74091917, 60.97339885, 51.22996817, 59.94545546, 
52.04453505, 68.02628358, 48.12415548, 48.18616981, 37.5333715, 
64.33887715, 34.83553793, 75.55990757, 61.47790315, 212.5252993, 
46.15653236, 51.66814997, 51.54448708, 48.15796924, 52.58811705, 
50.338595, 134.6565855, 92.947763, 113.2892524, 184.4882583, 
76.79211793, 78.02186687, 100.2169775, 134.192963, 172.3710276, 
135.007646, 77.34181146, 53.02946565, 51.45063901, 57.85099998, 
54.8408258, 51.00771395, 46.71004622, 50.99536823, 49.2562012, 
57.08851145, 89.29683815, 117.1361854, 114.0466909, 43.0627624, 
47.24608742, 48.61695044, 46.09568072, 156.5393615, 158.7677842, 
48.24319222, 60.8088045, 57.07667784, 39.95476442, 79.9685275, 
46.00768579, 36.10225934, 37.84199415, NA, 51.61021371, NA, 55.35805342, 
NA, NA, NA, NA, NA, NA, NA, 64.89165235, 36.39697421, 80.23448401, 
NA, NA, NA, 39.64508533, 43.2842462, 40.19597469, 40.34891239, 
43.24592537, NA, 78.20795608, 84.26777086, 82.66617621, NA, 52.78491697, 
46.39238199, 45.9445694, 44.6364773, 68.43642232, 47.11145011, 
44.51102427, 47.3619957, 42.77557491, 73.06434007, 70.84859867, 
NA, NA, NA, NA, 46.45195742, 117.5935617, NA, NA, 55.98705693, 
42.50930456, 51.66444439, 53.17215605, 35.68116766, NA, 62.7451114, 
51.08682513, NA, NA, NA, NA, 47.7346021, 46.59835742, 48.2035047, 
45.10152968, 60.15167948, 63.06433866, 45.02669254, 43.90926208, 
43.5724, 43.52469388, 50.90621174, 43.68814516, 43.17605811, 
43.20115733, 45.53042434, 72.50807737, 54.52220614, 59.02903389, 
74.85277294, NA, 30.20761354, 47.86751982, 63.5109485, NA, NA, 
80.17260595, 67.52261705, 68.89230992, 68.5740699, 54.70130851, 
41.90772568, 54.16163286, 61.3124205, 51.9215854, NA, NA, 41.64250609, 
69.71743068, NA, NA, 69.62073248, 64.77743919, 67.85062936, NA, 
62.46540787, 62.22649962, NA, NA, NA, 55.837478, 58.36762505, 
58.49723944, 53.90150304, 41.30913732, 43.46321644, 45.92690062, 
57.38622077, 60.96372166, NA, NA, NA, NA, 57.34177914, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 70.2099179, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 73.83368391, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
NA, NA, NA, NA, NA, NA, NA), .Dim = c(124L, 8L), .Dimnames = list(
    NULL, c("Peak1", "Peak2", "Peak3", "Peak4", "Npeak1", "Npeak2", 
    "Npeak3", "Npeak4")))

Thanks in advance
Elisa

> Date: Wed, 27 Feb 2013 16:14:28 -0800
> From: smartpink111 at yahoo.com
> Subject: Re: a similar question
> To: eliza_botto at hotmail.com
> CC: r-help at r-project.org
> 
> Dear Elisa,
> 
> Try this:
> dat1<- read.table(text="
>  V1    V2
> 1 1.939064 1.791503
> 2   1.670409 1.508740
> 3   2.641346       NA
> 4   2.772702       NA
> 5   3.401987       NA
> 6   2.323481 1.965391
> 7   2.135438 1.724173
> 8   2.404821       NA
> 9   3.273616       NA
> 10  2.058073       NA
> 11  2.711462       NA
> 12  1.871119 1.599446
> 13  2.599939       NA
> 14  2.174406 1.831875
> 15  2.377755       NA
> 16  1.735191 1.713195
> 17  1.968154 1.836061
> 18  3.455566       NA
> 19  2.443566       NA
> 20  1.865691 1.511079
> 21  2.219164 1.872533
> 22  3.011284       NA
> 23  1.573508 1.312573
> 24  2.645591       NA
> 25  1.781628 1.563686
> 26  4.223059       NA
> 27  2.617404       NA
> 28  2.910262       NA
> 29  2.621719       NA
> 30  2.084573       NA
> 31  2.779583       NA
> 32  2.925445       NA
> 33  2.768464       NA
> ",sep="",header=TRUE)
>  mat1<-as.matrix(dat1)
>  attr(mat1,"dimnames")<- NULL
>  mat2<- head(mat1)
> #mat2
> resNew2<-do.call(cbind,lapply(seq_len(nrow(mat2)),function(i) do.call(rbind,lapply(split(rbind(mat2[i,],mat2[-i,]),1:nrow(rbind(mat2[i,],mat2[-i,]))),function(x) {x1<- rbind(mat2[i,],x); x2<-x1[order(x1[,2]),];ifelse(sum(is.na(x2[,2]))==0, sum(abs(diff(x2[,1])),abs(diff(x2[,2]))),ifelse(sum(is.na(x2[,2]))==1,sum(abs(diff(x2[,1])),abs(x2[1,2]-x2[2,1])),abs(x2[1,1]-x2[2,1]))) }))))
>  resNew3<- do.call(cbind,lapply(seq_len(ncol(resNew2)),function(i) c(c(tail(resNew2[seq(1,i,1),i],-1),0),resNew2[-c(1:i),i])))
>  attr(resNew3,"dimnames")<-NULL
>  resNew3
> #         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]
> #[1,] 0.000000 0.551418 1.552125 1.814837 3.073407 0.558305
> #[2,] 0.551418 0.000000 2.103543 2.366255 3.624825 1.109723
> #[3,] 1.552125 2.103543 0.000000 0.131356 0.760641 0.993820
> #[4,] 1.814837 2.366255 0.131356 0.000000 0.629285 1.256532
> #[5,] 3.073407 3.624825 0.760641 0.629285 0.000000 2.515102
> #[6,] 0.558305 1.109723 0.993820 1.256532 2.515102 0.000000
> 
> 
> #mat1
> res<-do.call(cbind,lapply(seq_len(nrow(mat1)),function(i) do.call(rbind,lapply(split(rbind(mat1[i,],mat1[-i,]),1:nrow(rbind(mat1[i,],mat1[-i,]))),function(x) {x1<- rbind(mat1[i,],x); x2<-x1[order(x1[,2]),];ifelse(sum(is.na(x2[,2]))==0, sum(abs(diff(x2[,1])),abs(diff(x2[,2]))),ifelse(sum(is.na(x2[,2]))==1,sum(abs(diff(x2[,1])),abs(x2[1,2]-x2[2,1])),abs(x2[1,1]-x2[2,1]))) }))))
>  res2<- do.call(cbind,lapply(seq_len(ncol(res)),function(i) c(c(tail(res[seq(1,i,1),i],-1),0),res[-c(1:i),i])))
>  attr(res2,"dimnames")<-NULL
>  head(res2,3)
> #         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
> #[1,] 0.000000 0.551418 1.552125 1.814837 3.073407 0.558305 0.263704 1.079075
> #[2,] 0.551418 0.000000 2.103543 2.366255 3.624825 1.109723 0.680462 1.630493
> #[3,] 1.552125 2.103543 0.000000 0.131356 0.760641 0.993820 1.423081 0.236525
> #         [,9]    [,10]    [,11]    [,12]    [,13]    [,14]    [,15]    [,16]
> #[1,] 2.816665 0.385579 1.692357 0.260002 1.469311 0.275714 1.024943 0.282181
> #[2,] 3.368083 0.936997 2.243775 0.291416 2.020729 0.827132 1.576361 0.269237
> #[3,] 0.632270 0.583273 0.070116 1.812127 0.041407 1.276411 0.263591 1.834306
> #        [,17]    [,18]    [,19]    [,20]    [,21]    [,22]    [,23]    [,24]
> #[1,] 0.073648 3.180565 1.156565 0.353797 0.361130 2.292001 0.844486 1.560615
> #[2,] 0.625066 3.731983 1.707983 0.197621 0.912548 2.843419 0.293068 2.112033
> #[3,] 1.478477 0.814220 0.197780 1.905922 1.190995 0.369938 2.396611 0.004245
> #        [,25]    [,26]    [,27]    [,28]    [,29]    [,30]    [,31]    [,32]
> #[1,] 0.385253 4.715551 1.504241 2.089957 1.512871 0.438579 1.828599 2.120323
> #[2,] 0.166165 5.266969 2.055659 2.641375 2.064289 0.989997 2.380017 2.671741
> #[3,] 1.937378 1.581713 0.023942 0.268916 0.019627 0.556773 0.138237 0.284099
> #        [,33]
> #[1,] 1.806361
> #[2,] 2.357779
> #[3,] 0.127118
> 
> A.K.
> 
> 
> 
> 
> 
> 
> 
> 
> 
> ________________________________
> From: eliza botto <eliza_botto at hotmail.com>
> To: "smartpink111 at yahoo.com" <smartpink111 at yahoo.com> 
> Sent: Wednesday, February 27, 2013 6:49 PM
> Subject: RE: a similar question
> 
> 
> 
> Dear Arun,
> Sorry for late replying. i only want to be applied for two column as in case of "res" matrix i sent to you.
> Here it is again
>  [,1]       [,2]
>  1   1.939064 1.791503
>  2   1.670409 1.508740
>  3   2.641346       NA
>  4   2.772702       NA
>  5   3.401987       NA
>  6   2.323481 1.965391
>  7   2.135438 1.724173
>  8   2.404821       NA
>  9   3.273616       NA
>  10  2.058073       NA
>  11  2.711462       NA
>  12  1.871119 1.599446
>  13  2.599939       NA
>  14  2.174406 1.831875
>  15  2.377755       NA
>  16  1.735191 1.713195
>  17  1.968154 1.836061
>  18  3.455566       NA
>  19  2.443566       NA
>  20  1.865691 1.511079
>  21  2.219164 1.872533
>  22  3.011284       NA
>  23  1.573508 1.312573
>  24  2.645591       NA
>  25  1.781628 1.563686
>  26  4.223059       NA
>  27  2.617404       NA
>  28  2.910262       NA
>  29  2.621719       NA
>  30  2.084573       NA
>  31  2.779583       NA
>  32  2.925445       NA
>  33  2.768464       NA
> 
> 
> > Date: Wed, 27 Feb 2013 14:00:09 -0800
> > From: smartpink111 at yahoo.com
> > Subject: Re: a similar question
> > To: eliza_botto at hotmail.com
> > 
> > HI Elisa,
> > 
> > Your question is not clear.  
> > 
> >  dim(dat) #dput data
> > #[1]  33 124
> > 
> > 
> > 
> > Between row 1 and 2
> > ===============
> > abs[(1,1)-(2,1)]+abs[(1,2)-(2,2)]
> > 
> > Between row 1 and 3
> > ===============
> > abs[(1,1)-(3,1)]+abs[(1,2)-(3,1)]
> > 
> > You saw that when there is an NA in the row then both the values are to subtracted from the same number.
> > 
> > Now if there is NA in both the rows then
> > 
> > Between row 3 and 4
> > ===============
> > abs[(3,1)-(4,1)]
> > 
> > Here, there are many columns and how do you want this to be applied?
> > For example in the simple case of 2 columns, this is possible, but with multiple columns and multiple rows with NA or numbers, how is the 2nd and 3rd possibility to be applied?
> > Arun
> > 
> > 
> > 
> > 
> > 
> > 
> > 
> > 
> > 
> > ________________________________
> > From: eliza botto <eliza_botto at hotmail.com>
> > To: "smartpink111 at yahoo.com" <smartpink111 at yahoo.com> 
> > Sent: Wednesday, February 27, 2013 12:48 PM
> > Subject: a similar question
> > 
> > 
> > 
> > Dear Arun,
> > 
> > I have a much problem to the one i sent you earlier...
> > 
> >  i have two columns of a matrix with 124 rows like the followings
> > >res
> >    [,1]       [,2]
> > 1   1.939064 1.791503
> > 2   1.670409 1.508740
> > 3   2.641346       NA
> > 4   2.772702       NA
> > 5   3.401987       NA
> > 6   2.323481 1.965391
> > 7   2.135438 1.724173
> > 8   2.404821       NA
> > 9   3.273616       NA
> > 10  2.058073       NA
> > 11  2.711462       NA
> > 12  1.871119 1.599446
> > 13  2.599939       NA
> > 14  2.174406 1.831875
> > 15  2.377755       NA
> > 16  1.735191 1.713195
> > 17  1.968154 1.836061
> > 18  3.455566       NA
> > 19  2.443566       NA
> > 20  1.865691 1.511079
> > 21  2.219164 1.872533
> > 22  3.011284       NA
> > 23  1.573508 1.312573
> > 24  2.645591       NA
> > 25  1.781628 1.563686
> > 26  4.223059       NA
> > 27  2.617404       NA
> > 28  2.910262       NA
> > 29  2.621719       NA
> > 30  2.084573       NA
> > 31  2.779583       NA
> > 32  2.925445       NA
> > 33  2.768464       NA
> > 
> > 
> > i want to the following operations
> > 
> > Between row 1 and 2
> > ===============
> > abs[(1,1)-(2,1)]+abs[(1,2)-(2,2)]
> > 
> > Between row 1 and 3
> > ===============
> > abs[(1,1)-(3,1)]+abs[(1,2)-(3,1)]
> > 
> > You saw that when there is an NA in the row then both the values are to subtracted from the same number.
> > 
> > Now if there is NA in both the rows then
> > 
> > Between row 3 and 4
> > ===============
> > abs[(3,1)-(4,1)]
> > 
> > > dput(head(res,33))
> > structure(c(0, 0.55141840341246, NA, NA, NA, 0.21052966128458, 
> > 0.12904397545065, NA, NA, NA, NA, 0.26000212187111, NA, 0.19497039591301, 
> > NA, 0.28218135777027, 0.0154679811559602, NA, NA, 0.35379744900785, 
> > 0.19907070014522, NA, 0.84448533806427, NA, 0.38525327569244, 
> > NA, NA, NA, NA, NA, NA, NA, NA, 0.0141081484682, 0, NA, NA, NA, 
> > 0.19642151281638, 0.24959601164407, NA, NA, NA, NA, 0.11000417752067, 
> > NA, 0.18086224744481, NA, 0.13967186326883, 0.0295761296241601, 
> > NA, NA, 0.19294355018347, 0.18496255167702, NA, 0.29306693465181, 
> > NA, 0.0562739047873999, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, 0.55830502624402, 1.10972342965648, NA, NA, NA, 0, 0.42926105079337, 
> > NA, NA, NA, NA, 0.81830714811513, NA, 0.28259146007687, NA, 0.84048638401429, 
> > 0.48465725879666, NA, NA, 0.91210247525187, 0.1971749861944, 
> > NA, 1.40279036430829, NA, 0.94355830193646, NA, NA, NA, NA, NA, 
> > NA, NA, NA, 0.26370416011227, 0.68046237886311, NA, NA, NA, 0.0531744988276899, 
> > 0, NA, NA, NA, NA, 0.38904609732176, NA, 0.0687337641992598, 
> > NA, 0.41122533322092, 0.27917214126823, NA, NA, 0.4828414244585, 
> > 0.0646334599670497, NA, 0.97352931351492, NA, 0.51429725114309, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 0.12411232598887, 
> > 0.29141628154135, NA, NA, NA, 0.0864173352957098, 0.1395918341234, 
> > NA, NA, NA, NA, 0, NA, 0.0708580699241399, NA, 0.2496760407895, 
> > 0.13958030714483, NA, NA, 0.09379532713674, 0.07495837415635, 
> > NA, 0.58448321619316, NA, 0.12525115382133, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, 0.27571356616715, 0.82713196957961, NA, NA, NA, 
> > 0.0155592653715699, 0.1466695907165, NA, NA, NA, NA, 0.53571568803826, 
> > NA, 0, NA, 0.55789492393742, 0.21043837706897, NA, NA, 0.629511015175, 
> > 0.00410030423221008, NA, 1.12019890423142, NA, 0.66096684185959, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, 0.12556371480063, 0.26923704564219, 
> > NA, NA, NA, 0.33609337608521, 0.3892678749129, NA, NA, NA, NA, 
> > 0.0221792358991602, NA, 0.32053411071364, NA, 0, 0.11009573364467, 
> > NA, NA, 0.0716160912375798, 0.32463441494585, NA, 0.562303980294, 
> > NA, 0.10307191792217, NA, NA, NA, NA, NA, NA, NA, NA, 0.0736477674473601, 
> > 0.62506617085982, NA, NA, NA, 0.22599764244054, 0.0553962080032897, 
> > NA, NA, NA, NA, 0.33364988931847, NA, 0.20206579871979, NA, 0.35582912521763, 
> > 0, NA, NA, 0.42744521645521, 0.21453868130118, NA, 0.91813310551163, 
> > NA, 0.4589010431398, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > 0.20705169865167, 0.19762095440461, NA, NA, NA, 0.00347796263291, 
> > 0.0566524614605999, NA, NA, NA, NA, 0.0829393726627998, NA, 0.0120813027386599, 
> > NA, 0.3326154134523, 0.22251967980763, NA, NA, 0, 0.00798099850644984, 
> > NA, 0.49068788905642, NA, 0.13666964539607, NA, NA, NA, NA, NA, 
> > NA, NA, NA, 0.36113004004962, 0.91254844346208, NA, NA, NA, 0.0114589611393598, 
> > 0.23208606459897, NA, NA, NA, NA, 0.62113216192073, NA, 0.0854164738824699, 
> > NA, 0.64331139781989, 0.28748227260226, NA, NA, 0.71492748905747, 
> > 0, NA, 1.20561537811389, NA, 0.74638331574206, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, 0.11337473035895, 0.0992665818907501, NA, 
> > NA, NA, 0.0971549309256299, 0.15032942975332, NA, NA, NA, NA, 
> > 0.01073759562992, NA, 0.0815956655540599, NA, 0.23893844515958, 
> > 0.12884271151491, NA, NA, 0.0936769682927199, 0.08569596978627, 
> > NA, 0, NA, 0.0429926771033502, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, 0.0703820532555999, 0.16616512772002, NA, NA, NA, 0.14014760802898, 
> > 0.19332210685667, NA, NA, NA, NA, 0.0537302727332702, NA, 0.12458834265741, 
> > NA, 0.19594576805623, 0.08585003441156, NA, NA, 0.0314558266845901, 
> > 0.12868864688962, NA, 0.45923206237183, NA, 0, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 0.2646099845617, 
> > 0.81602838797416, NA, NA, NA, 0.0909196620326802, 0.13556600911105, 
> > NA, NA, NA, NA, 0.52461210643281, NA, 0.0111035816054501, NA, 
> > 0.54679134233197, 0.19096221711434, NA, NA, 0.61840743356955, 
> > 0.0794607008933204, NA, 1.10909532262597, NA, 0.64986326025414, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, 0.14935730289419, 0.70077570630665, NA, NA, NA, 0.19938780049321, 
> > 0.0203133274435403, NA, NA, NA, NA, 0.4093594247653, NA, 0.12635626327296, 
> > NA, 0.43153866066446, 0.07570953544683, NA, NA, 0.50315475190204, 
> > 0.18792883935385, NA, 0.99384264095846, NA, 0.53461057858663, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, 0.04990651671688, 0.0500236721055403, 
> > NA, NA, NA, 0.1606231445677, 0.21379764339539, NA, NA, NA, NA, 
> > 0.0742058092719902, NA, 0.14506387919613, NA, 0.17547023151751, 
> > 0.0653744978728401, NA, NA, 0.14759728229907, 0.14916418342834, 
> > NA, 0.34309060675735, NA, 0.0204755365387199, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 0.0451834562530797, 
> > 0.0310753077848798, NA, NA, NA, 0.1653462050315, 0.21852070385919, 
> > NA, NA, NA, NA, 0.0789288697357904, NA, 0.14978693965993, NA, 
> > 0.17074717105371, 0.0606514374090399, NA, NA, 0.16186824239859, 
> > 0.15388724389214, NA, 0.0681912741058703, NA, 0.0251985970025201, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 0.45350165923706, 
> > 1.00492006264952, NA, NA, NA, 0.0759252471688001, 0.32445768378641, 
> > NA, NA, NA, NA, 0.71350378110817, NA, 0.17778809306991, NA, 0.73568301700733, 
> > 0.3798538917897, NA, NA, 0.80729910824491, 0.0923716191874402, 
> > NA, 1.29798699730133, NA, 0.8387549349295, NA, NA, NA, NA, NA, 
> > NA, NA, NA, 0.20869063663638, 0.30476352983468, NA, NA, NA, 0.00183902464820007, 
> > 0.0550135234758899, NA, NA, NA, NA, 0.0845783106475098, NA, 0.0137202407233699, 
> > NA, 0.33425435143701, 0.22415861779234, NA, NA, 0.10714257543007, 
> > 0.00961993649115978, NA, 0.59783046448649, NA, 0.13859840211466, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 0.15548483984177, 0.14137669137357, 
> > NA, NA, NA, 0.0550448214428099, 0.1082193202705, NA, NA, NA, 
> > NA, 0.0313725138529, NA, 0.0394855560712399, NA, 0.2810485546424, 
> > 0.17095282099773, NA, NA, 0.0515668588098999, 0.04358586030345, 
> > NA, 0.04211010948282, NA, 0.0851027865861702, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, 1.77257861901135, 2.32399702242381, NA, NA, NA, 1.21427359276733, 
> > 1.6435346435607, NA, NA, NA, NA, 2.03258074088246, NA, 1.4968650528442, 
> > NA, 2.05475997678162, 1.69893085156399, NA, NA, 2.1263760680192, 
> > 1.41144857896173, NA, 2.61706395707562, NA, 2.15783189470379, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 0.5768660317613, 
> > 1.12828443517376, NA, NA, NA, 0.01856100551728, 0.44782205631065, 
> > NA, NA, NA, NA, 0.83686815363241, NA, 0.30115246559415, NA, 0.85904738953157, 
> > 0.50321826431394, NA, NA, 0.93066348076915, 0.21573599171168, 
> > NA, 1.42135136982557, NA, 0.96211930745374, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 0.19394432473824, 0.17983617627004, 
> > NA, NA, NA, 0.01658533654634, 0.0697598353740299, NA, NA, NA, 
> > NA, 0.0698319987493699, NA, 0.00102607117477005, NA, 0.31950803953887, 
> > 0.2094123058942, NA, NA, 0.01310737391343, 0.00512637540698013, 
> > NA, 0.45793766238951, NA, 0.12356227148264, NA, NA, NA, NA, NA, 
> > NA, NA, NA, 0.09762102156823, 0.21834144893267, NA, NA, NA, 0.30815068285281, 
> > 0.3613251816805, NA, NA, NA, NA, 0.0730748326086801, NA, 0.29259141748124, 
> > NA, 0.0279426932324001, 0.0821530404122699, NA, NA, 0.0207204945280599, 
> > 0.29669172171345, NA, 0.51140838358448, NA, 0.05217632121265, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, 0.17881748976533, 0.73023589317779, NA, 
> > NA, NA, 0.15549234768719, 0.0497735143146802, NA, NA, NA, NA, 
> > 0.43881961163644, NA, 0.0968960764018201, NA, 0.4609988475356, 
> > 0.10516972231797, NA, NA, 0.53261493877318, 0.14403338654783, 
> > NA, 1.0233028278296, NA, 0.56407076545777, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, 0.53092836559544, 1.0823467690079, NA, NA, NA, 
> > 0.0273766606485797, 0.40188439014479, NA, NA, NA, NA, 0.79093048746655, 
> > NA, 0.25521479942829, NA, 0.81310972336571, 0.45728059814808, 
> > NA, NA, 0.88472581460329, 0.16979832554582, NA, 1.37541370365971, 
> > NA, 0.91618164128788, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> > NA, NA, NA, NA, NA, NA, NA), .Dim = c(33L, 124L))
> > 
> > 
> > Thankyou very much iindeed in advance
> > 
> > Elisa



More information about the R-help mailing list