similar to: Package for multivariate binary logistic regression?

I've been using R to do logistic regresssion, and that's working well, but there are two things I haven't figured out how to do. (1) Is there some pre-existing function that will let you compute the odds ratios and confidence intervals for them for a specific fit. I know how to do this manually or even write a function that I can call with the coefficients and se, but

I want to get the logLik to calculate McFadden.R2 ,ML.R2 and Cragg.Uhler.R2, but the value from multinom does not have logLik.So my quetion is : is logLik meaningful to multinomial logistic model from multinom?If it does, how can I get it? Thank you! ps: I konw VGAM has function to get the multinomial logistic model with logLik, but I prefer use the function from "official" R

Hi, We?re looking at options for logging agents into the system programmatically via Perl/PHP and I was wondering if anyone else is doing this and if so, how. We're using AgentCallbackLogin now but would like to set up a web interface instead. I've been looking at Asterisk::Manager and didn't see anything relevant and wanted to ask the group before we dove into the Asterisk source.

> > >From: robert-mcfadden w o2.pl > > >Date: 2008/04/02 Wed AM 09:58:28 CDT > > >To: r-help w r-project.org > > >Subject: [R] Nonlinear equation > > > > hi: you need to give an example and details or > > you won't get much response, if any. Equation e.g. (A, B are known constants): 3log(gamma(x))-log(gamma(x)*gamma(2x))+(x-1)*A+B=0

Hi all I know this must be an easy one so sorry for the trouble. I would like to select a list of variables within a factor The following example is given in help for subset: subset(airquality, Temp > 80 select = c(Ozone, Temp)) So how do I select all temperatures of 90 and 80 ie Temp = c(80,90) I would appreciate your help. Kind regards andy Andrew McFadden MVS BVSc Incursion

Dear all, can somebody please help me how to calculate "McFadden R-square" for a LOGIT model? Corresponding definition can be found here: http://publib.boulder.ibm.com/infocenter/spssstat/v20r0m0/index.jsp?topic=%2Fcom.ibm.spss.statistics.help%2Falg_plum_statistics_rsq_mcfadden.htm Here is my data: Data <- structure(c(1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1,

I all I am just trying to get the idea of a loop if statement for another purpose. I was going to label a new variable based on date ie first week 1, second week 2 etc. My code in bold is wrong but seems logical to me, could someone help. x <- rep(c(2660156,2663703,2658165,2659303,2661531,2660914),c(2,2,2,2,1,1)) y <- rep(c(6476767,6475013,6475487,6479659,6477004,6476388),c(2,2,2,2,1,1))

Hi I have estimated parameters of my data with mlogit and the following commands. I would like to know also the McFadden R^2 and the intercept, could soweone tell me how that can be done? library(RODBC) library(mlogit) library(foreign) z<-odbcConnectExcel("D:\\MALLI11ARVOT.xls") y<-sqlFetch(z,"Taul1")

A couple of questions: How to be sure that gls allowes errors to be correlated and/or have unequal variances? (is this on auto or is there a switch?) How to calculate confidence limits for a linear regresssion? -------------- next part -------------- A non-text attachment was scrubbed... Name: dthompson.vcf Type: text/x-vcard Size: 303 bytes Desc: Card for David Thompson Url :

Hi, I'm a new user of R and I'm trying to make a linear model from this kind of dataset x [1] 16.87 19.93 25.85 20.94 17.06 19.49 19.93 25.45 27.74 20.15 25.81 21.06 17.17 20.03 25.50 27.79 20.44 16.88 19.93 25.79 z<-x-10 y [1] 0.80 1.27 2.22 1.32 0.90 1.18 1.84 2.41 2.97 1.25 2.07 1.41 1.14 1.66 2.59 3.51 1.53 0.81 1.26 2.30 plot(x,y) I want to be able to force the line of

Dear all: Does anyone has harmonic regresssion analysis package written in R (to be used in Windows platform) ? Thanks _______________________________________ YM - Â÷½u°T®§ ´Nºâ§A¨S¦³¤Wºô¡A§AªºªB¤Í¤´¥i¥H¯d¤U°T®§µ¹§A¡A·í§A¤Wºô®É´N¯à¥ß§Y¬Ý¨ì¡A¥ô¦ó»¡¸Ü³£ÉN¨«¥¢¡C http://messenger.yahoo.com.hk [[alternative HTML version deleted]]

Hi R users I apologise in advance for this question as I suspect it is simple and perhaps others have had this problem. I am struggling to sort out how to fix the x axes so that the labels don't overlap. I have put the following example together to show my problem. library(lattice) titre <- as.factor(rep(c(10999,20999,30999,40999,50999,60999,

Hi all The code of trying to write relates to selecting properties (given by x and y co-ordinates) spatially (distance "X" from "infected" properties identified by date) over a certain time period. i.e. what properties are within 3 km from properties infected on "2008-01-01" over the last 14 days. Is any one able to give me some clues on how to write code to solve

I am hoping I can get a little light she on a problem I am having trying to bring up a guest domain. I suspect that the problem is in my supplied file system First the environment. VT enabled dual core EM64T processor, Asus Mb (I can get the exact chipset infor if relevant) 2 Gig of DRAM SATA drives. I have build and booted the domain 0 using 3.0.2-2 version of xen. This seems ok.

Hi I am pretty new to this, and this may be an easy one, but I have not found the answer in the docs. The kernels produced from a build seem to have the same name for both the dom0 and domU instances. I get one kernel named vmlnuz-2.16.6-xen Is this expected? Although the dom0 boots, I ma having a crappy time getting a virtualized OS to boot and am trying to see if this is the

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 Paczkowski

Hi all I am trying to do a simple recode which I am stumbling on. I figure there must be any easy way but haven't come across it. Given data of A","B","C","D","E","A" it would be nice to recode this into say three categories ie A and B becomes "Treat1", C becomes "Treat 2" and E becomes "Treat 3". I tried

Hi R users I am trying to create a spatial join between two datasets. The first data set is large and contains descriptive data including x and y co-ordinates. The second dataset is small and has been selected spatially. The only data contained within the second dataset is the x and y coordinates only i.e. no descriptive data. The aim of a join made between the two datasets is to select

I all I am trying to combine columns from two dataframes to make a completely new dataframe consisting of one column of dates (ie a combination of dates from the two dataframes). >From the following dataframes a b 1 2008-07-27 1 2 2008-10-01 2 3 2008-08-15 3 4 2008-08-14 4 5 2008-08-14 5 6 2008-09-20 6 c d 1 2008-07-27 1 2 2008-10-01 2 I would like to get: z

Hi all I have created a lattice graph but the text appears too big. Just wondering if there was a way I could change the font size. tp1.nDai=xyplot(value ~ yr | Name, data = nDai,type = 'o') update(tp1.nDai, aspect = c(0.8), layout = c(4, 4)) I would appreciate your help. Kind regards andy Andrew McFadden MVS BVSc Incursion Investigator Investigation & Diagnostic Centres -