示例#1
0
gsl_vector *cooks_distance(apop_model *in){
    apop_data  *c = apop_data_copy(in->data);
    apop_ols->prep(in->data, in);
    ols_data = in->data;
    predicted = project(in->data, in);
    p_dot_mse  = c->matrix->size2 * sum_squared_diff(in->data->vector, predicted);
    return jack_iteration(c->matrix, cook_math);
}
示例#2
0
static gsl_vector *cooks_distance(apop_model *in){
  apop_data  *c = apop_data_copy(in->data);
  apop_ols->prep(in->data, in);
  ols_data = in->data;
  apop_data * t = project(in->data, in);
  predicted = t->vector;
  p_dot_mse  = c->matrix->size2 * sum_squared_diff(in->data->vector, t->vector);
  
  gsl_vector* out = jack_iteration(c->matrix, cook_math);
  apop_data_free(t);
  return out;
}