R的诡异错误

A pitfall in R

先看下面的代码：

freq <- 0.29
altAllele <- freq * 200 * 2
print(altAllele) ## this will output 116
print(length(rep(1, altAllele))) ## what is the output of this??

最后一行输出的是115， 不是116。为什么？
答案是：

altAllele是一个numeric类型的数，尽管打印出来是116，但在电脑里存储的是一个比整数116稍小的浮点数。因此rep(1, altAllele)是一个长度是115的vector。

R的程序有很多陷阱。这个陷阱是因为R不需要定义数据类型，但它的底层还是有类型的。因此写R程序一定要注意这点。解决方法如下：

freq <- 0.29
altAllele <- as.integer(freq * 200 * 2)
print(altAllele) ## this will output 116
print(length(rep(1, altAllele))) ## what is the output of this??

Python里有没有类似的陷阱？结论是一样有 (Python 2.7.5 NumPy 1.8.0)：

freq = .29
altAllele = freq * 200 * 2
len(np.repeat(1, altAllele)) ## output 115 here.

但是Python里打印altAllele的值不是116，而是115.99999999999999。这总比R打印出16要好一点。