I'm working on the type 3 sums of squares and have a question
about the covariance struct.
First: my data look like this:
data list list / v0 to v2.
begin data
3.2 1 1
3.1 1 1
3.3 1 2
3.4 1 2
3.2 1 3
3.3 1 3
3.3 1 4
3.2 1 4
2.8 2 1
2.9 2 1
3.3 2 2
3.0 2 2
3.1 2 3
3.2 2 3
3.2 2 4
3.1 2 4
end data
GLM v0 by v1 v2
/INTERCEPT = include.
Currently in run_glm we have:
struct covariance *cov = covariance_2pass_create (cmd->n_dep_vars,
cmd->dep_vars,...);
Then later we have
gsl_matrix *cm = covariance_calculate_unnormalized (cov);
Later, when I need to find dimensions for another matrix, cov tells me
there is only 1 variable represented in the struct, but cm->size1 is
5. Shouldn't it be 1, based on the above call to
covariance_2pass_create?
Anyway, I created a list of all the variables (dependent and explanatory) and
did this:
struct covariance *cov = covariance_2pass_create (n_vars, vars,...);
vars contains v0, v1 and v2. And now, cm->size1 is 7. But shouldn't it be
3?
Also,
If I'm missing something, please let me know. I need to find a ay to
correctly count the number of dimensions in the covariance matrix,
either from the struct or cm.
-Jason
_______________________________________________
pspp-dev mailing list
[email protected]
https://lists.gnu.org/mailman/listinfo/pspp-dev
