[R] how to overlook the zero in the denominator

Christoph Buser buser at stat.math.ethz.ch
Wed Jul 27 11:09:53 CEST 2005


What about

x<-c(0,0,0.28,0.55,1.2,2,1.95,1.85,
     1.6,0.86,0.78,0.6,0.21,0.18)
y<-c(0,0,0,0.53,1.34,1.79,2.07,1.88,
    1.52,0.92,0.71,0.55,0.32,0.19)

sum(((x-y)^2/x)[x!=0])

Regards,

Christoph Buser

--------------------------------------------------------------
Christoph Buser <buser at stat.math.ethz.ch>
Seminar fuer Statistik, LEO C13
ETH (Federal Inst. Technology)	8092 Zurich	 SWITZERLAND
phone: x-41-44-632-4673		fax: 632-1228
http://stat.ethz.ch/~buser/
--------------------------------------------------------------


Chun-Ying Lee writes:
 > Dear R users:
 > 
 > I have two set of data, as follow:
 > x<-c(0,0,0.28,0.55,1.2,2,1.95,1.85,
 >      1.6,0.86,0.78,0.6,0.21,0.18)
 > y<-c(0,0,0,0.53,1.34,1.79,2.07,1.88,
 >     1.52,0.92,0.71,0.55,0.32,0.19)
 > i<-1:length(x)
 > 
 > I want to sum each (x[i]-y[i])^2/x[i] together, 
 > like:
 > >Sum <-sum((x[i]-y[i])^2/x[i])
 > >Sum
 > [1] NaN
 > 
 > Because the denominator shoud not be zero.
 > So I want to overlook those when x[i]=0,
 > and just to sum those x[i] not equal to 0.
 > What should I do?
 > Any suggestion.
 > Thanks in advance !!
 > 
 > ______________________________________________
 > R-help at stat.math.ethz.ch mailing list
 > https://stat.ethz.ch/mailman/listinfo/r-help
 > PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html




More information about the R-help mailing list