From: unknown <hskwk@.(none)>
Date: Wed, 18 Aug 2010 08:33:17 +0000 (+0900)
Subject: solver2
X-Git-Url: http://git.osdn.jp/view?a=commitdiff_plain;h=1b8c92d90625dad0dcee96ad2db68d47b4efe73f;p=psychlops%2Fsilverlight.git

solver2
---

diff --git a/dev4/Pages/BinomialSolver.xaml b/dev4/Pages/BinomialSolver.xaml
new file mode 100644
index 0000000..c2be161
--- /dev/null
+++ b/dev4/Pages/BinomialSolver.xaml
@@ -0,0 +1,108 @@
+<controls:ChildWindow x:Class="PsychlopsSilverlight4.Pages.BinomialSolver"
+           xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation" 
+           xmlns:x="http://schemas.microsoft.com/winfx/2006/xaml" 
+           xmlns:sdk="http://schemas.microsoft.com/winfx/2006/xaml/presentation/sdk"
+           xmlns:controls="clr-namespace:System.Windows.Controls;assembly=System.Windows.Controls"
+           Width="400" Height="300" 
+           Title="Binomial Likelihood Solver" xmlns:toolkit="http://schemas.microsoft.com/winfx/2006/xaml/presentation/toolkit">
+
+
+	<controls:ChildWindow.Resources>
+		<Style TargetType="toolkit:Chart">
+			<!-- To Remove Legend-->
+			<Setter Property="Padding" Value="10 0" />
+			<Setter Property="BorderBrush" Value="Transparent" />
+			<Setter Property="PlotAreaStyle">
+				<Setter.Value>
+					<Style TargetType="Grid">
+						<Setter Property="Margin" Value="0"/>
+					</Style>
+				</Setter.Value>
+			</Setter>
+			<Setter Property="TitleStyle">
+				<Setter.Value>
+					<Style TargetType="toolkit:Title">
+						<Setter Property="Width" Value="0"/>
+						<Setter Property="Height" Value="0"/>
+					</Style>
+				</Setter.Value>
+			</Setter>
+			<Setter Property="LegendStyle">
+				<Setter.Value>
+					<Style TargetType="toolkit:Legend">
+						<Setter Property="Width" Value="0"/>
+						<Setter Property="Height" Value="0"/>
+					</Style>
+				</Setter.Value>
+			</Setter>
+		</Style>
+	</controls:ChildWindow.Resources>
+
+	<Grid x:Name="LayoutRoot" Margin="2">
+        <Grid.RowDefinitions>
+            <RowDefinition />
+            <RowDefinition Height="Auto" />
+        </Grid.RowDefinitions>
+
+		<StackPanel Orientation="Horizontal">
+			<Canvas Width="250" Height="Auto">
+				<toolkit:Chart x:Name="chart1" Title="" Width="240" Height="220" Canvas.Left="0" Canvas.Top="0"
+							   FontSize="10">
+					<toolkit:Chart.Series>
+						<toolkit:ScatterSeries x:Name="Fitted"
+										Title="Fitted"
+										ItemsSource="{Binding}"
+										IndependentValueBinding="{Binding c_abscissa}"
+										DependentValueBinding="{Binding c_ratio}"
+										>
+
+							<toolkit:ScatterSeries.DependentRangeAxis>
+								<toolkit:LinearAxis Title="" ShowGridLines="False" Orientation="Y" Minimum="0" Maximum="1">
+								</toolkit:LinearAxis>
+							</toolkit:ScatterSeries.DependentRangeAxis>
+						</toolkit:ScatterSeries>
+						<toolkit:ScatterSeries x:Name="Raw"
+										Title="Fitted"
+										ItemsSource="{Binding}"
+										IndependentValueBinding="{Binding abscissa}"
+										DependentValueBinding="{Binding ratio}"
+										>
+							<!--
+						<toolkit:ScatterSeries.IndependentAxis>
+							<toolkit:LinearAxis Title="" Orientation="X" Maximum="10" Minimum="0">
+							</toolkit:LinearAxis>
+						</toolkit:ScatterSeries.IndependentAxis>
+						-->
+
+							<toolkit:ScatterSeries.DependentRangeAxis>
+								<toolkit:LinearAxis Title="" ShowGridLines="False" Orientation="Y" Minimum="0" Maximum="1">
+								</toolkit:LinearAxis>
+							</toolkit:ScatterSeries.DependentRangeAxis>
+						</toolkit:ScatterSeries>
+					</toolkit:Chart.Series>
+				</toolkit:Chart>
+				<Slider x:Name="OrdinateSlider" Canvas.Left="0" Canvas.Top="10" Orientation="Vertical"  Height="200"></Slider>
+				<Slider x:Name="AbscissaSlider" Canvas.Left="20" Canvas.Top="210" Width="200" MouseLeftButtonUp="AbscissaSlider_MouseLeftButtonUp"></Slider>
+			</Canvas>
+
+
+			<StackPanel Orientation="Vertical">
+				<sdk:DataGrid x:Name="dataGrid" Margin="0,5,0,10" RowHeight="40" AutoGenerateColumns="False" >
+					<sdk:DataGrid.Columns>
+						<sdk:DataGridTextColumn Header="Abscissa" Width="40" Binding="{Binding abscissa}" />
+						<sdk:DataGridTextColumn Header="Positive Cases" Width="40"  Binding="{Binding positive}" />
+						<sdk:DataGridTextColumn Header="Negative Cases" Width="40" Binding="{Binding negative}" />
+					</sdk:DataGrid.Columns>
+				</sdk:DataGrid>
+				<Button x:Name="AddConditionButton" Content="Add Condition" Click="AddConditionButton_Click" HorizontalAlignment="Right" FontSize="10" Margin="0,12,0,0" Grid.Row="1" />
+			</StackPanel>
+		</StackPanel>
+
+		<TextBlock Text="Mean" x:Name="Parameter0Name" Width="40" Height="23" HorizontalAlignment="Left" Margin="0,12,0,0" Grid.Row="1" />
+		<TextBox x:Name="Parameter0" IsReadOnly="True" Width="50" Height="23" HorizontalAlignment="Left" Margin="40,12,0,0" Grid.Row="1" />
+		<TextBlock Text="S.D." x:Name="Parameter1Name" Width="40" Height="23" HorizontalAlignment="Left" Margin="110,12,0,0" Grid.Row="1" />
+		<TextBox x:Name="Parameter1" IsReadOnly="True" Width="50" Height="23" HorizontalAlignment="Left" Margin="150,12,0,0" Grid.Row="1" />
+		<Button x:Name="CalcButton" Content="Calculate" Click="CalcButton_Click" Width="75" Height="23" HorizontalAlignment="Right" Margin="0,12,79,0" Grid.Row="1" />
+        <Button x:Name="OKButton" Content="OK" Click="OKButton_Click" Width="75" Height="23" HorizontalAlignment="Right" Margin="0,12,0,0" Grid.Row="1" />
+    </Grid>
+</controls:ChildWindow>
\ No newline at end of file
diff --git a/dev4/Pages/BinomialSolver.xaml.cs b/dev4/Pages/BinomialSolver.xaml.cs
new file mode 100644
index 0000000..d589370
--- /dev/null
+++ b/dev4/Pages/BinomialSolver.xaml.cs
@@ -0,0 +1,216 @@
+using System;
+using System.Collections.Generic;
+using System.Collections.ObjectModel;
+using System.Linq;
+using System.Net;
+using System.Windows;
+using System.Windows.Controls;
+using System.Windows.Documents;
+using System.Windows.Input;
+using System.Windows.Media;
+using System.Windows.Media.Animation;
+using System.Windows.Shapes;
+
+namespace PsychlopsSilverlight4.Pages
+{
+	public partial class BinomialSolver : ChildWindow
+	{
+		int func_selected;
+		Psychlops.Solver.Constants.Func1 func1;
+		Psychlops.Solver.Constants.Func2 func2;
+		Psychlops.Solver.Constants.Func3 func3;
+
+
+		struct WholeData
+		{
+			public ObservableCollection<Data> data {get; set;}
+			public ObservableCollection<Data> fitted {get; set;}
+		}
+		WholeData data;
+
+		public BinomialSolver(Psychlops.Solver.Constants.Func1 f)
+		{
+			InitializeComponent();
+			initialize();
+			func_selected = 1;
+			func1 = f;
+		}
+		public BinomialSolver(Psychlops.Solver.Constants.Func2 f)
+		{
+			InitializeComponent();
+			initialize();
+			func_selected = 2;
+			func2 = f;
+		}
+		public BinomialSolver(Psychlops.Solver.Constants.Func3 f)
+		{
+			InitializeComponent();
+			initialize();
+			func_selected = 3;
+			func3 = f;
+		}
+		public void initialize()
+		{
+			dataGrid.RowHeight = 20;
+			data.data = Data.GetDataList();
+			dataGrid.ItemsSource = data.data;
+
+			chart1.DataContext = data.data;
+			solver = new Psychlops.Solver.BinomialLikelihood();
+		}
+		public class Data
+		{
+			public double c_abscissa { get; set; }
+			public double c_ratio { get; set; }
+			public double abscissa { get; set; }
+			public int positive { get; set; }
+			public int negative { get; set; }
+			public double ratio { get { return (double)(positive) / (positive+negative); } }
+
+			// A null value for IsSubscribed can indicate 
+			// "no preference" or "no response".
+			public Boolean? IsSubscribed { get; set; }
+
+			public Data(double absc, int posi, int nega)
+			{
+				this.abscissa = absc;
+				this.positive = posi;
+				this.negative = nega;
+				c_abscissa = 0.0;
+				c_ratio = 0.0;
+			}
+
+			public static ObservableCollection<Data> GetDataList()
+			{
+				return new ObservableCollection<Data>(new Data[] {
+					new Data(-5, 0, 10),
+					new Data(0, 1, 9),
+					new Data(1, 3, 7),
+					new Data(2, 5, 5),
+					new Data(3, 7, 3),               
+					new Data(5, 9, 1),                
+					new Data(10, 10, 0)                
+				});
+			}
+		}
+
+		private void OKButton_Click(object sender, RoutedEventArgs e)
+		{
+			this.DialogResult = true;
+		}
+
+		double minimum, maximum;
+		Psychlops.Solver.BinomialLikelihood solver;
+		private void CalcButton_Click(object sender, RoutedEventArgs e)
+		{
+			var ber = new Psychlops.Solver.BernoulliProcess();
+			ber.elems = new Psychlops.Solver.BernoulliProcess.Data[data.data.Count];
+			int i = 0;
+
+			minimum = Double.PositiveInfinity;
+			maximum = Double.NegativeInfinity;
+			foreach (var d in data.data)
+			{
+				if (d.abscissa < minimum) { minimum = d.abscissa; }
+				if (d.abscissa > maximum) { maximum = d.abscissa; }
+				ber.elems[i].x = d.abscissa;
+				ber.elems[i].pos = d.positive;
+				ber.elems[i].neg = d.negative;
+				i++;
+			}
+
+			solver.data = ber;
+
+
+			Psychlops.Interval itvl = new Psychlops.Interval();
+			solver.itvl[0] = (-100 <= itvl <= 100);
+			solver.itvl[1] = (-100 <= itvl <= 100);
+			solver.itvl[2] = (-100 <= itvl <= 100);
+
+			switch (func_selected)
+			{
+				case 1:
+					solver.begin(func1);
+					break;
+				case 2:
+					solver.begin(func2);
+					break;
+				case 3:
+					solver.begin(func3);
+					break;
+			}
+
+			reFunc();
+		}
+		private void reFunc()
+		{
+			data.fitted = new ObservableCollection<Data>();
+			foreach (var d in data.data)
+			{
+				data.fitted.Add(d);
+			}
+			Data dd;
+			switch (func_selected)
+			{
+				case 1:
+					for (double j = minimum; j < maximum; j += (maximum - minimum)/ 100 )
+					{
+						dd = new Data(0, 0, 0);
+						dd.c_abscissa = j;
+						dd.c_ratio = func1(j, solver.champ[0]);
+						data.fitted.Add(dd);
+					}
+					break;
+				case 2:
+					for (double j = minimum; j < maximum; j += (maximum - minimum)/ 100 )
+					{
+						dd = new Data(0, 0, 0);
+						dd.c_abscissa = j;
+						dd.c_ratio = func2(j, solver.champ[0], solver.champ[1]);
+						data.fitted.Add(dd);
+					}
+					break;
+				case 3:
+					for (double j = minimum; j < maximum; j += (maximum - minimum)/ 100 )
+					{
+						dd = new Data(0, 0, 0);
+						dd.c_abscissa = j;
+						dd.c_ratio = func3(j, solver.champ[0], solver.champ[1], solver.champ[2]);
+						data.fitted.Add(dd);
+					}
+					break;
+			}
+
+
+			chart1.DataContext = null;
+			chart1.DataContext = data.fitted;
+
+			OrdinateSlider.Maximum = 1.0;
+			OrdinateSlider.Minimum = 0.0;
+			OrdinateSlider.Value = 0.5;
+
+			AbscissaSlider.Maximum = maximum;
+			AbscissaSlider.Minimum = minimum;
+			AbscissaSlider.Value = solver.champ[0];
+
+			Parameter0.Text = String.Format("{0:D4}", solver.champ[0].ToString());
+			Parameter1.Text = String.Format("{0, 4:F}", solver.champ[1].ToString());
+
+			//chart1.UpdateLayout();
+		}
+
+		private void AddConditionButton_Click(object sender, RoutedEventArgs e)
+		{
+			data.data.Add(new Data(0, 0, 0));
+			this.UpdateLayout();
+		}
+
+		private void AbscissaSlider_MouseLeftButtonUp(object sender, MouseButtonEventArgs e)
+		{
+			solver.champ[0] = AbscissaSlider.Value;
+			reFunc();
+		}
+
+	}
+}
+