Merge branch 'master' of https://github.com/xieguigang/sciBASIC

Author: xieguigang

Data_science/Mathematica/Math/Math.Statistics/ShapleyValue/CharacteristicFunction.vb

@@ -0,0 +1,98 @@
+Imports Microsoft.VisualBasic.ComponentModel.Collection
+
+Namespace ShapleyValue
+
+    ''' <summary>
+    ''' represent the following function often written by v
+    ''' 2^N to R
+    ''' for each subset of the set {1.. N} is associated a value
+    ''' per convention v(empty set) = 0
+    ''' 
+    ''' @author Franck Benault
+    ''' 
+    ''' </summary>
+    Public Class CharacteristicFunction
+
+        Dim v As IDictionary(Of ISet(Of Integer), Double)
+
+        Public Overridable ReadOnly Property NbPlayers As Integer
+
+        Private Sub New(builder As CharacteristicFunctionBuilder)
+            NbPlayers = builder.nbPlayers
+            v = New Dictionary(Of ISet(Of Integer), Double)()
+            v = builder.v
+        End Sub
+
+        Public Overridable Function getValue(ParamArray coalition As Integer?()) As Double
+            Dim coalitionSet As ISet(Of Integer) = New HashSet(Of Integer)()
+
+            For Each player As Integer In coalition
+                coalitionSet.Add(player)
+            Next
+
+            Return v(coalitionSet)
+        End Function
+
+        Public Overridable Function getValue(coalitionSet As ISet(Of Integer)) As Double
+
+            Return v(coalitionSet)
+        End Function
+
+        Public Overrides Function ToString() As String
+            Return "CharacteristicFunction [nbPlayers=" & NbPlayers.ToString() & ", v=" & v.ToString() & "]"
+        End Function
+
+        Public Class CharacteristicFunctionBuilder
+            Friend ReadOnly nbPlayers As Integer
+            Friend v As IDictionary(Of ISet(Of Integer), Double)
+
+            Public Sub New(nbPlayers As Integer)
+                Me.nbPlayers = nbPlayers
+                v = New Dictionary(Of ISet(Of Integer), Double)()
+            End Sub
+
+            Public Overridable Function addCoalition(value As Double, ParamArray coalition As Integer?()) As CharacteristicFunctionBuilder
+                Dim [set] As ISet(Of Integer) = New HashSet(Of Integer)()
+                For Each player As Integer In coalition
+                    [set].Add(player)
+                Next
+                v([set]) = value
+                Return Me
+            End Function
+
+            Public Overridable Function build() As CharacteristicFunction
+                Return New CharacteristicFunction(Me)
+            End Function
+
+        End Class
+
+        Public Overridable Sub addDummyUser()
+            Dim newV As IDictionary(Of ISet(Of Integer), Double) = New Dictionary(Of ISet(Of Integer), Double)()
+
+
+            Dim coalitions = v.Keys
+            Dim coalitionSet As HashSet(Of Integer)
+
+            For Each coalition In coalitions
+
+                coalitionSet = New HashSet(Of Integer)()
+                coalitionSet.AddAll(coalition)
+                coalitionSet.Add(NbPlayers + 1)
+
+                newV(coalitionSet) = v(coalition)
+            Next
+
+            coalitionSet = New HashSet(Of Integer)()
+            coalitionSet.Add(NbPlayers + 1)
+            v(coalitionSet) = 0.0
+
+            For Each coalition In newV.Keys
+                v(coalition) = newV(coalition)
+            Next
+
+            _NbPlayers += 1
+        End Sub
+
+    End Class
+
+End Namespace

Data_science/Mathematica/Math/Math.Statistics/ShapleyValue/CoalitionStrategy.vb

@@ -0,0 +1,71 @@
+Imports std = System.Math
+
+Namespace ShapleyValue
+
+    ''' 
+    ''' <summary>
+    ''' @author Franck Benault
+    ''' 
+    ''' @version	0.0.2
+    ''' @since 0.0.2
+    ''' 
+    ''' </summary>
+    Public NotInheritable Class CoalitionStrategy
+
+        Public Shared ReadOnly SEQUENTIALField As CoalitionStrategy = New CoalitionStrategy("SEQUENTIAL", InnerEnum.SEQUENTIAL)
+
+        Public Shared ReadOnly RANDOM As CoalitionStrategy = New CoalitionStrategy("RANDOM", InnerEnum.RANDOM)
+
+        Private Shared ReadOnly valueList As List(Of CoalitionStrategy) = New List(Of CoalitionStrategy)()
+
+        Shared Sub New()
+            valueList.Add(SEQUENTIALField)
+            valueList.Add(RANDOM)
+        End Sub
+
+        Public Enum InnerEnum
+            SEQUENTIAL
+            RANDOM
+        End Enum
+
+        Public ReadOnly innerEnumValue As InnerEnum
+        Private ReadOnly nameValue As String
+        Private ReadOnly ordinalValue As Integer
+        Private Shared nextOrdinal As Integer = 0
+
+        Private Sub New(name As String, innerEnum As InnerEnum)
+            nameValue = name
+            ordinalValue = nextOrdinal
+            nextOrdinal += 1
+            innerEnumValue = innerEnum
+        End Sub
+
+        Public ReadOnly Property Sequential As Boolean
+            Get
+                Return Equals(SEQUENTIALField)
+            End Get
+        End Property
+
+        Public Shared Function values() As CoalitionStrategy()
+            Return valueList.ToArray()
+        End Function
+
+        Public Function ordinal() As Integer
+            Return ordinalValue
+        End Function
+
+        Public Overrides Function ToString() As String
+            Return nameValue
+        End Function
+
+        Public Shared Function valueOf(name As String) As CoalitionStrategy
+            For Each enumInstance In valueList
+                If Equals(enumInstance.nameValue, name) Then
+                    Return enumInstance
+                End If
+            Next
+            Throw New ArgumentException(name)
+        End Function
+    End Class
+
+End Namespace

Data_science/Mathematica/Math/Math.Statistics/ShapleyValue/FactorialUtil.vb

@@ -0,0 +1,16 @@
+Namespace ShapleyValue
+
+    Public Class FactorialUtil
+
+        Public Shared Function factorial(n As Long) As Long
+            If n > 20 Then
+                Throw New ArithmeticException("Capacity exceded for factorial " & n.ToString())
+            End If
+            Return If(n > 1, n * factorial(n - 1), 1)
+        End Function
+
+
+
+    End Class
+
+End Namespace

Data_science/Mathematica/Math/Math.Statistics/ShapleyValue/Node.vb

@@ -0,0 +1,59 @@
+
+Namespace ShapleyValue
+
+    Public Class Node
+
+        Public Overrides Function ToString() As String
+            Return "Node [value=" & valueField.ToString() & ", nextNode=" & nextNodeField.ToString() & "]"
+        End Function
+
+        Private valueField As NodeValue
+
+        Private nextNodeField As Node
+        Private prevNodeField As Node
+
+        Public Sub New(nodeValue As NodeValue, prevNode As Node)
+            valueField = nodeValue
+            prevNodeField = prevNode
+
+            If valueField.NextValues.Count > 0 Then
+                Dim values = valueField.NextValues
+                Dim nextValue As NodeValue = New NodeValue(values)
+                nextNodeField = New Node(nextValue, Me)
+            End If
+        End Sub
+
+        Public Overridable ReadOnly Property Value As NodeValue
+            Get
+                Return valueField
+            End Get
+        End Property
+
+        Public Overridable Property NextNode As Node
+            Get
+                Return nextNodeField
+            End Get
+            Set(value As Node)
+                nextNodeField = value
+
+            End Set
+        End Property
+
+        Public Overridable Sub updateValue()
+            valueField.updateValue()
+        End Sub
+
+
+        Public Overridable Property PrevNode As Node
+            Get
+                Return prevNodeField
+            End Get
+            Set(value As Node)
+                prevNodeField = value
+            End Set
+        End Property
+
+
+    End Class
+
+End Namespace

Data_science/Mathematica/Math/Math.Statistics/ShapleyValue/NodeValue.vb

@@ -0,0 +1,45 @@
+Namespace ShapleyValue
+
+    Public Class NodeValue
+
+        Public Overrides Function ToString() As String
+            Return "NodeValue [value=" & valueField.ToString() & ", nextValues=" & nextValuesField.ToString() & "]"
+        End Function
+
+        Private valueField As Integer
+
+        Private nextValuesField As IList(Of Integer)
+
+        Public Sub New(nextValues As IList(Of Integer))
+
+            nextValuesField = New List(Of Integer)()
+            CType(nextValuesField, List(Of Integer)).AddRange(nextValues)
+            valueField = nextValues(0)
+            nextValuesField.RemoveAt(valueField)
+        End Sub
+
+        Public Overridable ReadOnly Property Value As Integer
+            Get
+                Return valueField
+            End Get
+        End Property
+
+        Public Overridable ReadOnly Property NextValues As IList(Of Integer)
+            Get
+                Return New List(Of Integer)(nextValuesField)
+            End Get
+        End Property
+
+        Public Overridable Sub updateValue()
+            If nextValuesField.Count > 0 Then
+                valueField = nextValuesField(0)
+            End If
+            nextValuesField.RemoveAt(valueField)
+
+        End Sub
+
+
+
+    End Class
+
+End Namespace