Generates \(\beta\) estimates for MLE using a conditioning approach with patterning support.
bsm.fit(x, y, z, pattern)
An \(N \times (P + F)\) design matrix, where \(F\) is the number of columns conditioned on. This is equivalent to the multiplication of \(xyzb\).
The \(N \times (Q - F)\) matrix of observations, where \(F\) is the number of columns conditioned on. This is equivalent to the multiplication of \(Yz_a\).
A \((Q - F) \times L\) design matrix, where \(F\) is the number of columns conditioned on.
An optional \(N-F x F\) matrix of 0's and 1's indicating which elements of \(\beta\) are allowed to be nonzero.
A list with the following components:
The least-squares estimate of \(\beta\).
The \((P+F)\times L\) matrix with the \(ij\)th element being the standard error of \(\hat{\beta}_ij\).
The \((P+F)\times L\) matrix with the \(ij\)th element being the t-statistic based on \(\hat{\beta}_ij\).
The estimated covariance matrix of the \(\hat{\beta}_ij\)'s.
A \(p\)-dimensional vector of the degrees of freedom for the \(t\)-statistics, where the \(j\)th component contains the degrees of freedom for the \(j\)th column of \(\hat{\beta}\).
The \((Q - F) \times (Q - F)\) matrix \(\hat{\Sigma}_z\).
The \(Q \times Q\) residual sum of squares and crossproducts matrix.