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_bo...@hotmail.com> To: "smartpink...@yahoo.com" <smartpink...@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: smartpink...@yahoo.com > Subject: Re: a similar question > To: eliza_bo...@hotmail.com > CC: r-help@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_bo...@hotmail.com> > To: "smartpink...@yahoo.com" <smartpink...@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: smartpink...@yahoo.com > > Subject: Re: a similar question > > To: eliza_bo...@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_bo...@hotmail.com> > > To: "smartpink...@yahoo.com" <smartpink...@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 ______________________________________________ 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.