search for: logit

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,...

The following is a simple example with a poor solution that shows my difficulties with scope. The function /logit.test /has 3 arguments: /model.start,/ an initial model; /model.finish/, an all-inclusive model, /my.data/, a dataset, in this case trivial. There are 2 function calls in l/ogit.test,/ first to /glm/ to get an initial fit (local variable /logit/) using /model.start;/ then a call to /stepAIC/ u...

...;t suggest to use the quasibinomial family. This # works out, but when applied to the true data, the # variance function does not seams to be # appropriate. # # I couldn't see in the # theory why this does not work. # Is this a bug, or are there theoretical reasons? # One problem might be, that logit(0)=-Inf and logit(1)=Inf. # But I can't see how this disturbes the calculation of quasi-Likelihood. # # Thank you very much, # best, # # Bj??rn set.seed(0) y <- sample(c(0,1), size=100, replace=T) # the following models work: glm(y ~ 1) glm(y ~ 1, family=binomial(link=logit)) glm(y ~ 1, fa...

Hello dear r-gurus! I have a question about the logit-model. I think I have misunderstood something and I'm trying to find a bug from my code or even better from my head. Any help is appreciated. The question is shortly: why I'm not having same coefficients from the logit-regression when using a link-function and an explicite transformation o...

Hell all. I¡¯m trying to estimate mixed logit model using MSLE. In order to see that mixed logit model works better than simple logit model ( the logit model with fixed coefficient) I simulated a dataset with random coefficients and tried to fit the data with both mixed logit and simple logit model. Because my mixed logit model contai...

Hello, Could anybody provide me with codes (procedure) how to obtain Conditional Logit (McFadden) and Mixed Logit (say, assuming normal distribution) estimates in R? Thanks, David U. -- View this message in context: http://www.nabble.com/Conditional-Logit-and-Mixed-Logit-tf4489238.html#a12802959 Sent from the R help mailing list archive at Nabble.com.

...core matching. For each treated unit, I want to find all control units whose propensity scores lie within a certain distance from the treated unit. The sample code is as follows: > library(Matching) > x <- rnorm(100000) > y <- rnorm(100000) > z <- rbinom(100000,1,0.002) > logit.reg <- glm(z~x+y,family=binomial(link='logit')) > match <- > Match(Y=NULL,Tr=z,X=logit.reg$fitted,version='fast',ties=TRUE,M=1,distance.tolerance=1e-5) According to the function definition (http://sekhon.berkeley.edu/matching/Match.html): "distance.tolerance: Thi...

I'm running a multinomial logit in R using the Zelig packages. According to str(trade962a), my dependent variable is a factor with three levels. When I run the multinomial logit I get an error message. However, when I run 'model=logit' it works fine. any ideas on whats wrong? ## MULTINOMIAL LOGIT anes96two <- zelig(tr...

...l( p_{ij}, N_{ij} ) I've generated data using this model, and I'd like to fit it. My data is a data frame with 3 columns, "response" (0/1), "subject" (a factor), and "condition" (another factor). Here is my link definition: ############################# halflogit=function(){ half.logit=function(mu) qlogis(2*mu-1) half.logit.inv=function(eta) .5*plogis(eta)+.5 half.logit.deriv=function(eta) .5*(exp(eta/2)+exp(-eta/2))^-2 half.logit.inv.indicator=function(eta) TRUE half.logit.indicator=function(mu) mu>.5 & mu<1 link <- "half.logit" stru...

Hi I am carrying out some logit regressions and want to (a) make sure I'm taking the right approach and (b) work out how to carry out some additional analysis. So, to carry out a logit regression where the dependent variable is a factor db, I use something like: res1_l <- glm(formula = db ~ y1 + + y5, family = binomial...

Hi all, Thank you for taking the time to read my message. I'm trying to make a figure that plots p-values by a range of different adjustment values. (Using the **logit** function in package **car**) My Statistical analyses were conducted on probability estimates ranging from 0% to 100%. As it's not ideal to run linear models on percentages that are bounded between 0 and 1, these estimates were logit transformed. However, this introduces a researcher degree...

Hi I am trying to make glm() work to analyze a toy logit system. I have a dataframe with x and y independent variables. I have L=1+x-y (ie coefficients 1,1,-1) then if I have a logit relation with L=log(p/(1-p)), p=1/(1+exp(L)). If I interpret "p" as the probability of success in a Bernouilli trial, and I can observe the result (...

Good morning, I'd like to estimate a simple multinomial logit model in R (not a McFadden conditional logit). For instance, I'd like to estimate the probability of someone having one of eight titles in a company with the independent variables being the company characteristics. A binary logit is well documented. What about the multinomial? Thanks, Walt...

Hi, I'm new to the glm and logit world... and I'm reading some lecture notes and examples. I would like to try and generate the same result in R.. but I don't seem to be able to find the proper way to specify the formula.... let's say i have Desire Using Drugs Not Using Drugs Yes 219...

Dear R users this may be a simple question - but i would appreciate any thoughts does anyone know how you would get one lower and one upper confidence interval for a set of data that consists of proportions. i.e. taking a usual confidence interval for normal data would result in the lower confidence interval being negative - which is not possible given the data (which is constrained between

...mily=binomial like so: >yy succ fail [1,] 76 23 [2,] 32 67 [3,] 56 43 ... [24,] 81 18 >xx c1219 c643 X1 0.04545455 0.64274145 X2 0.17723669 0.90392792 ... X24 0.80629054 0.12239320 >test.gl1ce <- gl1ce(yy ~ xx, family = binomial(link=logit), bound = 0.5 ) or >omit <- c(2,3) >test.gl1ce <- gl1ce(yy[-omit,] ~ xx[-omit,], family = binomial(link=logit), bound = 1 ) this seems to work fine and as i change the shrinkage parameter everything behaves as expected. if i try to get the fitted values (y-hat) using predict i ha...

I have a logit model with about 10 predictors and I am trying to plot the probability curve for the model. Y=1 = 1 / 1+e^-z where z=B0 + B1X1 + ... + BnXi If the model had only one predictor, I know to do something like below. mod1 = glm(factor(won) ~ as.numeric(bid), data=mydat, family=binomial(link="l...

...e cron job is hitting some sort of resource limit?? If I run the perl script from the command line it appears to dsync all the accounts OK.?? The relevant part of my Perl code is below.?? Any ideas? Thanks, Murray foreach (@userlist) { ??? my $email = $_; ??? chomp ($email); # get rid of CR ??? logit("Started sync of mailbox $email") if $debug; ??? my $dsync_command="dsync -v -D -u $email mirror ssh -C root@$mailstore dsync -u $email"; ??? logit("dsync_command=$dsync_command"); ??? # Run dsync of mailbox ??? logit("Started secondary sync of $mailbox to $mails...

Hello, I'd like to fit a logit function to my data. The data is distributed like a logit (like in this plot on wikipedia http://en.wikipedia.org/wiki/File:Logit.png) but the values on the x-axis are not between 0 and 1. I don't think using a glm is the solution because I simply want to infer the parameters of the logit fun...

...t;not realized", "realized") ~ v1 + v2 , family="binomial") However, sometimes I don't have the number of cases (realized, not realized), but only the proportion and thus cannot compute the binomial model. I just found out that the package car contains a function "logit" which allows for logit transformation. Would it be possible to transform the proportion data with this function and analyze the transformed data with a glm with family="gaussian"? Thank you very much -- View this message in context: http://r.789695.n4.nabble.com/function-logit-v...