similar to: Ploting betas and p-values together

[Please use the subject line!] In the help page for summary.lm, the "Value" section says that the returned object has a component called "fstatistic", which has the F-statistic and the associated numerator and denominator degrees of freedom. You can get the p-value by something like: fstat <- summary(speciallinearmodel)$fstatistic pval <- pf(fstat[1], fstat[2],

Hi All, I have checked everything I could find abot graphics, but still cannot solve the problem. Are there any ways to make a graph that plots two lines and two different y-axes, each of them has a scale that is related to the respective line. For example, y1 has a range 1:50 and y1 ranges 0:1. The x-axe is the same for both. Thank you in advance. --- Gregor Gawron

Dear all, I am applying the StructTS function in ts-package. For some time series the program terminates and the following error appears: Error in optim(init[mask], getLike, method = "L-BFGS-B", lower = rep(0, : L-BFGS-B needs finite values of fn Do someone know what do I have to adjust in the original time series to avoid this error? It works fine for some subsets of the

Dear R-users, I am facing a trivial problem when trying to parameterise the start date of a time series object. I am working with monthly data (104) performing n-steps-ahead (6) forecasts and using a fixed window size (36). At the end of calculations I have a list that contains 69 forecasts. I have no problems in fixing the window size by parametrization, e.g. k<- control variable in a for

Hi to all, First of all, I prefer to tell that I am a R-newbie, so I apologize if this is a silly question (I have tried looking in the manuals, but without luck). I have two variables, y and z, that I want to plot against x in the same plot. I have done this before, using points() after plot(). But now the problem is that y and z are in different units of measurement, and their ranges are very

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

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

Hi: I been having a lot of trouble trying to plot multiple time series on the same plot. What I want to do is the following: I have a table with sales per day on different states and what i would like to do is plot time series for all the diferents states in the same plot (divided in small squares, not all the time series together) so I can visualize the data; but I do not know how to do this, im

Hi, I have trouble with the heatmap function (package stats). The row labels are wrongly ordered and don't correspond to the Rowv dendrogram. I know there is a bug with the heatmap fonction. Emmanuel Paradis (http://tolstoy.newcastle.edu.au/R/e2/help/07/05/16227.html )suggested a modification to fix it but in my case the row labels are still wrongly ordered. Heatmaps with 2 phylograms have

I'm trying to understand how to set graphical parameters for trees created with the party package.  For example take the following code:   library(party)     data(airquality)     airq <- subset(airquality, !is.na(Ozone))     airct <- ctree(Ozone ~ ., data = airq,                    controls = ctree_control(maxsurrogate = 3))     plot(airct)   My problem is, I've got a ctree that has

I used in previous versions of R (R 1.1) to plot axes on projections of factorial coordinates using the function plaxes. This doesn't exist any more. Why this change ? No alternative ? --------------- Charles RAUX, Laboratoire d'Economie des Transports CNRS-Universit? Lumi?re Lyon 2-ENTPE email : charles.raux at let.ish-lyon.cnrs.fr http://www.ish-lyon.cnrs.fr/let

Hi, It's my first post on this group. I've just started learning & using R and I like it ;-) I have I think simple question. I'm trying to plot a histogram for my data set. My data set is defined as follows: Class N 12.5 3 17.5 10 22.5 12 27.5 8 32.5 7 37.5 3 42.5 4 47.5 2 Class means middle of set of my ranges I define. N column stores number of measurements counted to

Dear all, I wanted to plot the mean and standard error on the regression equation (instead of individual data points) in the following code, but I could not find the right code in the help files. Could someone please show how to do this. Thank you very much. temp <- c(16,16,16,16,16, 20,20,20,20,20, 24,24,24,24,24, 28,28,28,28,28, 32,32,32,32,32) dev1hr <- c(36.2, 34, 32.2, 36.4, 36,