Exemple #1
0
        internal bool Propose(out IterationState iterationState)
        {
            RequireFalse(IsComplete);

            if (_currentIteration.IsntFound())
            {
                for (_currentIteration = 0; _currentIteration < _iterations; ++_currentIteration)
                {
                    if (IsNaN(_errorData.Rows[_currentIteration].Field <double>(0)))
                    {
                        break;
                    }
                }
            }

            RequireFalse(IsComplete);

            if (_currentIteration == 0)
            {
                ProposeInitial(out iterationState);
            }
            else
            {
                ResumeIteration(out iterationState);
            }

            return(true);
        }
Exemple #2
0
        public void PartitionBooleanGrouping()
        {
            var(evens, odds) =
                Enumerable.Range(0, 10)
                .GroupBy(x => x % 2 == 0)
                .Partition((t, f) => Tuple.Create(t, f));

            Assert.That(evens, Is.EquivalentTo(new[] { 0, 2, 4, 6, 8 }));
            Assert.That(odds, Is.EquivalentTo(new[] { 1, 3, 5, 7, 9 }));
        }
Exemple #3
0
        public void PartitionNullableBooleanGrouping()
        {
            var xs = new int?[] { 1, 2, 3, null, 5, 6, 7, null, 9, 10 };

            var(lt5, gte5, nils) =
                xs.GroupBy(x => x != null ? x < 5 : (bool?)null)
                .Partition((t, f, n) => Tuple.Create(t, f, n));

            Assert.That(lt5, Is.EquivalentTo(new[] { 1, 2, 3 }));
            Assert.That(gte5, Is.EquivalentTo(new[] { 5, 6, 7, 9, 10 }));
            Assert.That(nils, Is.EquivalentTo(new int?[] { null, null }));
        }
