在 Julia 中运行 Wald 测试的问题

Question

我在运行以下函数以对我的数据执行 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, T, 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,n);