Displaying 2 results from an estimated 2 matches for "p_1_2".
Did you mean:
p_12
2010 Jan 26
3
Problem with "nls" function
...c_1_F))+exp(gamma_2_F*stdS*(s_F-c_2_F)))/(1+exp(-gamma_1_F*stdS*(s_F-c_1_F))+exp(gamma_2_F*stdS*(s_F-c_2_F)))*D_3+a_4_2*(exp(-gamma_1_F*stdS*(s_F-c_1_F))+exp(gamma_2_F*stdS*(s_F-c_2_F)))/(1+exp(-gamma_1_F*stdS*(s_F-c_1_F))+exp(gamma_2_F*stdS*(s_F-c_2_F)))*D_4+p_1_1*Y_1+p_2_1*Y_2+p_3_1*Y_3+p_4_1*Y_4+p_1_2*(exp(-gamma_1_G*stdG*(s_G-c_1_G))+exp(gamma_2_G*stdG*(s_G-c_2_G)))/(1+exp(-gamma_1_G*stdG*(s_G-c_1_G))+exp(gamma_2_G*stdG*(s_G-c_2_G)))*Y_1+p_2_2*(exp(-gamma_1_G*stdG*(s_G-c_1_G))+exp(gamma_2_G*stdG*(s_G-c_2_G)))/(1+exp(-gamma_1_G*stdG*(s_G-c_1_G))+exp(gamma_2_G*stdG*(s_G-c_2_G)))*Y_2+p_3_2*(exp(-g...
2011 Nov 01
1
condition has length > 1 for LL denominator
...ed. I'm trying to
construct the denominator 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 <= departu...