thr3ads.net - R help - [R] a similar question [Feb 2013]

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

arun

2013-Feb-28 00:14 UTC

[R] a similar question

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

arun

2013-Mar-26 23:02 UTC

head link

[R] a similar question

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

arun

2013-Mar-27 18:04 UTC

head link

[R] a similar question

Hi Elisa,
Try this:
dat1<- read.csv("DQV.csv",sep="\t")

res<- do.call(cbind,lapply(seq_len(nrow(dat1)),function(i)
do.call(rbind,lapply(split(rbind(dat1[i,],dat1[-i,]),1:nrow(rbind(dat1[i,],dat1[-i,]))),function(x)
{x1<- rbind(dat1[i,],x); colnames(x1)<-
gsub("[.]","",colnames(x1));
if({indx<- colSums(is.na(x1)); indx[2]==0 & indx[3]==1}) # 2 peaks 1 peak
comparison
{x2<- x1[rev(order(x1$Peak2v, x1$Npeak2v)),];
with(x2,{abs(Peak1v[1]-Peak1v[2])+abs(Peak1v[1]-Peak2v[2]) +
abs(Npeak1v[1]-Npeak1v[2])+abs(Npeak1v[1]-Npeak2v[2]) })
}

else {x1[is.na(x1)]<-0; x1; #cases where peaks are similar
with(x1,{abs(Peak1v[1]-Peak1v[2])+ abs(Peak2v[1]-Peak2v[2])+
abs(Npeak1v[1]-Npeak1v[2])+ abs(Npeak2v[1]-Npeak2v[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.0000000 0.8024471 2.3537210 2.1929718 3.7746302
#2 0.8024471 0.0000000 3.1233237 2.9954189 4.5442329
#3 2.3537210 3.1233237 0.0000000 0.4094026 0.8108279
#4 2.1929718 2.9954189 0.4094026 0.0000000 0.8571442
#5 3.7746302 4.5442329 0.8108279 0.8571442 0.0000000

A.K.






________________________________
From: eliza botto <eliza_botto at hotmail.com>
To: "smartpink111 at yahoo.com" <smartpink111 at yahoo.com> 
Sent: Wednesday, March 27, 2013 12:07 PM
Subject: RE: a similar question



Dear Arun,
The last question regarding this distance measurement is that how to modify the
following codes for calculation of distance if we have 2 peaks atmost. The data
file is attached

Thanks in advance
Elisa

> Date: Tue, 26 Mar 2013 16:02:01 -0700
> From: smartpink111 at yahoo.com
> Subject: Re: a similar question
> To: eliza_botto at hotmail.com
> CC: r-help at r-project.org
> 
> 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.
> 
> 
> 
>

R help - Feb 2013 - a similar question

[R] a similar question

[R] a similar question

[R] a similar question