First let me say thank you for all the help everyone gave me.
Below is the code I came up with to look at Rainfall pH reading for
Kentucky and Tennessee.
Also there is a more *generic* ver of the code on my web page if anybody
ever needs it in the future.
You can also use range( MC.pH, MV.pH, na.rm=TRUE).
On Tue, Apr 24, 2012 at 1:29 PM, David Doyle kydaviddo...@gmail.com wrote:
Hi Greg,
Sloved my own problem.
I had some missing data NA in the datasets. So I manually entered the
ylim=range(4,6)
and it worked!!!
Thanks!!
David
On Tue,
Assuming that you want event as the x-axis (horizontal) you can do
something like (untested without reproducible data):
par(mfrow=c(2,1))
scatter.smooth( event, pH1 )
scatter.smooth( event, pH2 )
or
plot( event, pH1, ylim=range(pH1,pH2) , col='blue')
points( event, pH2, col='green' )
lines(
Hi Michael,
What I run into is that I have a lot of time-series data from groundwater
wells but when plot it..it becomes way to busy to see what is going on.
For example, here is the code to load and plot my Mammoth Cave rainfall pH
data
#Load data from web page
data -
Hi Greg,
Thanks,
I got the 1st example to work using the following code:
data - read.csv(http://doylesdartden.com/Monthly-pH-example.csv;, sep=,)
attach(data)
par(mfrow=c(2,1))
scatter.smooth( Year, MC.pH )
scatter.smooth( Year, MV.pH )
This is good but what I'm really looking for is to
Hi Greg,
Sloved my own problem.
I had some missing data NA in the datasets. So I manually entered the
ylim=range(4,6)
and it worked!!!
Thanks!!
David
On Tue, Apr 24, 2012 at 1:55 PM, David Doyle kydaviddo...@gmail.com wrote:
Hi Greg,
Thanks,
I got the 1st example to work using the
Hi folks.
If I have the following in my data
eventpH1pH2
14.0 6.0
24.3 5.9
34.1 6.1
44.0 5.9
and on and on. for about 400 events
Is there a way I can get R to plot event vs. pH1 and event vs. pH2 and
then do a loess or
The scatter plot is easy:
plot(pH1 ~ pH2, data = OBJ)
When you say a loess for each -- how do you break them up? Are there
repeat values for pH1? If so, this might be hard to do in base
graphics, but ggplot2 would make it easy:
library(ggplot2)
ggplot(OBJ, aes(x = pH1, y = pH2)) + geom_point()
8 matches
Mail list logo