r - How to change na.action for zero-inflated regression model? -

-

i running zero-inflated negative binomial regression model using function zeroinfl pscl package.

i need exclude na's model in order able plot residuals against dependent variable later in analysis.

therefore, want set na.action="na.exclude". can without problem non-zero-inflated negative binomial regression model (using glm.nb glm package), eg. 

fm_nbin <- glm.nb(dv ~ factor(idv) + contr1                +contr2 + contr3, data=df,                 subset=(df$var<500), na.action="na.exclude") fm_nbin.res = resid(fm_nbin)  plot(fm_nbin.res~df$var)

works fine. however, when same zero-inflated model, not work: 

zinfl <- zeroinfl(dv ~ factor(idv) + contr1                +contr2 + contr3 | factor(idv) + contr1                +contr2 + contr3, data=df,                 subset=(df$var<500), na.action="na.exclude") zinfl.res = resid(zinfl)  plot(zinfl.res~df$var)

gives error 

error in function (formula, data = null, subset = null, na.action = na.fail,  :    variable lengths differ (found 'df$var')

is there other command should use exclude na's regression?

edit: this nearest of answer find. can in way applied problem? also, can naresid in way applied?

as 1 finds following trail of documentation zeroinfl glm.fit: "the ‘factory-fresh’ default na.omit." notice have not put quotes around since supposed function rather function accept name doesn't matter if quoted. admit don't know how na.omit , na.exclude differ (something residuals read), go default setting first, since delivers want regression functions. try leaving out:

zinfl <- zeroinfl(dv ~ factor(idv) + contr1            +contr2 + contr3 | factor(idv) + contr1            +contr2 + contr3, data=df,             subset=(df$var<500) )

Comments

Post a Comment

Popular posts from this blog

c++ - Function signature as a function template parameter -

-
i'd use function signature template argument. works great classes, when try same trick function templates, msvc throws error: error c2768: 'func' : illegal use of explicit template arguments here's code: template <typename signature> void func(); template <typename r, typename a1> void func<r(a1)>(); what should make work? you cannot partially specialize function template, not supported language. can create partially specialized class template static member function, , possibly trampoline function instantiate class template , invoke static function. something this: namespace detail { template<typename signature> struct helper; template<typename r, typename a1> struct helper<r(a1)> { static void call() { // stuff r , a1... } }; } template<typename signature> void func() { detail::helper<signature>::call(); }
Read more

algorithm - What are some ways to combine a number of (potentially incompatible) sorted sub-sets of a total set into a (partial) ordering of the total set? -

-
i have sorting problem i'm pretty sure doesn't have canonical "answer" - has number of different approaches, each pros/cons. i'm interested in hearing different approaches. suppose bunch of people {p_i | i=1,...,m} ranking favorite ice cream flavors. in total there n different flavors, single person p_i ranks n_i << n flavors or familiar with. i'm interested in combining these sub-rankings plausible overall ranking of n flavors. for concrete situation: m , n both 1000 (by coincidence), , each n_i 20 . can assume there's enough overlap in flavors people rank such no flavor "isolated". again, i'm interested in hearing different ways approach this, if there isn't single clear answer. thanks! one way approach problem construct directed cyclic graph representing of constraints on ordering of elements. nodes in graph represent different elements compared , have edge first node second node if had ranked first e...
Read more

How to call a javascript function after the page loads with a chrome extension? -

-
there's function on page want call after loads. not own page, can't insert html it. can please provide specific examples how chrome extension or greasemonkey? (function submitaction() ). both chrome extensions , greasemonkey offer ways specify when script inject if executed. in chrome extension, you set in manifest : "content_scripts": [ { "run_at" : "document_idle", in greasemonkey compatible userscripts, set in metadata // @run-at document-end in both cases, have add script if want ensure resources loaded when function executed : document.addeventlistener('load', submitaction);
Read more