[R] Comparing density plots using same axes or same axes scale

2012-09-26 Thread Meredith Ballard LaBeau
Good Evening-
  I have a set of nine scenarios I want to plot to see how the distribution
is changing, if one tail is getting larger in certain scenario, currently I
am using this code


colnames<-dimnames(sag_pdfs)[[2]]

par(mfrow=c(3,3))

for(i in 1:9) {

d<-density(sag[,i])

plot(d,type="n", main=colnames[i])

polygon(d,col="red",border="grey")}

where sag is a 7305x9 double matrix and 9 different scenarios. I want to be
able to compare the distribution using the same axes scale.
Can anyone help?

Thanks
Meredith LaBeau

-- 
Doctoral Candidate
Department of Civil and Environmental Engineering
Michigan Technological University

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Merging multiple columns into one column

2012-09-28 Thread Meredith Ballard LaBeau
Good Evening-
 I have a dataframe that has 10 columns that has a header and 7306 rows in
each column, I want to combine these columns into one. I utilized the stack
function but it only returned 3/4 of the data...my code is:
where nfcuy_bw is the dataframe with 7305 obs. and 10 variables
Once I apply this code I only receive a data frame with 58440 obs. of 2
variables, of which there should be 73,050 obs. of 2 variables, just
wondering what is happening here?

 View(nfcuy_bw)

attach(nfcuy_bw)

cuyahoga_nf<-data.frame(s5,s10,s25,s27,s33,s41,s51,his_c)

cuy_nf<-stack(cuyahoga_nf)

Thanks
Meredith

-- 
Doctoral Candidate
Department of Civil and Environmental Engineering
Michigan Technological University

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Removing lower whisker in boxplot to see the effects of the high values

2012-09-29 Thread Meredith Ballard LaBeau
Good Afternoon-
  I was wanting to alter the boxplot to remove the lower whisker, both the
whisker line and staple just on the lower end. Is there a way to do this?
As my code is currently:
boxplot(log_loads~ind,data=nfmaum, horizontal=TRUE, notch=T, outline=FALSE,
whisker=0, main="Maumee River Near Future Climate Scenarios", ylab="Log
Load",xlab="Climate Scenarios")

I just want to better see the medians and high end tail.

Thanks
Meredith

-- 
Doctoral Candidate
Department of Civil and Environmental Engineering
Michigan Technological University

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.