similar to: Ploting different values with different colors

library (ellipse) shape1 = c (1, 0, 0,1) dim(shape1) = c(2,2) ellipse (center = c(0,0), shape = shape1, radius = 1) ============================= Error in plot.xy(xy.coords(x, y), type = type, col = col, lty = lty, ...) : plot.new has not been called yet It is really frustrating. Also what do the shape matrix, radius correspond to an ellipse function (x-x0)^2/a + (y-y0)^2/b = 1 ?

Hi, I'm a SAS user trying to convert myself to R but I still have problems with some pretty simple commands. First I wanted to add up a number of red-tailed hawks seen per day (julian day) per year. So I tried: RTyrday <- tapply(RThr,list(year,julian),sum) And then I tried the following regression: mod1 <- glm(RTyrday~julian+year, family=gaussian (link=identity),data=RT) Wich

I can't seem to be able to change the font size in an effect display. I've tried the following: > par(cex.lab=4) > plot(effect ("alti",reg8), ylab="detection probability") and > plot(effect ("alti",reg8), ylab="detection probability", cex=4) but nothing changes. Can anyone help me? thanks. Emilie Berthiaume graduate student

Hi list! I did the following regression: reg1 <- glm(alti~sp + ovent + vivent + nuage, family=gaussian, data=meteo1) I was interested in knowing the effect of the species (sp) in reg1 and so I used the function «effect»: effect.sp <- effect ("sp", reg1, se=TRUE) with this output: sp AK BW NH OS RT SS 2.730101 2.885363 2.753774 2.750311

I'm trying to use multiple plotting colors in my code. My first "ifelse" statement successfully does what I want. However, now I want anything less than -4.5 to be green and the rest black. I want another "col" argument but can only use one. How could I go about getting separate colors for anything above 4.5 and less than -4.5? plot(three, type="h",

Hi, I'm just getting started with R and I'm having problems with some simple operations: I want to get the the sum of the column "SStot" for each year using the function by. The data set is named "SS". I've tried this: by (SS, year, sum(SStot)) and it's not working. Is it because there's a different number of rows for each year? How else can I do

I tried the following ordered logistic regression in R: mod1 <- polr(altitude~sp + wind_dir + wind_speed + hr, data=altioot) But when I asked The summary of my regression I got the folloing error message: > summary (mod1) Re-fitting to get Hessian Error in optim(start, fmin, gmin, method = "BFGS", hessian = Hess, ...) : the initial value of 'vmin' is not

Greetings, I have what I hope is a simple question. I would like to change my contour interval on the vis.gam( plot.type="contour") in the mgcv package. Is this a situation where I need to modify the function or is there a default value I can change? Thanks

Hello, I am facing a problem in drawing heat map using R. I have a 70X3 matrix and I want to draw a heat map with 1 coloumn on X axis another on Y axis and want to show the value of 3rd coloumn using heat colours. Could you please help me with this . Thank you . Regards, Sumit [[alternative HTML version deleted]]

I was wondering if someone had some X11 fonts that look like windows font? or the configuration libes to add into .wine/config to get them. I tryed to add "New Times Roman" = "-adobe-times-" but wine just crashed.

http://r.789695.n4.nabble.com/file/n2260087/%E6%8D%95%E8%8E%B7.png I am going to plot my data set like this, with means and 25% & 75% quentiles. I've tried "boxplot", but the output is not what I want. Should I use other functions? Thanks -- View this message in context: http://r.789695.n4.nabble.com/ploting-dots-with-quentiles-tp2260087p2260087.html Sent from the R help

I'm trying to do something very simple... I wan to plot a histogram where the y axis represent the percentage of the total sample that each bin represents. I know how to plot a histogram with the counts and density... but can't find anything that gives me perenct of sample on the y axis. Any help is appriciated Below is the script I'm working with par(mfrow=c(1,2))

hi: Yesterday I post a message about hoy to plot a time series, but someone told me to post more information about the file so here it is: the file was read using read.table and the name is list. When I use str(list) it tells the following variables: YEAR int: 2003,2003,2003....2004 MONTH int:1,1,1,1,.... DAY:int 1,1,1,... STATE: factor with 51 levels. SALES: int (sales per day) The reason why

Hi, I have spatial data for plants that I would like to plot according to its radius using sysmbols(). Also I would like to color the circles according to an individual's age using heat.colors(17) [ages within the whole data set range from 0:16]. I have not found an easy way to assign ages (=a) a specific heat.color. Any suggestions? Thanks, Patrick x <- c(1,3,4,6,7,3,2,8,3,6,5,5) y

Dear all, I would like to draw a 3-D horizontal cylinder preferably in RGL device (because this gives the look from different angles). Basic idea is from http://www.tau.ac.il/cc/pages/docs/sas8/insight/chap18/sect3.htm. Below is the description exactly what I want to do. Please see at figure 18.6, 1st plot. Here it is an confidence ellipsoid. Suppose now you put another 4 same ellipsoids on top

En innebygd og tegnsett-uspesifisert tekst ble skilt ut... Navn: ikke tilgjengelig Nettadresse: https://stat.ethz.ch/pipermail/r-help/attachments/20080208/5aae5afd/attachment.pl

this is the command i made for a normal distribution, but when i try to plot the histograms, i dont know why the bars don't stick on the line... nsamples<-1000 sampsize<-15 Samples<-matrix(rnorm(nsamples*sampsize,0,1),nrow=nsamples) a<-apply(Samples,1,var) NC14<-a*14 x<-0:40 plot(x,dchisq(x,14),type='h') hist(NC14,freq=F,add=T) -- View this message in context:

Hi, I would like to be able to develop a function to plot an igraph object with lattice (trellis type displays will be usefull for grouping etc). Anyway, I mostly feeble Igraph requires that you convert two columns of data two an igraph object and to be able to plot the graph...I have tried a very, very simplistic (if not naive) approach and surprise, surprise, it didn't work. First of all,

Hello! I have such a problem... Estimated a model based on common data (you can find it in R library), and I wanted to plot the orginal values with the estimated one. Unfortunately I can only see the original values. Below is the code with data library: / library(forecast) data(AirPassengers) AP <- AirPassengers class(AP) start(AP) end(AP) frequency(AP) lgAP <- log(AP) t<-2:length(AP)

Hi, I've written a function: myfun <- function(x, y) { // blah blah } and I want to graph it. My plan is to use persp(), and my question is: how do I create the array of values? One possibility is: x <- seq(0, 10, by=.1) y <- seq(0, 10, by=.1) inputs <- <somehow create an array of x,y pairs> outputs <- apply(A, c(1,2), myfun) The "inputs" array would