search for: ltheta

...post the package somewhere -- with the package installed, it's only a few steps to reproduce the example. I end up, deep in the process of trying to compute a likelihood profile, with the following situation: I want to evaluate "call": call mle2(minuslogl = function (lmu = NULL, ltheta = NULL) { if (!is.null(parameters)) { pars <- unlist(as.list(match.call())[-1]) for (i in seq(along = parameters)) { assign(vars[i], mmats[[i]] %*% pars[vpos[[i]]]) } } arglist1 <- lapply(arglist1, eval, envir = data, enclos = sys.frame(sys.nfra...

...a 3d plot. It defaults to the [0,1] interval and when I try to change it I get errors. Example: This works: ------------ D <- c(1,2,3,4,5,6,7,8,9,10) M <- c(11,12,13,14,15,16,17,18,19,20) DM <- cbind(D,M) persp(DM, theta = 40, phi = 30, expand = 0.5, col = "lightblue", ltheta = 120, shade = 0.75, ticktype = "detailed", xlab = "X", ylab = "Y", zlab = "Z") --------------------- But I want the axis to count 1 by ones. So I try: ----------------- D <- c(1,2,3,4,5,6,7,8,9,10) M <- c(11,12,13,14,15,16,17,18,19,20) DM <...

...(m1,m2), dim=c(6,5,2)) ## colors colo <- c("red", "blue") ## axis x <- 1:6 y <- 1:5 z <- mm[,,1] z1 <- mm[,,2] ## surface with heights and colors ## related to the first matrix (no good) persp(x, y, z, theta = 30, phi = 30, expand = 0.5, col = colo, ltheta = 120, ticktype = "detailed", xlab = "X", ylab = "Y", zlab = "values" ) ## surface with heights and colors ## related to the second matrix (no good as well) persp(x, y, z1, theta = 30, phi = 30, expand = 0.5, col = colo, ltheta = 120, tickty...

...lcano) z <- 2 * volcano # Exaggerate the relief x <- 10 * (1:nrow(z)) # 10 meter spacing (S to N) y <- 10 * (1:ncol(z)) # 10 meter spacing (E to W) ## Don't draw the grid lines : border = NA persp(x, y, z, theta = 135, phi = 30, col = "green3", scale = FALSE, ltheta = -120, shade = 0.75, border = NA, box = FALSE) and now I would like to include another surface. I was hoping for a possibility like add=TRUE such as in contour: persp(x, y, z+10, theta = 135, phi = 30, col = "red", add=TRUE) Can someone point out to me how it can be accomplished (mayb...

...tor(y)); # Data z <- f(x,y); fit <- mars(X,as.vector(z),nk=200,penalty=2,thresh=1e-3,degree=2); # Plotting par(mfrow=c(1,2),pty="s") lims <- c(min(c(min(z),min(fit$fitted))),max(c(max(z),max(fit$fitted)))) persp(z=z,ticktype='detailed',col='lightblue',shade=.75,ltheta=50, xlab='x',ylab='y',zlab='z',main='true',phi=25,theta=55,zlim=lims) persp(z=matrix(fit$fitted.values,nrow=nrow(x),byrow=F),ticktype='detailed', col='lightblue', xlab='x',ylab='y',zlab='z',shade=.75,ltheta=50,mai...

...<- x f <- function(x,y) { r <- sqrt(x^2+y^2); 10 * sin(r)/r } z <- outer(x, y, f) z[is.na(z)] <- 1 # 'bg' works but 'tck' is not showing any effect... par(bg="gray90", tck=0.01) persp(x, y, z, theta = 30, phi = 30, expand = 0.5, col = "lightblue", ltheta = 120, shade = 0.75, ticktype = "detailed", xlab = "X", ylab = "Y", zlab = "Sinc( r )" ) -> res Any help would be highly appreciated! Thanks, -- Kim. [[alternative HTML version deleted]]

...t;- (x1 + x2 < 1) term1*term2*term3 } z <- outer(x1, x2, f) z[z<=0] <- NA persp(x1, x2, z, main = "Dirichlet Distribution", col = "lightblue", theta = 50, phi = 20, r = 50, d = 0.1, expand = 0.5, ltheta = 90, lphi = 180, shade = 0.75, ticktype = "detailed", nticks = 5) #<end code> It works fine (I guess), except for a1=a2=a3=1. In that case I get the error: Error in persp.default... invalid 'z' limits. The z matrix has only elements 2 and NA....

...<-data.normal[!is.na(data.normal[,2]),] x<-new.data.normal[,2] y<-new.data.normal[,3] op <- KernSur(x,y, xgridsize=50, ygridsize=50, correlation=0.4968023, xbandwidth=1, ybandwidth=1) #3D density plot persp(op$xvals, op$yvals, op$zden, theta=30,phi=10,expand=0.5,ltheta=120, xlab="TECH3661.A",ylab="TECH3661.B",zlab="Prob",col="pink", , main="3D DENSITY PLOT-TECH3661 ", sub=" TECH3661.A AND TECH3661.B", box = T, axes = TRUE,ticktype = "detailed", ) #countour plot image(op$x...

