Displaying 20 results from an estimated 10000 matches similar to: "interval regression"
2018 Mar 19
4
Struggling to compute marginal effects !
Dear Oscar,
and any other R-project person,
Can you please help me to figure out the meaning of the following error
message in red ?
Error in eval(predvars, data, env) :
numeric 'envir' arg not of length one
I computed ordered logit models using 'polr' in R (I just followed the
guidance a handout I found on princeton.edu about logit, probit and
multinomial logit models) . The
2003 Oct 11
1
boot statictic fn for dual estimation of 2 stats?
Hi,
I am trying to use boot() to refit an ordinal logit (polr in MASS) model.
(A very basic bootstrap which samples from the data frame without
replacement and updates the model.)
I need to extract two statistics per run (the coefficients and zeta) and I
tried concatenating them into a single vector after fitting, but I get the
following error:
Error in "[<-"(*tmp*, r, ,
2018 Mar 20
0
Struggling to compute marginal effects !
In that case, I can't work out why the first model fails but not the
second. I would start looking at "Data" to see what it contains. if:
object2 <- polr(Inc ~ Training ,Data,Hess = T,method = "logistic" )
works, the problem may be with the "Adopt" variable.
Jim
On Tue, Mar 20, 2018 at 10:55 AM, Willy Byamungu
<wmulimbi at email.uark.edu> wrote:
>
2004 Mar 24
2
Ordered logit/probit
Hello everyone
I am trying to fit an ordered probit/logit model for bank rating
prediction.
Besides polr() in MASS package which is not written especially for this as
far as I know, do you know how else I can do this?
I already found the modified polr () version on the
Valentin STANESCU
Enrst and Young
Tel. 402 4000
----------------------------------------------------------
The information
2010 Nov 03
2
bugs and misfeatures in polr(MASS).... fixed!
In polr.R the (several) functions gmin and fmin contain the code
> theta <- beta[pc + 1L:q]
> gamm <- c(-100, cumsum(c(theta[1L], exp(theta[-1L]))), 100)
That's bad. There's no reason to suppose beta[pc+1L] is larger than
-100 or that the cumulative sum is smaller than 100. For practical
datasets those assumptions are frequently violated, causing the
2005 Jun 10
1
problem with polr ?
I want to fit a multinomial model with logit link.
For example let this matrix to be analyzed:
male female aborted factor
10 12 1 1.2
14 14 4 1.3
15 12 3 1.4
(this is an example, not the true data which are far more complex...)
I suppose the correct function to analyze these data is polr from MASS library.
The data have been
2003 Feb 19
1
Multiple Logit/Probit
Does anybody know how to do multiple logit/probit analysis with R?
Thanks in advance!
2005 Nov 21
2
Multinomial Nested Logit package in R?
Dear R-Help,
I'm hoping to find a Multinomial Nested Logit package in R. It would
be great to find something analogous to "PROC MDC" in SAS:
> The MDC (Multinomial Discrete Choice) procedure analyzes models
> where the
> choice set consists of multiple alternatives. This procedure
> supports conditional logit,
> mixed logit, heteroscedastic extreme value,
2004 Dec 03
3
multinomial probit
Hello All,
I'm trying to run a multinomial probit on a dataset with 28 data
points and five levels (0,1,2,3,4) in the latent choice involving
response variable.
I downloaded the latest mnp package to run the regression. It starts
the calculation and then crashes the rpogram. I wish I could give the
error message but it literally shuts down R without a warning.
I'm using the R
2004 Jun 12
2
ordered probit or logit / recursive regression
> I make a study in health econometrics and have a categorical
> dependent variable (take value 1-5). I would like to fit an ordered
> probit or ordered logit but i didn't find a command or package who
> make that. Does anyone know if it's exists ?
R is very fancy. You won't get mundane things like ordered probit off
the shelf. (I will be very happy if someone will show
2011 Feb 16
1
error in optim, within polr(): "initial value in 'vmmin' is not finite"
Hi all. I'm just starting to explore ordinal multinomial regression. My dataset is 300,000 rows, with an outcome (ordinal factor from 1 to 9) and five independent variables (all continuous). My first stab at it was this:
pomod <- polr(Npf ~ o_stddev + o_skewness + o_kurtosis + o_acl_1e + dispersal, rlc, Hess=TRUE)
And that worked; I got a good model fit. However, a variety of other
2007 Apr 18
3
Problems in programming a simple likelihood
As part of carrying out a complicated maximum likelihood estimation, I
am trying to learn to program likelihoods in R. I started with a simple
probit model but am unable to get the code to work. Any help or
suggestions are most welcome. I give my code below:
************************************
mlogl <- function(mu, y, X) {
n <- nrow(X)
zeta <- X%*%mu
llik <- 0
for (i in 1:n) {
if
2004 Nov 11
1
polr probit versus stata oprobit
Dear All,
I have been struggling to understand why for the housing data in MASS
library R and stata give coef. estimates that are really different. I also
tried to come up with many many examples myself (see below, of course I
did not have the set.seed command included) and all of my
`random' examples seem to give verry similar output. For the housing data,
I have changed the data into numeric
2002 May 03
3
Regression models for ordinal responses ??
Hello list,
Is there any mean to fit models for ordinal response other than multinomial
polytomous ("multinom" from nnet ) and cumulative logit ("polr" from MASS)?
I am particularly interested in continuation-ratio model and
adjacent-category logit model. It is for the sake of epidemiology in
wild-living populations!
Many thanks,
Emmanuelle Fromont
2006 Aug 17
1
Setting contrasts for polr() to get same result of SAS
Hi all,
I am trying to do a ordered probit regression using polr(), replicating a
result from SAS.
>polr(y ~ x, dat, method='probit')
suppose the model is y ~ x, where y is a factor with 3 levels and x is a
factor with 5 levels,
To get coefficients, SAS by default use the last level as reference, R by
default use the first level (correct me if I was wrong),
The result I got is a
2007 Jun 04
2
How to obtain coefficient standard error from the result of polr?
Hi - I am using polr. I can get a result from polr fit by calling
result.plr <- polr(formula, data=mydata, method="probit");
However, from the 'result.plr', how can I access standard error of the estimated coefficients as well as the t statistics for each one of them?
What I would like to do ultimately is to see which coefficients are not significant and try to refit the
2008 Mar 15
1
again with polr
hello everybody
solved the problem with summary, now I have another one
eg I estimate
> try.op <- polr(
> as.ordered(sod.sit.ec.fam) ~
> log(y) +
> log(1 + nfiglimin) +
> log(1 + nfiglimagg) +
> log(ncomp - nfiglitot) +
> eta +
> I(eta^2) +
>
2006 Aug 15
1
help: cannot allocate vector of length 828310236
Hi all,
I was trying a probit regression using polr() and got this message,
Error in model.matrix.default(Terms, m, contrasts) :
cannot allocate vector of length 828310236
The data is about 20M (a few days ago I asked a question about large file,
thank you for responses, then I use MS Access to select those columns I
would use).
R is 2.3.1, Windows XP, 512M Ram.
I am going to read
2010 Dec 22
1
tests on polr object
Using ordered probit model, I get errors from dwt and bptest.
dwt:
Error in durbinWatsonTest.default(...) : requires vector of residuals
bptest:
Error in storage.mode(y) <- "double" :
invalid to change the storage mode of a factor
I imagine I have to restate as an individual probit model for each category,
but is there an easier way?
thanks,
bp
[[alternative HTML version
2000 Oct 24
2
multinominal probit & logit
Dear everybody!
Are there algorithms for multinominal logit/probit available for R? Is it my
fault that I cannot find these in CRAN? Has somebody programmed these?
with best wishes
Ott Toomet
Ott.Toomet at mail.ee
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info",