/
BetaPosterior.cs
135 lines (130 loc) · 5.45 KB
/
BetaPosterior.cs
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
using System.Windows.Forms;
using System.Windows.Forms.DataVisualization.Charting;
namespace VISAP商科应用
{
public partial class BetaPosterior : Form
{
int Alpha = 0;
int Beta = 0;
double MaxValue = 0;
public BetaPosterior()
{
InitializeComponent();
label4.Text = "输入回答\"" + Binomial.Classification[0] + "\"的次数(用逗号分隔):";
Alpha = Binomial.Alpha;
Beta = Binomial.Beta;
chart_PosteriorDensity.Series.Clear();
DrawBeta(Alpha - Binomial.CountTimes[0], Beta - Binomial.CountTimes[1],"先验分布");
DrawBinomial(Binomial.CountTimes[0] + Binomial.CountTimes[1], Binomial.CountTimes[0], "似然函数");
DrawBeta(Alpha, Beta, "后验分布");
}
private void DrawBeta(int Alpha,int Beta,string Name){
Series series = new Series(Name);
series.BorderWidth = 3;
series.ChartType = SeriesChartType.Spline;
double Temp = 0;
for (double i = 0.01; i < 1; i = i + 0.05)
{
Temp = Stat.BetaPDF(i, Alpha, Beta);
if (MaxValue < Temp)
{
MaxValue = Temp;
}
series.Points.AddXY(i,Temp);
}
chart_PosteriorDensity.Series.Add(series);
var XAxis = chart_PosteriorDensity.ChartAreas[0].AxisX;
XAxis.MajorGrid.Enabled = false;
var YAxis = chart_PosteriorDensity.ChartAreas[0].AxisY;
YAxis.MajorGrid.Enabled = false;
XAxis.Maximum = 1;
XAxis.Minimum = 0;
YAxis.Maximum = Math.Ceiling(MaxValue);
}
private void DrawBinomial(int n, int y, string Name)
{
Series series = new Series(Name);
series.BorderWidth = 3;
series.ChartType = SeriesChartType.Spline;
double MaxValue = 0;
double Temp = 0;
for (double i = 0.01; i < 1; i = i + 0.05)
{
Temp = n * Stat.BinomialPDF(i,n,y);
if (MaxValue < Temp)
{
MaxValue = Temp;
}
series.Points.AddXY(i,Temp );
}
chart_PosteriorDensity.Series.Add(series);
var XAxis = chart_PosteriorDensity.ChartAreas[0].AxisX;
XAxis.MajorGrid.Enabled = false;
var YAxis = chart_PosteriorDensity.ChartAreas[0].AxisY;
YAxis.MajorGrid.Enabled = false;
XAxis.Maximum = 1;
XAxis.Minimum = 0;
YAxis.Maximum = Math.Ceiling(MaxValue);
}
private void button_Percentile_Click(object sender, EventArgs e)
{
char[] separator = { ',' };
string[] Percentiles = textBox_Percentiles.Text.Split(separator);
double Temp = 0;
StringBuilder Str = new StringBuilder();
foreach (string EachPer in Percentiles)
{
if (double.TryParse(EachPer, out Temp))
{
if (Temp > 0 && Temp < 1)
{
Str.Append(StrManipulation.PadLeftX(Temp.ToString() + "分位数: ", ' ', 12));
Str.Append(StrManipulation.PadLeftX(MathV.NumberPolish(Stat.BetaUa(1 - Temp, Alpha, Beta).ToString()),' ',12));
Str.Append("\r\n");
}
}
}
MainForm.S.richTextBox1.AppendText(Str.ToString());
MainForm.S.richTextBox1.Select();//让RichTextBox获得焦点
MainForm.S.richTextBox1.Select(MainForm.S.richTextBox1.TextLength, 0);//将插入符号置于文本结束处
MainForm.S.richTextBox1.ScrollToCaret();
}
private void button_predict_Click(object sender, EventArgs e)
{
char[] separator = { ',' };
string[] Prediction = textBox_success.Text.Split(separator);
int Temp = 0;
int n = Convert.ToInt32(textBox_NewSampleSize.Text);
StringBuilder Str = new StringBuilder();
Str.Append(StrManipulation.PadLeftX("次数",' ',12));
Str.Append("\t");
Str.Append(StrManipulation.PadLeftX("预测概率值", ' ', 12));
Str.Append("\r\n");
foreach (string EachPred in Prediction)
{
if (int.TryParse(EachPred, out Temp))
{
if (Temp >= 0 && Temp <= n)
{
Str.Append(StrManipulation.PadLeftX(Temp.ToString(), ' ', 12));
Str.Append("\t");
Str.Append(StrManipulation.PadLeftX(MathV.NumberPolish(Stat.BetaPrediction(Alpha,Beta,n,Temp).ToString()), ' ', 12));
Str.Append("\r\n");
}
}
}
MainForm.S.richTextBox1.AppendText(Str.ToString());
MainForm.S.richTextBox1.Select();//让RichTextBox获得焦点
MainForm.S.richTextBox1.Select(MainForm.S.richTextBox1.TextLength, 0);//将插入符号置于文本结束处
MainForm.S.richTextBox1.ScrollToCaret();
}
}
}