在 Julia 中运行 Wald 测试的问题

如何解决在 Julia 中运行 Wald 测试的问题

我在运行以下函数以对我的数据执行 Wald 测试时遇到问题。它总是向我显示未定义变量 theta 的错误消息。我试图在 struct WaldTestFun{F,T,Z} 中定义它，但这不起作用。我得到的错误信息是
错误消息 = undefvarerror：theta 未定义
这是代码的一部分：

 using Linearalgebra
 using Optim
 using prettytables
 using Printf
 using Statistics
 using StatsBase
 using StatsFuns


 struct WaldTest
 tbl::NamedTuple
 rankmin::Int64
 rankₘₐₓ::Int64
 end

 struct WaldTestFun{F,Z}
 f::F
 r::Int64
 vecsigma::T
 Vhat::Z
 end

 (wf::WaldTestFun)(theta) = wf.f(theta,wf.r,wf.vecsigma,wf.Vhat) #Here the error occurs

 function waldobjfun(th,r,vecsigma,Vhat)
 r,k = size(theta)                                                 #Here the error occurs
 theta = reshape(theta,r+1,length(th)÷(r+1))
 sigmamat = diagm(0=>theta[1,:].^2) .+ theta[2:r+1,:]'*theta[2:r+1,:]
 tempsigma = sigmamat[findall(tril(ones(size(sigmamat))).==1)]
 (vecsigma -tempsigma)' /Vhat *(vecsigma - tempsigma)
 end

 X = randn(100,10);
 fm = Factotum.FactorModel(X,3)

 function waldtest(fm::FactorModel,minrank::Int = 0,maxrank::Int = 2)
 X = copy(fm.X)
 T,n = size(X)
 ## normalize factor
 Xs = X / diagm(0=>sqrt.(diag(cov(X))))
 covX = cov(Xs)
 meanX = mean(Xs,dims=1)
 vecsigma = vech(covX)
 bigN = length(vecsigma)
 Vhat = Array{Float64}(undef,bigN,bigN)
 varvecsig = zeros(n,n,n);