Hello guys

I have a melted dataset that has 36 different time series included. With 
stat.summary() I added the mean. The problem is that these 36 time series do 
not have the same length. Therefore I would like to plot the mean only for the 
length that the shortest of these 36 time series has. Example. First time 
series; 28 years long; 2. time series 21 years long and so so. So the mean in 
this example should not be longer than these 21 years. I searched the internet 
but no return so far. Below is the code:

melted_SWT <- melt(DF_filled_SWT,id.vars="X37")

mean_SW <- melt(apply(DF_filled_SWT,1,mean))

stat_summary(aes(shape="mean",group=1),fun.y = "mean",size = 1, geom = "line", 
show.legend = T)+

XIntercept <- data.frame( Intercept=c(1))
YIntercept <- data.frame( Intercept=c(0))

ggplot(data=melted_SWT,aes(x=X37,y=value,colour=variable)) +  geom_line() +
    ggtitle("� GMST during Hiatuses")+ylab("(�K)") + xlab("timelag")+
  theme(panel.background = element_rect(fill = 'grey87')) +  
stat_summary(aes(shape="mean",group=1),fun.y = "mean",size = 1, geom = "line", 
show.legend = T)+
  theme(axis.title.x = element_text( size=24, face="bold")) + theme(title = 
element_text( size=24, face="bold"))+
  theme(axis.text.x= element_text( size=24, face="bold") ) +  
theme(axis.text.y= element_text( size=24, face="bold") ) +
  geom_vline(aes(xintercept=XIntercept),data=Intercept) + 
geom_hline(aes(yintercept=YIntercept),data=Intercept)+
  theme(legend.text = element_text( size=24, face="bold")) + 
theme(legend.margin = unit(0.7, "line")) +scale_x_continuous(breaks = 
c(-10,1,10,20), labels = c(-10,1,10,20))+
  theme(legend.key.height=unit(2.6,"line")) + theme(legend.key.size = unit(2.3, 
"cm"))



data :  X37 variable      value
1 -10       X1 -0.1940203
2  -9       X1 -0.1578422
3  -8       X1 -0.1658450
4  -7       X1 -0.1873830
5  -6       X1 -0.1471910
6  -5       X1 -0.1160171


Thank for the help

        [[alternative HTML version deleted]]

______________________________________________
[email protected] mailing list -- To UNSUBSCRIBE and more, see
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.

Reply via email to