similar to: lattice package - wireframe plot : adding more than one surface and addiding a curve overlaid on the plot

Dear All, I know this may be a trivial question. In the past I have used glm to make logistic regressions on data. The output creates an object with the results of the logistic regression. This object can then be used to make predictions. Great. I have a different problem. I need to make predictions from a logistic regression taken from a paper. Thus I need to (by hand) enter the reported odds

Dear users, I can access an database oracle database using sql developer. This was done by importing an xml file that contains the login details - username, password and specifies that it uses the KERBEROS_AUTHENTICATION. I'm trying to connect R - so that it can access this database - do sql queries and convert the resulting tables into dataframes. I am a novice in SQL and database access

All, How does one extract the level-1 variance from a model fit via lmer()? In the code below the level-2 variance component may be obtained via subscripting, but what about the level-1 variance, viz., the 3.215072 term? (actually this term squared) Didn't see anything in the archives on this. Cheers, David > fm <- lmer( dv ~ time.num*drug + (1 | Patient.new), data=dat.new )

Is anyone using the setup described at the Rails blog? http://weblog.rubyonrails.org/2006/7/3/pound-makes-lighty-and-mongrel- play-nice I''m using it in production, and the only problem I have is the same as people have described in the comments over there, namely that Rails thinks all requests come from localhost. See http://isabont.com/requestinfo for the headers. Any

Buenas tengo una matriz de contingencia de 5x4 donde en la última fila y última columna tengo las frecuencias marginales, en función de estas quiero sacar las esperadas, pero cuando hago por ejemplo el producto,, Esperadas<-Estrategia["Suma",]*Estrategia[,"Suma_e"] E1 E2 E3 E4 Suma 5550 14160 52650 89600 20720 Warning message: In Estrategia["Suma", ]

Let's see I've got a network about 10 PCs running XP and 1 running Linux. I had setup Samba and it was running like a charm. XP pcs could access linux box through network neighborhood and access 2 shared dir and edit/change files. Really not sure what has gone wrong. I am using Samba 2.2.8a I run both smbd and nmbd as daemons. When I run ps -aux they show up. root 1553 0.0

Rails Core Weekly July 16 - August 6 Rails Core Weekly summarizes the rails-core mailing list, a list dedicated to Ruby on Rails internals and its development. RCW is brought to you by Rodney Ramdas and proofread by the good people of The Caboose (http://blog.caboo.se/) for accuracy. RSS feed and RCW home: http://www.pinupgeek.com/articles/category/rails-core-weekly-news The ruby on rails

I'd like to plot n (say n = 10) semi-transparent lines in such a way that if all n happen to exactly overlay each other, the resulting line is completely opaque. In principle, I believe that this is what setting alpha = 1/n should accomplish. In practice, different values of n produce different results. Consider the following function, which overlays n semi-transparent red lines: mylines

Hello all-- I'm attempting to produce overlaid histograms with partially transparent columns. Whether this display will end up being useful, I can't say. But I do want to get it right. I've already got one solution (shown below), but I tried some other versions and had questions about my results. (Note: I'm using a quartz device, so transparency shows up correctly. You might

I am new to R and am trying to solve the following problem: I have a data file containing tick-by-tick, millisecond level prices for some stocks. I have another file or two containing orders and trades, again, with millisecond time-stamps. Both of these files are irregularly spaced and the time stamps are in an iso format (<date> <time>.<millisecond>) I would like to create a

In the following code, the only difference between the two plots is the order the variables are plotted. In this case, the plot of "cdata.den" in plot #1 is different from its plot in #2. Specifically, "cdata.den" spans the x-axis from -5 to 30 in plot #1 and from 0 to 20 in plot #2. Does anyone understand why these two plots do not yield the same result? #Make density

Hi, I am facing the problem that I want to plot a histogram chart set freq to true and overlay with normal or weibull or exponential distribution curve. The sample code is shown as below: >samp<-c(-8.2262,-8.2262,-8.2262,-8.20209,-8.09294,-8.07321,-8.07321, -8.07321,-8.07175,-8.04948,-8.04948,-8.04948,-8.03848,-8.03848,

I'm trying to make a set of contour plots of bivariate kernel density estimates, showing three such plots overlaid, similar to this plot http://euclid.psych.yorku.ca/SCS/Private/Test/ridge-boot2.pdf except that I would like to have the contours *filled* (using transparent colors). To make this reproducible, I've saved the results of KernSmooth::bkde2D() in the following file:

Dear R-Users, I am learning how to use the package VGAM to do multinomial regressions. I have worked through the example provided by UCLA stats group: http://www.ats.ucla.edu/stat/r/dae/mlogit.htm However - what I would like to learn is how to use the constraint option. So in the example given in the link i would like to introduce a constraint that does not allow individuals of age less than

Please, could you tell me how to enter a mixed or a purely conditional multinomial logit model in NNET. I know how to do a multinomial logit in NNET but I don't know how to do conditional or mixed models using this package. I do know how to do this with VGAM - but would like to compare my results to those obtained by NNET. Thanks. Raffaele.

Dear Help, I am using the library RODBC to export data into an excel file with different sheets. Here is the code: library(RODBC) save2excel <- function(x, t.name,channel) { sqlSave(channel, x, tablename = t.name, rownames = FALSE) } ## Write Parameters to Excel channel <- odbcConnectExcel(paste(data.dir,"Calibrated_Parameters_Full_Set",".xls",sep=""),