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

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

新規登録して質問してみよう
ただいま回答率
85.50%
C

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

Q&A

解決済

1回答

1464閲覧

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

jackie687456

総合スコア17

C

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

0グッド

0クリップ

投稿2017/06/04 13:29

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

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

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

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

###該当のソースコード

C言語

1#include<stdio.h> 2#include<math.h> 3#include<stdlib.h> 4 5#define N_MAX 100 6 7void Random(double I, int n, double *R){ 8 int a, k; 9 double m; 10 11 a = 69621; 12 m = pow(2, 32); 13 14 for (k = 0; k<n; k++){ 15 I = a*I - floor(a*I / m)*m; 16 R[k] = I / (m - 1); 17 } 18} 19 20double f(double t, double a, double b){ 21 double y; 22 y = pow(t, a - 1)*pow(1 - t, b - 1); 23 return y; 24} 25 26double Simpson(double R, double a, double b){ 27 int n; 28 double h, I, sum1 = 0, sum2 = 0; 29 30 h = (R - 0) / 100000.0; 31 32 for (n = 1; n <= 100000 / 2 - 1; n++){ 33 sum1 += f(0 + 2 * n*h, a, b); 34 } 35 36 for (n = 1; n <= 100000 / 2; n++){ 37 sum2 += f(0 + (2 * n - 1)*h, a, b); 38 } 39 40 I = (h / 3)*(f(0, a, b) + f(R, a, b) + 2 * sum1 + 4 * sum2); 41 42 return I; 43} 44 45double beta(double x, double a, double b){ 46 double be; 47 be = Simpson(x, a, b) / Simpson(1.0, a, b); 48 49 return be; 50} 51 52int main(void){ 53 int i, n = N_MAX; 54 double alpha, beta, numer, denom, total, v, averx, avery, disp, t, w, a, b, p; 55 double x[N_MAX], y[N_MAX]; 56 57 Random(1.0, n, x); 58 Random(100.0, n, y); 59 60 total = 0; 61 62 for (i = 0; i < n; i++){ 63 total += x[i]; 64 averx = total / n; 65 } 66 67 total = 0; 68 69 for (i = 0; i < n; i++){ 70 total += y[i]; 71 avery = total / n; 72 } 73 74 numer = 0; 75 denom = 0; 76 77 for (i = 0; i < n; i++){ 78 numer += (y[i] - avery)*(x[i] - averx); 79 denom += (x[i] - averx)*(x[i] - averx); 80 } 81 82 beta = numer / denom; 83 alpha = avery - beta*averx; 84 85 v = 0; 86 87 for (i = 0; i<n; i++){ 88 v += (y[i] - alpha - beta*x[i])*(y[i] - alpha - beta*x[i]); 89 disp = v / n; 90 } 91 92 t = beta / (sqrt(disp) / sqrt(denom)); 93 94 w = (t + sqrt(t*t + n - 2)) / (2 * sqrt(t*t + n - 2)); 95 a = (n - 2) / 2; 96 b = (n - 2) / 2; 97 98 p = 1 - beta(w, a, b); //エラー該当箇所 99 100 printf("t = %f, p = %f\n", t, p); 101 102 printf("\n"); 103 104 105 return 0; 106} 107

気になる質問をクリップする

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

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

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

ベストアンサー

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

投稿2017/06/04 13:43

SaitoAtsushi

総合スコア5437

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

jackie687456

2017/06/04 13:55

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

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

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

ただいまの回答率
85.50%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問