Функция Boxcoxfit в R без intecept

Question

У меня проблема с функцией boxcoxfit.

Я смоделировал некоторые данные, и теперь я хочу создать оценки для параметров регрессии и параметра в преобразовании Бокса-Кокса.

Я использую пакет GeoR

У меня есть матрица X с 2 столбцами и Y с неотрицательными значениями (которые я получаю путем обратного преобразования box-cox).

Я использую boxcoxfit(Y~X) и ответ имеет 4 параметра (один дополнительный для перехвата). Когда я добавляю перехват к матрице mz и снова запускаю boxcosfit, для лямбда =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)

И вывод:

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

• Могу ли я как-то вычесть пересечение из boxcoxfit?

• Могу ли я получить оценку стандартного отклонение для коэффициентов?

Спасибо за ваши ответы

PS: извините за мой плохой английский