[R] rounding problem

Peterko lanikpeter at gmail.com
Fri Feb 27 11:55:34 CET 2009


all.equal is what i need, many thanks to help me

baptiste auguie-2 wrote:
> 
> Hi,
> 
> you probably want to use ?all.equal instead of "=="
> 
> I couldn't run your example, though
> 
> Hope this helps,
> 
> baptiste
> 
> On 27 Feb 2009, at 10:32, Peterko wrote:
> 
>>
>> hi i am creating some variables from same data, but somewhere is  
>> different
>> rouding.
>> look:
>> P = abs(fft(d.zlato)/480)^2
>> hladane= sort(P,decreasing=T)[1:10]/480
>>
>> pozicia=c(0,0,0,0,0)
>> for (j in 1:5){ for (i in 2:239){
>>  if (P[i]/480==hladane[2*j-1]){pozicia[j]=i-1}}}
>> period=479/pozicia
>>
>>> P[2]/334
>> [1] 0.0001279107
>>> hladane[1]
>> [1] 0.0001279107
>>> P[2]/334==hladane[1]
>> [1] FALSE
>>> abs(P[2]/334 - hladane[1]) < 0.0000001
>> [1] TRUE
>>
>> It is possible to avoid it ?
>> I know in this exam i can use 2x if to eliminate this rouding, but i  
>> need to
>> fix it in general.
>> --
>> View this message in context:
>> http://www.nabble.com/rounding-problem-tp22243179p22243179.html
>> Sent from the R help mailing list archive at Nabble.com.
>>
>> ______________________________________________
>> 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.
> 
> _____________________________
> 
> Baptiste Auguié
> 
> School of Physics
> University of Exeter
> Stocker Road,
> Exeter, Devon,
> EX4 4QL, UK
> 
> Phone: +44 1392 264187
> 
> http://newton.ex.ac.uk/research/emag
> 
> ______________________________________________
> 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/rounding-problem-tp22243179p22243567.html
Sent from the R help mailing list archive at Nabble.com.




More information about the R-help mailing list