Question

现在我这样做:

d = approx(density(csvdata[, X ],n=5000),xout=csvdata[, X ])
dfact = 40/max(d$y)
for(i in 1:nrow(csvdata)) {
  d$y[i] = (d$y[i]*dfact)-20
}

我在这里做的是转轨密度功能,其频率总是高于0,从我20台的海图底下显示,它始终适合20周年,因此,我更容易发现该线中的任何不合规定之处。既然你们能够看到我这样做,但也许会有一个线上的人会这样做?

Answer 1

<代码>y是一种病媒,* 和-为矢量功能,因此,您不必随时间推展计算的次数。他们一劳永逸:

d$y <- (d$y * dfact) - 20

或更好(无<条码>d$)

d <- within(d, y <- (y * dfact) - 20)

<代码>dfact是一个轴心(R中一个长的1个矢量,因为它有单独的体质概念),但R将把<条码>dfact (<>release/em> 是R中正确的术语)扩展至必要的长度,使上述计算方法能够正常进行。