Displaying 1 result from an estimated 1 matches for "get_denominator".
2011 Nov 01
1
condition has length > 1 for LL denominator
...tor for a likelihood function using the following
function. For the first row in "results", for example, I want the
denominator to be the sum of all possible arrive/depart combinations an
interceptor in period 1 could observe: exp(P_1_1) + exp(P_1_2) + exp(P_1_3)
(i.e. P_arrive_depart).
get_denominator = function(intercept, periods_per_day)
{
denominator = array("(", nrow(results))
for (arrival in 1:periods_per_day)
{
for (departure in arrival:periods_per_day)
{
while (arrival <= intercept & intercept <= departure)
{
addition_to_denom = paste("P", arr...