C言語を使って最小二乗法で回帰直線を求めたい

Question

C言語を使って最小二乗法で回帰直線を求めたいのですが、分かりません。

今ここまで書けております。

#include <stdio.h>
#include <math.h>

double total(double *x, int n)
{
    int j;
    double ans=0.0;
    for (j=0;j<n;j++){
        ans+=x[j];
    }
    return ans;
}
double lsm(double *x, double *y, int n, double *a, double *b)
{
    //Least square method
    int j;
    double xy[n],x2[n],y2[n];
    double xsum,ysum,xysum,x2sum,y2sum;
    double denom,cc;
    for (j=0;j<n;j++){
        xy[j]=x[j]*y[j];
        x2[j]=x[j]*x[j];
        y2[j]=y[j]*y[j];
    }
    xsum=total(x,n);
    ysum=total(y,n);
    xysum=total(xy,n);
    x2sum=total(x2,n);
    y2sum=total(y2,n);
    denom=n*x2sum-xsum*xsum;
    (*a)=(n*xysum-xsum*ysum)/denom;
    (*b)=(x2sum*ysum-xysum*xsum)/denom;
    //Correlation coefficient
    denom=sqrt((n*x2sum-xsum*xsum)*(n*y2sum-ysum*ysum));
    cc=(n*xysum-xsum*ysum)/denom;
    return cc;
}

int main(void){
    // Your code here!
    // lwsson 11-1
    double x[10]={1,2,3,4,5,6,7,8,9,10};
    double y[10]={10.083,5.760,34.672,11.894,42.462,
        37.928,17.584,54.281,61.774,74.931};
        
        double a,b,c; // Regression line y=ax+b. c=cor. coef.
        c=lsm(x,y,10,&a,&b);
        printf("y = (%f)x + (%f)\n",a,b);
        printf("Correlation coefficient = %f\n" ,c);   
}

上記までのコードを変更して以下を求めたいです。

答えは


となります。

どこをどう変更して入力したらいいのか分からないので、教えていただけると助かります。答えまでよろしくお願いします。

Answer 1

これを入力すればいいと思います。

Comments

[fana]

> C言語を使って最小二乗法で回帰直線を求めたい

という話なので，
最小二乗法で回帰直線を求める部分を別で実施した結果を用いるのでは
質問に対する回答とはならないように思えます．

[TaroToyotomi]

fanaさんの回答と同じ意味で、yの列の中身を配列に設定すればという意味の回答でした。
グラフはプログラムを走らせた結果を確認する意味でダブルチェックになればいいなとの思いで適当に付けました。

Answer 2

main関数内の配列yの内容を，その指示通りの値にしてやればよいだけなのではないでしょうか．

Comments

[yuharu]

どのように記述したら良いのか教えて頂けないでしょうか？よろしくお願いします。

[episteme]

y[i] = (int)(53 * y[i-1] + 29) % 199;
コレ↑を i=1～9 でくりかえす。

[yuharu]

繰り返すところからの記述を教えていただけないでしょうか？お手数おかけします。

[episteme]

for ループがわからんのですか? そんなはずないよね、関数 total, lsm を書いたんだから。

[yuharu]

これはもともと指定されていたものになり、そこからの出題になります。

[episteme]

いや、自分で書いたものじゃないとはいえ、 関数 total のコード読めば for ループの回し方くらい理解できるでしょ。わからんなら教本読めばいい。

# てかアナタ以前に漸化式の質問し、僕の回答を理解したんじゃなかったけか? そこでも forループ使ってる。
https://teratail.com/questions/303347

[fana]

提示コードみたいに yの値を全部を直値で書く（手計算で求めた結果をハードコーディングする）のでも良いのでは？（ダメなの？）