similar to: logistic regression 3D-plot

Hello, I have a problem concerning logistic regressions. When I add a quadratic term to my linear model, I cannot draw the line through my scatterplot anymore, which is no problem without the quadratic term. In this example my binary response variable is "incidence", the explanatory variable is "sun": > model0<-glm(incidence~1,binomial) >

Hello! I have created a bplot-figure using this code: *file <- "2dcali_red.ttt" ux<-as.matrix(read.table(file, dec = ",")) mode(ux)<-'numeric' vel<-ux[,1] ang<-ux[,2] x<-ux[,3] y<-ux[,4] dat<- data.frame(ang=ang, x=x,y=y) require(rms) ddist2 <- datadist(dat) options(datadist="ddist2") fitn <- lrm(ang ~ rcs(x,4) +

Dear R- Community, to learn i reanalysed some data provided and analysed by Zuur et. al. in their book "Mixed effect models and Extensions in Ecology with R". When i run the last command i get a warning message i dont understand. Loyn<- read.table(file = "loyn.txt",header = TRUE) Loyn$L.AREA<- log10(Loyn$AREA) fGRAZE <-factor(Loyn$GRAZE) M0<- lm(ABUND~ L.AREA

Hello Dear, I have a function, like z=f(x,y), and try a surface plot with this function. But, on the reference of "wireframe" requires data option, so I generated x and y, and computed z with them. But, still I have a problem to draw a surface plot. The code and errors are ################################################## mle_beta0=64.43707; mle_beta1=-24365.16; # generating for

Hi, trying this again - any help would be great! Basically, I have two matrices of equal dimension, one of them should produce something similar to a heatmap.. The 2nd matrix should be the "heights" for each value on the heatmap - producing a sort of 3d surface plot. Viewing this seems like a problem too, as I need to scale the x axis so that ~140 labels are visible. I've tried

Dear all, Thank you for any help. I have a data.frame and would like to plot it in 3D. I have tried wireframe() and cloud(), I got scatterplot3d(xs) Error: could not find function "scatterplot3d" > wireframe(xs) Error in wireframe(xs) : no applicable method for "wireframe" > persp(x=x, y=y, z=xs) Error in persp.default(x = x, y = y, z = xs) : (list)

It is a trivial matter to adjust the linewidth of the axes and plot line in a 2D plot: x <- 1:10; y <- 1:10; par(lwd=2); # axis (and others) linewidth plot(x,y,"l",lwd=6); # plot linewidth Is there anything comparable to this for persp or wireframe? I've tried a number of ways, and the framebox/axes seem to be linked to the linewidth used in generating the

Please forgive a rather na??ve question... Could someone please give a quick explanation for the differences in conf intervals achieved via confint.glm (based on profile liklihoods) and the intervals achieved using the Design library. For example, the intervals in the following two outputs are different. library(Design) x = rnorm(100) y = gl(2,50) d = data.frame(x = x, y = y) dd = datadist(d);

Hi, I'm trying to plot an equation in two variables to get a feel for sensitivity to its parameters. I've run expand.grid to get made-up vectors of the combinations of the two independent variables, and am trying to plot the output of the dependent, M, against both the dependent in a 3d space. t <- 0:100 DBH <- 10:100 TxDBH<-expand.grid(t,DBH) time<-TxDBH[,1]

Hi everyone, I'm doing a logistic regression with an ordinal variable. I'd like to set the contrasts on the ordinal variable. However, when I set the contrasts, they work for ordinary linear regression (lm), but not logistic regression (lrm): ddist = datadist(bin.time, exp.loc) options(datadist='ddist') contrasts(exp.loc) = contr.treatment(3, base = 3, contrasts = TRUE) lrm.loc =

How do I rotate 3D plots/surfaces generated by either cloud or wireframe? Thanks. Marc Marc W. Zodet, MS Health Statistician Agency for Healthcare Research and Quality [[alternative HTML version deleted]]

I was unsure of what suitable goodness-of-fit tests existed in R for logistic regression. After searching the R-help archive I found that using the Design models and resid, could be used to calculate this as follows: d <- datadist(mydataframe) options(datadist = 'd') fit <- lrm(response ~ predictor1 + predictor2..., data=mydataframe, x =T, y=T) resid(fit, 'gof'). I set up a

Hi, I am trying to run a simple logistic regression using lrm() to calculate a odds ratio. I found a confusing output when I use summary() on the fit object which gave some OR that is totally different from simply taking exp(coefficient), see below: > dat<-read.table("dat.txt",sep='\t',header=T,row.names=NULL) > d<-datadist(dat) > options(datadist='d')

Dear expeRts, I would like to add two points to a wireframe plot. The points have (x,y,z) coordinates where z is determined to be on the wireframe [same z-value]. Now something strange happens. One point is perfectly plotted, the other isn't shown at all. It only appears if I move it upwards in z-direction by adding a positive number. So somehow it disappears in the wireframe-surface

Hi, I have the following problem. Calling wireframe within a loop results into an empty window(s) #generate some data temp = expand.grid(A = 1:3,B = 1:3) temp = cbind(temp, y1 = rnorm(9)) temp = cbind(temp, y2 = runif(9)) #plot y1 and y2 in two different windows for(i in 1:2) { wireframe(y1 ~ A*B, temp, shade =T) windows() wireframe(y2 ~ A*B, temp, shade =T) } #However, calling it twice

Hi, I followed the installation instruction step by step with the wine 1.3.13 package. I stick at point 3.1 user guide: using winecfg. If I input /usr/local/bin/winecfg or only winecfg in the wine-1.3.13 directory, the error message -bash: winecfg: command not found appears. There is a directory winecfg under wine-1.3.13/programs/ but no executable files. The wine package is installed under

Hallo! ? I just want to make a 3D plot of a surface of a cone and want to plot some single points around. ? I tried wireframe but cannot find how to plot single points ? I tried scatterplot3d but there the surface is not simple to plot. And: How can I rotate the point of view by the z-axis ? I tried persp3d but how can I add some single points? ? Example: ? library(lattice) library(scatterplot3d)

Per the message from Alexander Nervedi, 29 April 2006: > I have to be making a riddiculously silly ommission. > when I run the fillowing i get the cloud plot ok. But I cant figure > out what I am missing out when I call wireframe. > Any help would be appreciated. > x<-runif(100) > y<-rnorm(100) > z<-runif(100) > temp <-data.frame(x,y,z) >

Hello R-classicist won't like the following question due to there's no minimal reproducible example and the question was posted already. Anyway I'm not able understand this rather complicated version of "add points to wireframe". what I have is matrix 3x2000 with a dem<-(x,y,z, coordinates) (not provided). with follwing code I want to produce the wireframe and add some

Hello, I try to increase the distance between tick labels and ticks in a lattice wireframe plot. Here's a minimal example: ## Minimal example x <- y <- z <- c(1,2,3) df <- data.frame(x, y, z) wireframe(z ~ x*y, df, scales = list(arrows = FALSE, col = "black", font = 1, tck=0.6)) I tried the axis.components option