每天一个R语言小技巧（1）——浮点数比较

代码

root<-function(x){
	a=0
	b=1
	while( a != b ){
		a = b
		b=0.5*(b + x/b)
	}
	b
}

问题

上面这段代码用于求x的root值。在这里有一个十分不专业的问题， 那就是a != b的写法。
对于numeric的变量来说，使用 ！= 或者是 = = 是非常不专业的。针对整型数，这会返回正确的结果，但是由于机器精度，对于浮点数来说，他就“可能”会返回一个错误结果，于是在上述代码中，该问题则会导致程序陷入While的无限循环。

解决方式

abs(a-b)<exp(-30)

对于向量形式的两个变量，则可使用 max(abs(a-b))<exp(-30) 来返回两个向量中的所有元素是否完全一致

PS：对于exp(n)中的n，可以按需选用。例如说，如果两个浮点数的差距只有10的负十几次方的话，n可以改为-40，-50，-60…，省事一点的话，可以直接选到-100之类的，省去麻烦。