Re: [R] geom_smooth with sd

2024-08-11 Thread SIBYLLE STÖCKLI via R-help
Thanks Erin Quite relevant. Yes now +sd and -sd are the same values. However they are about +/- 5 and not the values received by the simple code here. I still think it is as the length of y differs. Simple statistics > mean(MS2020[MS2020$Bio=="1",]$QI_A, na.rm=TRUE) [1] 26.81225 > sd(M

Re: [R] geom_smooth with sd

2024-08-11 Thread Erin Hodgess
Hi! This is probably completely off base, but your ymin and y max setup lines are different. One uses sqrt(y), while the second uses sqrt(length(y)). Could that play a part, please? Thank you Erin Hodgess, PhD mailto: erinm.hodg...@gmail.com On Sun, Aug 11, 2024 at 10:10 AM SIBYLLE STÖCKLI

[R] geom_smooth with sd

2024-08-11 Thread SIBYLLE STÖCKLI via R-help
Dear community Using after_stat() I was able to visualise ggplot with standard deviations instead of a confidence interval as seen in the R help. p1<-ggplot(data = MS1, aes(x= Jahr, y= QI_A,color=Bio, linetype=Bio)) + geom_smooth(aes(fill=Bio, ymax=after_stat(y+se*sqrt(len

Re: [R] Printing

2024-08-11 Thread Rui Barradas
Às 15:36 de 11/08/2024, Steven Yen escreveu: Thanks. Will try it. Have not tried it but I think the following may work: out$results<-NULL out$results$ei<-ap out$results$vi<-vap All I need is printing by returning out (unless I turn it off). And, retrieve ap and vap as needed as shown above.

Re: [R] Printing

2024-08-11 Thread Bert Gunter
Just for reference, R's "preferred" way to handle printing is to create "classed" objects from a function that produces all the necessary information and have an appropriate "print" method, i.e. function, for that class to print it. R has several class object systems that can do this, including the

Re: [R] Printing

2024-08-11 Thread Ivan Krylov via R-help
В Sun, 11 Aug 2024 22:36:08 +0800 Steven Yen пишет: > All I need is printing by returning out (unless I turn it off). And, > retrieve ap and vap as needed as shown above. Guess I need to read > more about invisible. Perhaps you could print(out) instead of returning it in the if (printing) branc

Re: [R] Printing

2024-08-11 Thread Steven Yen
Thanks. Will try it. Have not tried it but I think the following may work: out$results<-NULL out$results$ei<-ap out$results$vi<-vap All I need is printing by returning out (unless I turn it off). And, retrieve ap and vap as needed as shown above. Guess I need to read more about invisible.

Re: [R] Printing

2024-08-11 Thread Rui Barradas
Às 09:51 de 11/08/2024, Steven Yen escreveu: Hi In the following codes, I had to choose between printing (= TRUE) or deliver something for grab (ei, vi). Is there a way to get both--that is, to print and also have ei and vi for grab? Thanks. Steven ... out<-round(as.data.frame(cbind(ap,se,

[R] Printing

2024-08-11 Thread Steven Yen
Hi In the following codes, I had to choose between printing (= TRUE) or deliver something for grab (ei, vi). Is there a way to get both--that is, to print and also have ei and vi for grab? Thanks. Steven ... out<-round(as.data.frame(cbind(ap,se,t,p)),digits) out<-cbind(out,sig) out<-out[!gr