[R] R cycle to calculate distances
P.Branco
pjlbranco at yahoo.com
Tue Nov 18 12:37:53 CET 2008
It Works now.
Thanks to all.
P.Branco
Gavin Simpson wrote:
>
> On Mon, 2008-11-17 at 07:20 -0800, P.Branco wrote:
>> Sorry, it does not work.
>>
>> If I do a rnorm I lose the original values of my vectors, and the
>> equation
>> result must be attained by the use of the original values.
>
> Dimitris was generating some dummy data to test that the function
> worked. He had to do this because you didn't provide the data to make up
> r1, r2, etc.
>
> You don't need to do the calls to rnorm, just use your original data,
> stored in those variables/objects.
>
> If r1, r2 already contain your data, you just need to do:
>
> dd <- sqrt(0.723523 * outer(r1, s1, "-")^2 + 0.215518 * outer(r2, s2,
> "-")^2)
> dd
>
> Does that work?
>
> G
>
>>
>> Thanks,
>> P.Branco
>>
>>
>> Dimitris Rizopoulos-4 wrote:
>> >
>> > try this (presented only for two pairs):
>> >
>> > r1 <- rnorm(49)
>> > r2 <- rnorm(49)
>> > s1 <- rnorm(300)
>> > s2 <- rnorm(300)
>> >
>> > dd <- sqrt(0.723523 * outer(r1, s1, "-")^2 + 0.215518 * outer(r2, s2,
>> > "-")^2)
>> > dd
>> >
>> > I hope it helps.
>> >
>> > Best,
>> > Dimitris
>> >
>> > P.Branco wrote:
>> >> Dear all,
>> >>
>> >> I am for the first time trying to work with R, but I have bumped into
>> a
>> >> problem.
>> >>
>> >> I have four vectors:
>> >> r1 (49 values)
>> >> r21 (49 values)
>> >> r22 (49 values)
>> >> r3 (49 values)
>> >> s1 (300 values)
>> >> s21 (300 values)
>> >> s22 (300 values)
>> >> s3 (300 values)
>> >>
>> >> And I would like to calculate the distances from all references (r)
>> sites
>> >> to
>> >> sample (s) site, for that I wrote this equation:
>> >>
>> >>
>> drs=sqrt(((0.723523*(r1[i]-s1[ii]))^2)+(((-0.215518*(r21[i]-s21[ii]))^2)+((0.518612*(r22[i]-s22[ii]))^2))+((0.283655*(r3[i]-s3[ii]))^2))
>> >>
>> >> Now how can I run a cycle in order to calculate the distances?
>> >>
>> >> Could you please help me?
>> >>
>> >> P. Branco
>> >>
>> >
>> > --
>> > Dimitris Rizopoulos
>> > Assistant Professor
>> > Department of Biostatistics
>> > Erasmus Medical Center
>> >
>> > Address: PO Box 2040, 3000 CA Rotterdam, the Netherlands
>> > Tel: +31/(0)10/7043478
>> > Fax: +31/(0)10/7043014
>> >
>> > ______________________________________________
>> > R-help at 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.
>> >
>> >
>>
> --
> %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%
> Dr. Gavin Simpson [t] +44 (0)20 7679 0522
> ECRC, UCL Geography, [f] +44 (0)20 7679 0565
> Pearson Building, [e] gavin.simpsonATNOSPAMucl.ac.uk
> Gower Street, London [w] http://www.ucl.ac.uk/~ucfagls/
> UK. WC1E 6BT. [w] http://www.freshwaters.org.uk
> %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%
>
> ______________________________________________
> R-help at 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.
>
>
--
View this message in context: http://www.nabble.com/R-cycle-to-calculate-distances-tp20540340p20557906.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list