小弟是R的新手,请高手看看下面的程序有什么问题?
#box-cox变换
m=function(x,k)
{
n=length(x);y=rep(0,n)
for(i in 1:n)
{
if (k!=0) y[i]=(x[i]^k-1)/k else y[i]=log(x[i],10)
}
y
}
x=rnorm(100)
m(x,1) #m(x,1)等于前面x-1,但是y却不是x-1,搞不明白,请教高手指点一下。
y
m111=data.frame(x,y)
plot(y~x,data=m111,xlim=range(x),ylim=range(y),pch=16,type="p")##请问为什么不能用得到光滑的线条图像?
谢谢!