Exemple #4
0
        internal bool Propose(
            NumDataColumn independentData,
            Arr <NumDataColumn> outputData,
            out IterationState iterationState
            )
        {
            RequireTrue(_currentIteration >= 0);
            RequireTrue(_currentIteration == 0 || !IsNaN(_currentLogLikelihood));
            RequireFalse(IsComplete);

            double proposedLogLikelihood;

            try
            {
                proposedLogLikelihood = GetProposalLogLikelihood(
                    independentData,
                    outputData
                    );
            }
            catch (Exception ex)
            {
                FaultChain(ex);
                iterationState = default;
                return(false);
            }

            ProcessCurrentProposal(proposedLogLikelihood); // M-H + bookkeeping (parameters)

            var currentChainDataRow = _chainData.Rows[_currentIteration];

            var indexNextModelParameter = _currentIteration == 0
        ? NOT_FOUND
        : _modelParameters.FindIndex(
                mp => IsNaN(currentChainDataRow.Field <double>(ToAcceptColumnName(mp.Name)))
                );

            var isIterationComplete = indexNextModelParameter.IsntFound();

            if (isIterationComplete)
            {
                StorePosterior();                // model outputs collected
            }
            var isLastIteration = (_currentIteration + 1) == _iterations;

            if (isIterationComplete && isLastIteration)
            {
                ++_currentIteration;
                var currentValues = _modelParameters
                                    .Select(mp => (mp.Name, Value: currentChainDataRow.Field <double>(mp.Name)))
                                    .ToArr();
                iterationState = (_currentIteration, currentValues, default);
                return(false);
            }

            if (isIterationComplete)
            {
                UpdateModelOutputs(); // M-H + bookkeeping (error models)

                var doAdjustForAcceptRate =
                    !IsNaN(_targetAcceptRate) &&
                    _currentIteration > _burnIn &&
                    _currentIteration % ACCEPT_RATE_ADJUST_INTERVAL == 0;

                if (doAdjustForAcceptRate)
                {
                    AdjustForAcceptRate();
                }

                // done - advance and initialize next iteration using current
                ++_currentIteration;
                indexNextModelParameter = 0;
                var nextChainDataRow = _chainData.Rows[_currentIteration];
                _modelParameters.Iter(
                    mp => nextChainDataRow[mp.Name] = currentChainDataRow[mp.Name]
                    );
                _modelParameters.Iter(
                    mp => nextChainDataRow[ToLLColumnName(mp.Name)] = currentChainDataRow[ToLLColumnName(mp.Name)]
                    );
                currentChainDataRow = nextChainDataRow;
            }

            _currentProposal = _modelParameters[indexNextModelParameter].GetProposal();

            var proposedValues = _modelParameters
                                 .Select(mp => (
                                             mp.Name,
                                             Value: mp.Name == _currentProposal.Name
            ? _currentProposal.Value
            : currentChainDataRow.Field <double>(mp.Name)
                                             )
                                         )
                                 .ToArr();

            if (isIterationComplete) // meaning was completed...
            {
                // ...so include update for client
                var currentValues = _modelParameters
                                    .Select(mp => (mp.Name, Value: currentChainDataRow.Field <double>(mp.Name)))
                                    .ToArr();
                iterationState = (_currentIteration, currentValues, proposedValues);
                return(true);
            }

            iterationState = (_currentIteration, default, proposedValues);
 ToValueTuple <T1, T2, T3, T4, T5, T6, T7, T8, T9>(
     this Tuple <T1, T2, T3, T4, T5, T6, T7, Tuple <T8, T9> > value)
 {
     return(CreateLong(value.Item1, value.Item2, value.Item3, value.Item4, value.Item5, value.Item6, value.Item7,
                       ValueTuple.Create(value.Rest.Item1, value.Rest.Item2)));
 }
 ToValueTuple <T1, T2, T3, T4, T5, T6, T7>(
     this Tuple <T1, T2, T3, T4, T5, T6, T7> value)
 {
     return(ValueTuple.Create(value.Item1, value.Item2, value.Item3, value.Item4, value.Item5, value.Item6, value.Item7));
 }
Exemple #7
0
 ToTuple <T1, T2, T3, T4, T5, T6>(
     this ValueTuple <T1, T2, T3, T4, T5, T6> value)
 {
     return(Tuple.Create(value.Item1, value.Item2, value.Item3, value.Item4, value.Item5, value.Item6));
 }
 public static Func <T1, T2, Unit> ToFunc <T1, T2>(this Action <T1, T2> action)
 => (T1 t1, T2 t2) =>
 {
     action(t1, t2);
     return(Unit.Create());
 };
Exemple #9
0
 public static Unit ForEach <T>(this List <T> @this, Action <T> action)
 {
     @this.Map(action.ToFunc());
     return(Unit.Create());
 }
 ToValueTuple <T1>(
     this Tuple <T1> value)
 {
     return(ValueTuple.Create(value.Item1));
 }
 /// <summary>Returns a value that indicates whether the current <see cref="T:System.ValueTuple`4" /> instance is equal to a specified <see cref="T:System.ValueTuple`4" /> instance.</summary>
 /// <param name="other">The value tuple to compare with this instance.</param>
 /// <returns>
 ///     <see langword="true" /> if the current instance is equal to the specified tuple; otherwise, <see langword="false" />.</returns>
 // Token: 0x06000429 RID: 1065 RVA: 0x0000B230 File Offset: 0x00009430
 public bool Equals(ValueTuple <T1, T2, T3, T4> other)
 {
     return(EqualityComparer <T1> .Default.Equals(this.Item1, other.Item1) && EqualityComparer <T2> .Default.Equals(this.Item2, other.Item2) && EqualityComparer <T3> .Default.Equals(this.Item3, other.Item3) && EqualityComparer <T4> .Default.Equals(this.Item4, other.Item4));
 }
 // Token: 0x06000430 RID: 1072 RVA: 0x0000B544 File Offset: 0x00009744
 private int GetHashCodeCore(IEqualityComparer comparer)
 {
     return(ValueTuple.CombineHashCodes(comparer.GetHashCode(this.Item1), comparer.GetHashCode(this.Item2), comparer.GetHashCode(this.Item3), comparer.GetHashCode(this.Item4)));
 }
 /// <summary>Calculates the hash code for the current <see cref="T:System.ValueTuple`4" /> instance.</summary>
 /// <returns>The hash code for the current <see cref="T:System.ValueTuple`4" /> instance. </returns>
 // Token: 0x0600042E RID: 1070 RVA: 0x0000B4E8 File Offset: 0x000096E8
 public override int GetHashCode()
 {
     return(ValueTuple.CombineHashCodes(EqualityComparer <T1> .Default.GetHashCode(this.Item1), EqualityComparer <T2> .Default.GetHashCode(this.Item2), EqualityComparer <T3> .Default.GetHashCode(this.Item3), EqualityComparer <T4> .Default.GetHashCode(this.Item4)));
 }
Exemple #14
0
 public bool Equals(ValueTuple <T1, T2, T3> obj)
 {
     return(_t1Comparer.Equals(Item1, obj.Item1) &&
            _t2Comparer.Equals(Item2, obj.Item2) &&
            _t3Comparer.Equals(Item3, obj.Item3));
 }
Exemple #15
0
 ToTuple <T1, T2, T3, T4, T5, T6, T7, T8, T9, T10, T11, T12, T13, T14>(
     this ValueTuple <T1, T2, T3, T4, T5, T6, T7, ValueTuple <T8, T9, T10, T11, T12, T13, T14> > value)
 {
     return(CreateLongRef(value.Item1, value.Item2, value.Item3, value.Item4, value.Item5, value.Item6, value.Item7,
                          Tuple.Create(value.Rest.Item1, value.Rest.Item2, value.Rest.Item3, value.Rest.Item4, value.Rest.Item5, value.Rest.Item6, value.Rest.Item7)));
 }
 ToValueTuple <T1, T2>(
     this Tuple <T1, T2> value)
 {
     return(ValueTuple.Create(value.Item1, value.Item2));
 }
 ToValueTuple <T1, T2, T3>(
     this Tuple <T1, T2, T3> value)
 {
     return(ValueTuple.Create(value.Item1, value.Item2, value.Item3));
 }
 public static Func <T, Unit> ToFunc <T>(this Action <T> action)
 => t =>
 {
     action(t);
     return(Unit.Create());
 };
 ToValueTuple <T1, T2, T3, T4, T5>(
     this Tuple <T1, T2, T3, T4, T5> value)
 {
     return(ValueTuple.Create(value.Item1, value.Item2, value.Item3, value.Item4, value.Item5));
 }
 public static Func <Unit> ToFunc(this Action action)
 => () =>
 {
     action();
     return(Unit.Create());
 };
Exemple #21
0
 ToTuple <T1, T2, T3, T4>(
     this ValueTuple <T1, T2, T3, T4> value)
 {
     return(Tuple.Create(value.Item1, value.Item2, value.Item3, value.Item4));
 }