one addition ... > Hello Ricardo, > > another solution could be using package distr: > > library(distr) > A <- c(18,18,18,19,20,21,22,23,24,25,26,27,28) > DA <- DiscreteDistribution(A) > # maybe > # support(DA) > # plot (DA) > > B <- c(82,83,84,85,85,86,87,88,89,90,91,91,92) > DB <- DiscreteDistribution(B) > # support(DB) > # plot(DB) > > DC <- DB - DA# convolution with of DB with (-DA) > # support(DC) > # plot(DC) > to get the probabilities for the C-values: p(DC)(support(DC))
> hth, > Matthias > > [EMAIL PROTECTED] wrote: > >>> i have two vectors (A and B) and i need create another vector (C) >>> from the subtraction of A's values with the B's values. How can i >>> estimate the probability of C's values if i have differents values >>> combinations of A and B that can result in the same value? like 90 >>> -20 = 70 and 91 - 21 = 70. >>> >>> example: >>> B = {18 18 18 19 20 21 22 23 24 >>> 25 26 27 28....} >>> A = {82 83 84 85 85 86 87 88 89 >>> 90 91 91 92 ....} >>> >>> >> Try this: >> A = c(18,18,18,19,20,21,22,23,24,25,26,27,28) >> B = c(82,83,84,85,85,86,87,88,89,90,91,91,92) >> lenA = length(A) >> lenB = length(B) >> AA = rep(A, lenB) >> BB = rep(B, each=lenA) >> table(AA-BB)/(lenA*lenB) >> >> Regards, >> Richie. >> >> Mathematical Sciences Unit >> HSL >> >> >> ------------------------------------------------------------------------ >> ATTENTION: >> >> This message contains privileged and confidential inform...{{dropped:20}} >> >> ______________________________________________ >> 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. >> >> > > ______________________________________________ > 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. > ______________________________________________ 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.