similar to: Seasonal correlations

I'm new to R and have searched for help and consulted the the pdf manuals, but I can't seem to figure out how to plot more than one series on the same graph. I've tried using multiple par(new=TRUE) statements such as >plot(series1, ci.type="line", col="red", lwd=2, ci.lty=0, ci.col="red") >par(new=TRUE) >plot(series2, ci.type="line",

Hi all- I am stumped. The code in A. returns errors at lines 14 and 15 and fails to load series1 and series2. However, in B., if temp1 and temp2 are called again (which returns a "Table exists" error; see lines 14-17 in B.) series1 and series2 load correctly. Any ideas? Also-I am open to any suggestions to improve the code as I am a horrific programmer. Thanks A. 1 >

Hello everybody, I have two series of intervals, and I'd like to output the shared regions. For example: series1<-cbind(Start=c(10,21,40,300),End=c(20,26,70,350)) series2<-cbind(Start=c(25,60,210,500),End=c(40,100,400,1000)) > series1 Start End [1,] 10 20 [2,] 21 26 [3,] 40 70 [4,] 300 350 > series2 Start End [1,] 25 40 [2,] 60 100 [3,]

Howdy. Been running into a bit of trouble with plotting. Seems that axes=F is not "working". Whenever I plot (either a dataframe or xts/zoo series) and I set axes=F along with xlab/ylab="" I still get the default axes printed in my chart. Consider this: #Create some sample data, both 50 units of blah series2 = c(1:50) series1 = rep(25:74) testdf1 = as.data.frame(series1)

Hello users, I've noticed a problem when creating a stacked column plot when a column contains a negative and a positive value. e.g. series1<-c(-1,-2, 3, 4, 5) series2<-c( 5, -4,-3,-2, 1) data<-rbind(series1,series2) barplot(as.matrix(data), beside=FALSE) In these cases (i.e. first, third and fifth columns) the plotting is not handled correctly. Compare this output with that

Dear List, I would be pleased if someone can help me with the following issue: I'm about to plot two time series in one plot via ts.plot which looks like: ts.plot(series1, series2, main=main, xlab=xlab, ylab=ylab, col=c("green", "red", "blue"), lwd=2) The problem is, that R automatically sets the x axes labels in 5-year-intervalls. Every

Dear R heleprs   Suppose I have following data   Scenarios combination_names     series1    series2 Sc1 MAT2 GAU1 7.26554 8.409778 Sc2 MAT2 GAU2 7.438128 8.130275 Sc3 MAT3 GAU1 8.058422 8.06457 Sc4 MAT1 GAU2 8.179855 8.022071 Sc5 MAT3 GAU2 8.184033 8.191831 Sc6 MAT3 GAU2 7.50312 8.232425 Sc7 MAT1 GAU2 7.603291 8.200993 Sc8 MAT1 GAU1 8.221755 8.380097 Sc9 MAT3 GAU2 7.904908

uugh : i promise that this will be my last question of the day. i hate to constantly bother this group but it takes me time to get familar with all of these functions, tricks and and manipulations. i appreciate everyone's patience. i used too splus a lot but i've gotten rusty. i have a matrix of say 200 rows and 600 columns. i have a function "getprofit" that takes two series

I've created a dataframe containing multiple ACF lists through the command rev.acf<-apply(rev.matrix, 2, acf, na.action=na.contiguous, lag.max=12, plot=FALSE) where rev.matrix is an n by t matrix containing n time series in columns. I'd now like to pull out only the ACF information and store it in a seperate n by 12 matrix. So far the only way I can work out to access this is

Hi you guys: I have been wondering if there is any way to change the labeling in plot.ts( ), for example , if I plot two sequences, i always got y labels as "series1", "series2", I tried to use ylab=c((expression(mu_1)),(expression(mu_2))), but it doesn't work. thanks in advance for any help best.

I having problem with dendrogram leaf names when I read a tab delimited file into dataframe; I have a text file, tab delimited, using read.table into a data frame as follows: > test1<-read.table("c:\\R\\data\\Tremont4.txt", header=TRUE, sep="\t") When I do this the "test1" data frame is picking up my first column names as part of the data and not the case

Sadly you failed to set your email program to send plain text and the data is corrupted at my end. I also think you need to reduce the size of the data set... the intent here is to increase your understanding, not debug your particular analysis. I will say that I am having a very challenging time understanding what you are trying to accomplish though. What are the equations that you think need

I can't get my simultaneous equations to work using system fit. Please help. #Reproducible script Empdata<- read.csv("/Users/ngwinuiazenui/Documents/UPLOADemp.csv") View(Empdata) str(Empdata) Empdata$gnipc<-as.numeric(Empdata$gnipc) install.packages("systemfit") library("systemfit") pdata <- plm.data(Empdata,

OK, Let's try this again! Here is the reproducible script; it is long because I had to copy the panel dataset here. My question is related to systemfit; I don't know how to get the result for the entire panel. #Reproducible script Empdata<- read.csv("/Users/ngwinuiazenui/Documents/UPLOADemp.csv") View(Empdata) install.packages("systemfit")

Hi r-users,   I would like to extract the data that match.  Attached is my data: I'm interested in matchind the value in column 'intg' with value in column 'rand_no' > cbind(z=z,intg=dd,rand_no = rr)             z  intg rand_no    [1,]  0.00 0.000   0.001    [2,]  0.01 0.000   0.002    [3,]  0.02 0.000   0.002    [4,]  0.03 0.000   0.003    [5,]  0.04 0.000   0.003    [6,] 

I have a table like the following. I want to fit Cm to Vm like this: Cm ~ Cl+Q1*b1*38.67*exp(-b1*(Vm-Vp1)*0.03867)/(1+exp(-b1*(Vm-Vp1)*0.03867))^2+Q2*b2*38.67*exp(-b2*(Vm-Vp2)*0.03867)/(1+exp(-b2*(Vm-Vp2)*0.03867))^2 I use nls, with start=list(Q1=2e-3, b1=1, Vp1=-25, Q2=3e-3, b2=1, Vp2=200). But I always get 'singlular gradient' error like this. But in SigmaPlot I can get the result. How

Recently there was a question on using lme with large data sets. As an experiment I fit a linear mixed-effects model to a data set with about 350,000 observations on 6 predictors, a numerical response, and a single grouping factor. The timings shown below were on a 1.2 GHz Athlon with 1 GB of PC133 memory and 2 GB of swap. The operating system is Debian 3.0 GNU/Linux. The kernel is 2.4.14.

I am trying to learn how to format Ascii data files for scan or read into R. Precisely for a quickie project, I found some code (at end of this email) to do exactly what I need: To cluster and graph a dendrogram from package (stats). I am stuck on how to format a text file to run the script. I looked at the dataset USArrests (which would be replaced by my data and labels) using UltraEdit. That

How do I plot two time series plots on the same chart? Thanks, Jason [[alternative HTML version deleted]]

Hi R user I need to read in some values from a computer program output. I can't change the output format because the developer of the program doesn't allow to change the format of output. There are two formats. First one looks like this if I have 10 variables, ------------------------------------------------------------------------------------------------------ [ 1]