CoinFlipModel.cpp

#include "CoinFlipModel.hpp"
#include "ParameterDouble.hpp"
#include "Probability.hpp"
#include "Parameter.hpp"
#include "RandomVariable.hpp"
#include <iostream>
CoinFlipModel::CoinFlipModel(int numheads, int numtails): Model(), numberheads(numheads), numbertails(numtails){
    bias = new ParameterDouble(1.0, this, "bias", 0, 1);
    parameters.push_back(bias);
}
double CoinFlipModel::lnlikelihood(void) {
    return Probability::Binom::lnpmf(bias->getValue(), numberheads+numbertails, numberheads);
    
}
// CoinFlipModel (int numheads, int numtails);
//    double lnlikelihood (void);
double CoinFlipModel::lnpriorprobability (void){
    double lnpro = 0.0;
    for (Parameter* p:parameters){
        lnpro+=p->lnProbability();
    }
    return lnpro;
}

double CoinFlipModel::update (void){
    RandomVariable& rng = RandomVariable::randomVariableInstance();
    UpdatedParameter = parameters [(int)(rng.uniformRv()*parameters.size())];
    return UpdatedParameter -> update();
}
void CoinFlipModel::updateforacceptance (void){
    UpdatedParameter->UpdateForAcceptance();
}
void CoinFlipModel::updateforrejection (void){
    UpdatedParameter->UpdateForRejection();
}
void CoinFlipModel::print (void){
    std::cout<<"Not Implemented";
}

std::vector <std::string> CoinFlipModel::GetParameterNames (void){
    std::vector<std::string> names;
    for (Parameter*p:parameters){
        names.push_back(p->getName());
    }
    return names;
}
std::vector <double> CoinFlipModel::GetParameterString (void){
    std::vector<double> value;
    for (Parameter*p:parameters){
        ParameterDouble* pd = dynamic_cast<ParameterDouble*>(p);
        if (pd != nullptr)
            value.push_back(pd->getValue());
    }
    return value;
}