...gth(var1)] <- max(x[,,2],na.rm=TRUE) nbLines<-min(100,dim(x)[1]) if(missing(color)){color<-2:(nbLines+1)}else{} repeat{ res <- persp(x=time, y=var1, z=var2, theta = angle, phi = 10, expand = 0.5, col = "lightblue", ltheta = 120, shade = 0.75, ticktype = "detailed", xlab = "time", ylab = "var1", zlab = "var2") angle <- angle+1 for(i in 1:nbLines){ yy=x[i,,1] zz=x[i,,2] lines (trans3d(time, yy, zz, pmat...

Hi All, I am currently preparing some form of slideshow introducing R and its capabilities for some colleagues. The thing will be about 30 mins, and I'd like to have some "pretty pictures" and some "amazing facts" (I'm trying to sell, obviously :)). Can I ask if it's possible to easily retrieve a gross figure of the number of functions in R considering the

Dear all; Does anyone know how to add grid lines to a persp plot? I've tried using lines(trans3d..) but the lines of course are superimposed into the actual 3d surface and what I need is something like the plot shown in the following link: http://thermal.gg.utah.edu/tutorials/matlab/matlab_tutorial.html I'll appreciate any ideas Thanks PM

...s and the z-axis are the yields. The following code works fine: persp(y, x, as.matrix(ypper), xlab="Last 50 periods", ylab="Maturities (months)",zlab="Yields", zlim=c(3,6),theta=130, phi=15, col="Seagreen", box=T, ticktype="detailed", shade=0.3, ltheta=270, lphi=30, r=sqrt(3), expand=.75, main="Historical yield curves") Now, I want to replace the x- and y-axis with character strings, such that: temp<-seq(ISOdate(1999,1,1), by="week", length=(nrow(yall))) temp.ymd <- substr(as.character(temp),1,10)[(lastobs-index[length...

...equences of numbers with ordered lists of labels. if I pass levels(Columns) as the y-axis I get : non-numeric argument to binary operator - but I can't work out how to access the labels in any other way. This is the command I am using : persp(Rows,Columns,FRm,xlim=c(0,500),theta=28,phi=40,ltheta=20,lphi=135,col='gray',shade=0.7,r=1,box=TRUE,ticktype='detailed',d=1,expand=0.7,xlab='Frequency(Hz)',ylab="Source's Distance from Floor",zlab='Floor Related Transfer Function(dB)') with yaxis : > Columns [1] 1 2 3 4 5 6 7 8 x axis : > Rows...

Hi .. I would like to have math symbols in perspective plots i tried : persp(x,y,z,xlab=expression(phi)) but it plots it as phi. Thanks. Harsh __________________________________________________ [[alternative HTML version deleted]]

...rm = TRUE)+1), main = bquote(paste("Dirichlet eloszlás, ", alpha ,"=(",.(a1),",",.(a2),",",.(a3),")")), col = "lightblue", theta = 50, phi = 20, r = 50, d = 0.1, expand = 0.5, ltheta = 90, lphi = 180, shade = 0.75, ticktype = "detailed", nticks = 5) } plot.dirichlet<-function(){ refresh.code<-function(...){ a1<-slider(no=1); a2<-slider(no=2); a3<-slider(no=3) type= slider(obj.name="type") dirfelu...

...ed") HTMLon() showgraph() </pre> <p>Here is another graphic: </p> <pre dojoType="Rpad"> data(volcano) z <- 2 * volcano; x <- 10 * (1:nrow(z)); y <- 10 * (1:ncol(z)) persp(x, y, z, theta = 135, phi = 30, col = "green3", scale = FALSE, ltheta = -120, shade = 0.75, border = NA, box = FALSE) HTMLon() showgraph() </pre> </td></tr></table> <br><br><br>_{by Tom Short, tshort at epri.com, Copyright 2005 EPRI, license: GNU GPL v2 or greater} </body> </html>...

Deall all, I'm just learning R, but unfortunately I need to urgently do a rather more complex task so I need some help. I have just learnt the very basics a few days ago and am not ready yet to deal with panels and kernel densities, so a soft guidance would be most appreciated. I have a (very) large panel data set (400,000 individuals x50 time periods) and need to display the evolution of

Hello, I requested help a couple of weeks ago creating a dipole field in R but receieved no responses. Eventually I opted to create a 3d sinusoidal plot and concatenate this with its inverse as a means for a "next best" situation. It seems that this isn't sufficient for my needs and I'm really after creating a continuous 3d gaussian mesh with a "positive" and