[GenABEL-dev] [Genabel-commits] r658 - in pkg/VariABEL: R src

[Yurii Aulchenko]

Somewaht lost here:

-----------------------------------------------
> @@ -213,10 +216,15 @@
>  //     for(int i=0 ; i<p ; i++) {std::cout<<"betas["<<i<<"]="<<betas[i]<<", se["<<i<<"]="<<se[i]<<"\n";}
>
>
> -       *chi2 = (betas[p-1]/se[p-1])*(betas[p-1]/se[p-1]);
> +       std::cout<<"2\n";
> +       *chi2 = (betas_disp[1]/se_disp[1])*(betas_disp[1]/se_disp[1]);
> +       std::cout<<"3\n";
>        *df = 1;
> -       betas[p] = beta_snp;
> -  se[p]           = sebeta_snp;
> +       std::cout<<"4\n";
> +       betas[p] = betas_disp[1];
> +       std::cout<<"5\n";
> +  se[p]           = se_disp[1];
> +       std::cout<<"6\n";
>
>        }
-----------------------------------------------

where do you initialize p?

You say @line 88-89:

>	int p;
>	p = *p_;

but I can not quite track where you assign a value to it

is that something happening at line 185 (call to 'linear_regression')?

> 	        /*input variables:*/     trait, design_matrix, p_, &nids_nona,