質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

90.47%

  • C

    4659questions

    C言語は、1972年にAT&Tベル研究所の、デニス・リッチーが主体となって作成したプログラミング言語です。 B言語の後継言語として開発されたことからC言語と命名。そのため、表記法などはB言語やALGOLに近いとされています。 Cの拡張版であるC++言語とともに、現在世界中でもっとも普及されているプログラミング言語です。

プログラムのエラーの原因がわかりません

解決済

回答 1

投稿

  • 評価
  • クリップ 0
  • VIEW 329

jackie687456

score 10

前提・実現したいこと

データ(X_i,Y_i) (i=1,...,n)が与えられたとき、単回帰Y_i=β、X_i+αの回帰係数αとβを求めるプログラムを作成し、回帰係数βについてt検定を行うためのp値を求めよ。

という課題でプログラムを組んだのですが、以下のようなエラーが出てしまい実行できません。
簡単な何かを見落としているような気もするのですが、どうしても見つけられず、お力をお借りしたいです。
よろしくお願いいたします。

発生している問題・エラーメッセージ

98行目」で記述エラーを発見しました。
「function pointer」を付け忘れています。

該当のソースコード

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

#define N_MAX 100

void Random(double I, int n, double *R){
    int a, k;
    double m;

    a = 69621;
    m = pow(2, 32);

    for (k = 0; k<n; k++){
        I = a*I - floor(a*I / m)*m;
        R[k] = I / (m - 1);
    }
}

double f(double t, double a, double b){
    double y;
    y = pow(t, a - 1)*pow(1 - t, b - 1);
    return y;
}

double Simpson(double R, double a, double b){
    int n;
    double h, I, sum1 = 0, sum2 = 0;

    h = (R - 0) / 100000.0;

    for (n = 1; n <= 100000 / 2 - 1; n++){
        sum1 += f(0 + 2 * n*h, a, b);
    }

    for (n = 1; n <= 100000 / 2; n++){
        sum2 += f(0 + (2 * n - 1)*h, a, b);
    }

    I = (h / 3)*(f(0, a, b) + f(R, a, b) + 2 * sum1 + 4 * sum2);

    return I;
}

double beta(double x, double a, double b){
    double be;
    be = Simpson(x, a, b) / Simpson(1.0, a, b);

    return be;
}

int main(void){
    int i, n = N_MAX;
    double alpha, beta, numer, denom, total, v, averx, avery, disp, t, w, a, b, p;
    double x[N_MAX], y[N_MAX];

    Random(1.0, n, x);
    Random(100.0, n, y);

    total = 0;

    for (i = 0; i < n; i++){
        total += x[i];
        averx = total / n;
    }

    total = 0;

    for (i = 0; i < n; i++){
        total += y[i];
        avery = total / n;
    }

    numer = 0;
    denom = 0;

    for (i = 0; i < n; i++){
        numer += (y[i] - avery)*(x[i] - averx);
        denom += (x[i] - averx)*(x[i] - averx);
    }

    beta = numer / denom;
    alpha = avery - beta*averx;

    v = 0;

    for (i = 0; i<n; i++){
        v += (y[i] - alpha - beta*x[i])*(y[i] - alpha - beta*x[i]);
        disp = v / n;
    }

    t = beta / (sqrt(disp) / sqrt(denom));

    w = (t + sqrt(t*t + n - 2)) / (2 * sqrt(t*t + n - 2));
    a = (n - 2) / 2;
    b = (n - 2) / 2;

    p = 1 - beta(w, a, b);  //エラー該当箇所

    printf("t = %f, p = %f\n", t, p);

    printf("\n");


    return 0;
}
  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

回答 1

checkベストアンサー

+2

beta が変数としても関数としても定義されているからです。 この場合、 beta はローカル変数の beta と解釈され、もちろんそれは double 型であって関数型ではないので関数として呼び出すことができません。 名前を変えるだけで回避できます。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2017/06/04 22:55

    ご指摘の通りでした…。稚拙な質問で申し訳ありませんでした。ご回答ありがとうございました。

    キャンセル

15分調べてもわからないことは、teratailで質問しよう!

  • ただいまの回答率 90.47%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

同じタグがついた質問を見る

  • C

    4659questions

    C言語は、1972年にAT&Tベル研究所の、デニス・リッチーが主体となって作成したプログラミング言語です。 B言語の後継言語として開発されたことからC言語と命名。そのため、表記法などはB言語やALGOLに近いとされています。 Cの拡張版であるC++言語とともに、現在世界中でもっとも普及されているプログラミング言語です。