不带intecept的R中的Boxcoxfit函数

我对boxcoxfit函数有问题

我模拟了一些数据，现在我想为回归参数和box-cox变换中的参数创建估计器

我用的是一个包裹

矩阵X有两列，Y没有负值（通过box-cox的逆变换得到）

我使用了一个

boxcoxfit（Y~X）

，答案有4个参数（一个额外的参数用于截取）。当我把截距加到mz矩阵X上，然后再次运行boxcosfit，对于lambda=2，我得到了截距的非有效估计量

这是我的全部代码：

library(geoR)
#optional
set.seed(80974140)
XX=matrix(rnorm(2000,100,12),ncol=2,nrow=1000)
epsilon=rnorm(1000,0,1)
beta=c(0,2,3)
a=2

#invers transformation  
inverz=function(y,a){
  if (a==0) inverz<-exp(y)
  else inverz<-(y*a+1)^(1/a)
  return(inverz)
}

jedna=rep(1,1000)
X=cbind(jedna,XX) #intercept
TY=X%*%beta+epsilon #regression model
head(cbind(TY,X))
Y=inverz(TY,a) #Observed data
summary(Y)
head(cbind(Y,X,epsilon))
boxcoxfit(object=Y,xmat=X)

• 我可以从boxcoxfit中减去截距吗

• 我能得到系数的估计标准偏差吗

谢谢你的回答

PS：对不起，我的英语不好

Fitted parameters:
lambda     beta0     beta1     beta2   sigmasq 
1.9903028 2.0598958 1.9415787 2.8965162 0.9945854 `