[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
struct covariance
From: |
Jason Stover |
Subject: |
struct covariance |
Date: |
Wed, 1 Jun 2011 17:53:24 -0400 |
User-agent: |
Mutt/1.5.18 (2008-05-17) |
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
- struct covariance,
Jason Stover <=