Um texto embutido e sem conjunto de caracteres especificado associado... Nome: n?o dispon?vel Url: https://stat.ethz.ch/pipermail/r-help/attachments/20070508/1c7def0a/attachment.pl
You have to have a valid 'date' object on the x-axis. Try this:
simulation<-data.frame(cbind(my.year=c(rep(2000,8),rep(2001,12),rep(2002,12)),my.month=c(5:12,1:12,1:12)))
simulation$year.month<-paste(simulation$my.year,"_",ifelse(simulation$my.month>=10,simulation$my.month,paste("0",simulation$my.month,sep="")),sep="")
simulation$freq<-sample(1:40,32)
# create POSIXct time
simulation$time <- ISOdate(simulation$my.year, simulation$my.month,1)
attach(simulation)
plot(time, freq)
On 5/8/07, Milton Cezar Ribeiro <milton_ruser at yahoo.com.br>
wrote:> Dear all,
>
> I have two data frame, on with a complete list of my field survey with
frequency data of a sample species. This data frame looks like:
>
>
>
simulation<-data.frame(cbind(my.year=c(rep(2000,8),rep(2001,12),rep(2002,12)),my.month=c(5:12,1:12,1:12)))
>
simulation$year.month<-paste(simulation$my.year,"_",ifelse(simulation$my.month>=10,simulation$my.month,paste("0",simulation$my.month,sep="")),sep="")
> simulation$freq<-sample(1:40,32)
> attach(simulation)
> plot(year.month, freq)
>
> As you can see, I have a collumn with the year and month of my samples, and
a freq variable with simulated data. I would like to plot this data but when I
try to use the plot showed above, I get a error message.
>
> After bypass this problem, I would like add points in my graph with
simulated data for only a random number of survey month, but I need that the
full range of surveys be kept on the X-axis. Just to simulate a sample I am
using:
>
> simulation.sample<-simulation[sample(1:length(year.month),8,
replace=F),]
> simulation.sample$freq<-sample(1:40,8)
>
> Any ideas?
>
> Kind regards
>
> Miltinho
>
> __________________________________________________
>
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at stat.math.ethz.ch 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.
>
--
Jim Holtman
Cincinnati, OH
+1 513 646 9390
What is the problem you are trying to solve?
- use proper spacing to make it easier to read - start off with set.seed to make it reproducible - omit cbind and combine all the rep's into one rep in first line - make the date column a known date class (here "Date"), set.seed(1) sim <- data.frame( my.year = rep(2000:2002, c(8, 12, 12)), my.month = c(5:12, 1:12, 1:12), freq = sample(1:40, 32) ) sim$date <- as.Date(paste(sim$my.year, sim$my.month, 1, sep = "-")) plot(freq ~ date, sim) On 5/8/07, Milton Cezar Ribeiro <milton_ruser at yahoo.com.br> wrote:> Dear all, > > I have two data frame, on with a complete list of my field survey with frequency data of a sample species. This data frame looks like: > > > simulation<-data.frame(cbind(my.year=c(rep(2000,8),rep(2001,12),rep(2002,12)),my.month=c(5:12,1:12,1:12))) > simulation$year.month<-paste(simulation$my.year,"_",ifelse(simulation$my.month>=10,simulation$my.month,paste("0",simulation$my.month,sep="")),sep="") > simulation$freq<-sample(1:40,32) > attach(simulation) > plot(year.month, freq) > > As you can see, I have a collumn with the year and month of my samples, and a freq variable with simulated data. I would like to plot this data but when I try to use the plot showed above, I get a error message. > > After bypass this problem, I would like add points in my graph with simulated data for only a random number of survey month, but I need that the full range of surveys be kept on the X-axis. Just to simulate a sample I am using: > > simulation.sample<-simulation[sample(1:length(year.month),8, replace=F),] > simulation.sample$freq<-sample(1:40,8) > > Any ideas? > > Kind regards > > Miltinho > > __________________________________________________ > > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at stat.math.ethz.ch 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 understands only numerical and Date class values for axis. So either
a) plot them using the sequence 1, ..., 32 and then explicitly label
them. Here is an example:
n <- length(year.month)
plot( 1:n, freq, xaxt="n")
mtext( text=year.month, side=1, at=1:n, las=2 )
b) or create the dates in Date format. This option is preferable if the
dates were varying unequally.
x <- seq( as.Date("2000-05-01"), as.Date("2002-12-01"),
by="1 month" )
plot(x, simulation$freq)
BTW, you could also have created year.month via
paste( rep( 2000:2002, c(8,12,12) ),
formatC( c(5:12,1:12,1:12), width=2, flag="0" ) ,
sep="_" )
Regards, Adai
Milton Cezar Ribeiro wrote:> Dear all,
>
> I have two data frame, on with a complete list of my field survey with
frequency data of a sample species. This data frame looks like:
>
>
>
simulation<-data.frame(cbind(my.year=c(rep(2000,8),rep(2001,12),rep(2002,12)),my.month=c(5:12,1:12,1:12)))
>
simulation$year.month<-paste(simulation$my.year,"_",ifelse(simulation$my.month>=10,simulation$my.month,paste("0",simulation$my.month,sep="")),sep="")
> simulation$freq<-sample(1:40,32)
> attach(simulation)
> plot(year.month, freq)
>
> As you can see, I have a collumn with the year and month of my samples, and
a freq variable with simulated data. I would like to plot this data but when I
try to use the plot showed above, I get a error message.
>
> After bypass this problem, I would like add points in my graph with
simulated data for only a random number of survey month, but I need that the
full range of surveys be kept on the X-axis. Just to simulate a sample I am
using:
>
> simulation.sample<-simulation[sample(1:length(year.month),8,
replace=F),]
> simulation.sample$freq<-sample(1:40,8)
>
> Any ideas?
>
> Kind regards
>
> Miltinho
>
> __________________________________________________
>
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at stat.math.ethz.ch 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.
>
>
>