I had to place it right before the title and abline functions like so: for (i in 1:length(Participant_Word_Task)){ success_series <- xts(filter(Participant_Word_Task[[i]], GameEndReason == "TIMER_UP")$NumberOfSuccesfulWords , order.by=as.POSIXct(filter(Participant_Word_Task[[i]], GameEndReason == "TIMER_UP")$DeviceTime)) original_series <- xts(Participant_Word_Task[[i]]$NumberOfSuccesfulWords, order.by=as.POSIXct(Participant_Word_Task[[i]]$DeviceTime)) success_decomp <- ts(success_series, frequency = nweeks(success_series)) original_decomp <- ts(original_series, frequency = nweeks(success_series)) WordFrame[i,1] <- unique(Participant_Word_Task[[i]]$StudyId) WordFrame[i,2] <- "Boggle" WordFrame[i,3] <- MannKendall(success_decomp)$tau[1] WordFrame[i,4] <- MannKendall(success_decomp)$sl[1] WordFrame[i,5] <- sens.slope(success_decomp)$estimates WordFrame[i,6] <- sens.slope(success_decomp)$p.value WordFrame[i,7] <- pettitt.test(success_decomp)$estimate WordFrame[i,8] <- pettitt.test(success_decomp)$p.value WordFrame[i,10] <- difftime(time(success_series[WordFrame[i,7]]),time(original_series[1])) WordFrame[i,11] <- tail(which(grepl(success_series[WordFrame[i,7]], original_series)), n=1) WordFrame[i,12] <- sum(Participant_Word_Task[[i]]$TotalDuration[1:WordFrame[i,11]])-sum(Participant_Word_Task[[i]]$TotalTimePaused[1:WordFrame[i,11]]) WordFrame[i,13] <- difftime(time(original_series[length(original_series)]),time(original_series[1])) WordFrame[i,14] <- length(original_series) WordFrame[i,15] <- sum(Participant_Word_Task[[i]]$TotalDuration[1:length(original_series)])-sum(Participant_Word_Task[[i]]$TotalTimePaused[1:length(original_series)])
simplemovingaverage <- SMA(original_series, n = nweeks(original_series)) if (WordFrame[i,4] <= 0.05 & WordFrame[i,6] <= 0.05 & WordFrame[i,8] <= 0.05){ lines(simplemovingaverage,type='l',col='red'); title(paste("Word Task Acquisition for Subject", WordFrame[i,1])); abline(v = index(original_series[WordFrame[i,7]]),lty=2, col='green', lwd=3); dev.off() { pdf(file = paste(WordFrame[i,1],"_Word_Task_Acquisition.pdf",collapse = NULL, sep = "")) plout <- plot(original_series,type='l',col='blue',xlab="Date of Play",ylab="Number of Successful Words") lines(simplemovingaverage,type='l',col='red') print(plout) title(paste("Word Task Acquisition for Subject", WordFrame[i,1])) abline(v = index(original_series[WordFrame[i,7]]),lty=2, col='green', lwd=3) dev.off() } WordFrame[i,18] <- T WordFrame[i,16] <- (1-(WordFrame[i,10]/WordFrame[i,13])) WordFrame[i,17] <- (1-(WordFrame[i,12]/WordFrame[i,15])) } else { { pdf(file = paste(WordFrame[i,1],"_Word_Task_Acquisition.pdf",collapse = NULL, sep = "")) plout <- plot(original_series,type='l',col='blue',xlab="Date of Play",ylab="Number of Successful Words") lines(simplemovingaverage,type='l',col='red') print(plout) title(paste("Word Task Acquisition for Subject", WordFrame[i,1])) dev.off() } WordFrame[i,18] <- F WordFrame[i,16] <- 0 WordFrame[i,17] <- 0 } } But as soon as I did that it worked! Thank you so much! ________________________________ From: Joshua Ulrich <josh.m.ulr...@gmail.com> Sent: Wednesday, October 5, 2022 9:16 PM To: Deramus, Thomas Patrick <tdera...@partners.org> Cc: R-Help <r-help@r-project.org> Subject: Re: [R] Getting "Error in ect, plot.new has not been called yet" despite grouping plot call External Email - Use Caution Hi, My hunch is that you need to add print(plout) before you call dev.off(). See https://stackoverflow.com/a/39853861<https://secure-web.cisco.com/1hPT2Maoird0h5FIvJrkRlN2QFkkujjs4YnHgtB5wYLuhRcHZGZ5SmvPC6x4VQ-rT2HqRBreJ_YeWM8Ug8HtScdp7g_E7nS4BfkHFtjiAKRId_AjYMKIvoJa-YojoOECa97gbnLx2QU535CBNKtd0Z-dzQns8z-C4Du8Jh6K3rEvmlP25DscJvrNcRUbisHQ67um6b6a-H9bqYYN-7QtBp2OGBjnh4LBzDR3TnMNVEF0dbhYTuP5zgP2NsHtCvRrwRUUApS1ZqUT0iiiQCrto3zWrSVlQ6aVOqTBgneUBEk_-8SkwVLgOccgGZyw37JTo/https%3A%2F%2Fstackoverflow.com%2Fa%2F39853861> Try that and let me know if that works. If not, I'll take a closer look later. Best, Josh On Wed, Oct 5, 2022, 1:40 AM Deramus, Thomas Patrick <tdera...@partners.org<mailto:tdera...@partners.org>> wrote: Sorry to cross-post on Stackoverflow and here but I'm having some difficulty. https://stackoverflow.com/questions/73942794/still-getting-error-in-ect-plot-new-has-not-been-called-yet-despite-grouping<https://secure-web.cisco.com/1Xaq1EClv5yIy5MQnXyJBRP44YCvR9TkouPyPXmenkL1IktK5tDxQ_MynVlVIbnlxMlCViUOwZZPfM6GNcEI1NlAqV7y-6RhsR6Qy1m90ENv5mp19zbOX-f-IqNkbFHuGZZLCBFkq0sMS_waWZZSj3Su0bNo5Y1b3zCkVoVLl1Go2J2A1jWX3fvSvOP2ArzpHWt9qX8HZpoGVkadBfUmANamoUw5ucuBsXvvaWdl1skYiJzSVyLOaudIgsCzzRq71RtSGcI0wu7pqSjecwy87CNo2IwYi4hzHW45NiTjzbVUsvvb2AXC4G3Ctoh9c15CM/https%3A%2F%2Fstackoverflow.com%2Fquestions%2F73942794%2Fstill-getting-error-in-ect-plot-new-has-not-been-called-yet-despite-grouping> Trying to make a nested loop that produces PDFs off different graphs, one for ACF/PACF diagnostics and another for the actual data, based on some time-series analyses I'm doing. Unfortunately, I keep getting the dreaded: Error plot.new has not been called yet The code is meant to write a PDF containing the ACF and PACF graphs, then do some analyses on the timeseries, and then make a separate PDF containing a plot describing the timeseries based on the p-values of each test for each individual. library(plyr) library(dplyr) library(ggplot2) library(Kendall) library(lubridate) library(xts) library(TTR) library(trend) library(forecast) library(openxlsx) Game_Metrics_Word_Task <- read.xlsx("GamePack_Analytics_ALL_TIME_Short.xlsx", "Boggle") Game_Metrics_Word_Task <- Game_Metrics_Word_Task %>% filter(grepl('1440', StudyId)) Game_Metrics_Word_Task$DeviceTime <- ymd_hms(Game_Metrics_Word_Task$DeviceTime, tz = "America/New_York") Game_Metrics_Word_Task <- Game_Metrics_Word_Task[!duplicated(Game_Metrics_Word_Task[1:2,])] Participant_Word_Task <- split(arrange(Game_Metrics_Word_Task,StudyId,DeviceTime), arrange(Game_Metrics_Word_Task,StudyId,DeviceTime,StudyId,DeviceTime)$StudyId) WordFrame <- data.frame(Participant = c(0), Task = c(0), MannKendall_Tau = c(0), MannKendall_P = c(0), Sen_Slope_Value = c(0), Sen_Slope_Pval = c(0), Pettitts_CIV = c(0), Pettitts_Pval = c(0), ARIMA_Model = c(0), Time_to_Petit = c(0), Number_of_Trials_to_Pettitt = c(0), Playtime_to_Petit_seconds = c(0), Time_Start_to_end_days = c(0), Number_of_Total_Trials = c(0), Total_Playtime_seconds = c(0), Learning_rate_days = c(0), Learning_rate_seconds = c(0), Learned_Task = c(0)) for (i in 1:length(Participant_Word_Task)){ success_series <- xts(filter(Participant_Word_Task[[i]], GameEndReason == "TIMER_UP")$NumberOfSuccesfulWords , order.by<http://secure-web.cisco.com/10boIOil_1j37yFHh_eAQC6S-GAWujDCUxF8SFXlL05I4ABkjAjfCpv636hC2_D9ZcnRJNo07NJLs1husggxCVibrINE_e4DkPNcX0JYVosd5gut0FGPSMyIREh9RQoZW3GPaCVyAjIqzeAuKYHJz1C9DZqEdN-SFLGFqZaWzvaYwep785sKgQJCb4FlZA7G7jjr42MjSg97lgbRnFXj7vopsZ3gXb_EoSNtkIHNP7y5bMYQio505DE0WHrz1FFFbt_WehSBrNEAx5UXdrD5gkNwp1jzLMONBJA_2Dj3gZUhAP1AQOyHrI_qYEQKtFpiO/http%3A%2F%2Forder.by>=as.POSIXct(filter(Participant_Word_Task[[i]], GameEndReason == "TIMER_UP")$DeviceTime)) original_series <- xts(Participant_Word_Task[[i]]$NumberOfSuccesfulWords, order.by<http://secure-web.cisco.com/10boIOil_1j37yFHh_eAQC6S-GAWujDCUxF8SFXlL05I4ABkjAjfCpv636hC2_D9ZcnRJNo07NJLs1husggxCVibrINE_e4DkPNcX0JYVosd5gut0FGPSMyIREh9RQoZW3GPaCVyAjIqzeAuKYHJz1C9DZqEdN-SFLGFqZaWzvaYwep785sKgQJCb4FlZA7G7jjr42MjSg97lgbRnFXj7vopsZ3gXb_EoSNtkIHNP7y5bMYQio505DE0WHrz1FFFbt_WehSBrNEAx5UXdrD5gkNwp1jzLMONBJA_2Dj3gZUhAP1AQOyHrI_qYEQKtFpiO/http%3A%2F%2Forder.by>=as.POSIXct(Participant_Word_Task[[i]]$DeviceTime)) success_decomp <- ts(success_series, frequency = nweeks(success_series)) original_decomp <- ts(original_series, frequency = nweeks(success_series)) pdf(paste("Word_Task_Autocorrelation_plots_for_subject_",unique(Participant_Word_Task[[i]]$StudyId),".pdf" ,collapse = NULL, sep = "")) par(mfrow=c(1,2)) Autocorrelationplot <- acf(success_decomp, main=paste("")) PartialAutocorrelationplot <- pacf(success_decomp, main=paste("")) mtext(paste("Word Task Auto and Partialauto correlations for subject ",unique(Participant_Word_Task[[i]]$StudyId)), side = 3, line = -3, outer = TRUE) dev.off() AutomatedArimaoutputs <- auto.arima(success_decomp) p <- length(AutomatedArimaoutputs$model$phi) #AR component q <- length(AutomatedArimaoutputs$model$theta) #MA component d <- AutomatedArimaoutputs$model$Delta #order of difference WordFrame[i,1] <- unique(Participant_Word_Task[[i]]$StudyId) WordFrame[i,2] <- "Word" WordFrame[i,3] <- MannKendall(success_decomp)$tau[1] WordFrame[i,4] <- MannKendall(success_decomp)$sl[1] WordFrame[i,5] <- sens.slope(success_decomp)$estimates WordFrame[i,6] <- sens.slope(success_decomp)$p.value WordFrame[i,7] <- pettitt.test(success_decomp)$estimate WordFrame[i,8] <- pettitt.test(success_decomp)$p.value WordFrame[i,9] <- paste("ARIMA(",p,",",q,",",d,")", collapse = NULL, sep = "") WordFrame[i,10] <- difftime(time(success_series[WordFrame[i,7]]),time(original_series[1])) WordFrame[i,11] <- tail(which(grepl(success_series[WordFrame[i,7]], original_series)), n=1) WordFrame[i,12] <- sum(Participant_Word_Task[[i]]$TotalDuration[1:WordFrame[i,11]])-sum(Participant_Word_Task[[i]]$TotalTimePaused[1:WordFrame[i,11]]) WordFrame[i,13] <- difftime(time(original_series[length(original_series)]),time(original_series[1])) WordFrame[i,14] <- length(original_series) WordFrame[i,15] <- sum(Participant_Word_Task[[i]]$TotalDuration[1:length(original_series)])-sum(Participant_Word_Task[[i]]$TotalTimePaused[1:length(original_series)]) simplemovingaverage <- SMA(original_series, n = nweeks(original_series)) if (WordFrame[i,4] <= 0.05 & WordFrame[i,6] <= 0.05 & WordFrame[i,8] <= 0.05){ { pdf(paste(WordFrame[i,1],"_Word_Task_Acquisition.pdf",collapse = NULL, sep = "")) plout <- plot(original_series,type='l',col='blue',xlab="Date of Play",ylab="Number of Successful Words") lines(simplemovingaverage,type='l',col='red') title(paste("Word Task Acquisition for Subject", WordFrame[i,1])) abline(v = index(original_series[WordFrame[i,7]]),lty=2, col='green', lwd=3) dev.off() } WordFrame[i,18] <- T WordFrame[i,16] <- (1-(WordFrame[i,10]/WordFrame[i,13])) WordFrame[i,17] <- (1-(WordFrame[i,12]/WordFrame[i,15])) } else { { pdf(paste(WordFrame[i,1],"_Word_Task_Acquisition.pdf",collapse = NULL, sep = "")) plout <- plot(original_series,type='l',col='blue',xlab="Date of Play",ylab="Number of Successful Words") lines(simplemovingaverage,type='l',col='red') title(paste("Word Task Acquisition for Subject", WordFrame[i,1])) dev.off() } WordFrame[i,18] <- F WordFrame[i,16] <- 0 WordFrame[i,17] <- 0 } } It will work just fine if I run the lines individually (e.g. set i = 1, 2, ect), and if I comment out abline and title (lines seems to work fine). But it will throw the error everytime I try to run the loop without these commented. Have tried just about everything I could find on the Stack forums to run everything as a single argument and I'm just not sure what is wrong with it. dev.list() spits out: pdf 2 following the error. With abline and title commented out and lines run individually it's NULL. Happens in both RStudio 2022.07.2+576 "Spotted Wakerobin" Release (e7373ef832b49b2a9b88162cfe7eac5f22c40b34, 2022-09-06) for Ubuntu Bionic Mozilla/5.0 (X11; Linux x86_64) AppleWebKit/537.36 (KHTML, like Gecko) QtWebEngine/5.12.8 Chrome/69.0.3497.128 Safari/537.36 And R: platform x86_64-pc-linux-gnu arch x86_64 os linux-gnu system x86_64, linux-gnu status major 4 minor 2.1 year 2022 month 06 day 23 svn rev 82513 language R version.string R version 4.2.1 (2022-06-23) nickname Funny-Looking Kid My OS: PRETTY_NAME="Debian GNU/Linux 11 (bullseye)" NAME="Debian GNU/Linux" VERSION_ID="11" VERSION="11 (bullseye)" VERSION_CODENAME=bullseye ID=debian HOME_URL="https://www.debian.org/<https://secure-web.cisco.com/1ACDhN-pLD9mzYG0v4VTzMKxAyLmTIADEHLJIuzQBX8KrUHWm1giGPYPmZCHAPpLY42jxDLfqh6eM_dWl9xJg-NmLD-dzpJqBQ-Yy3WuhZONfcn8ZwoK19G46Dm_hBu-Hrl9eV2uzZsEDY_38eQoRMVILy98SEEufRocF9YyjPAW8Y-euSSmPOLDK2PF9CGMABP9qNF6O_sYhykJnPf2_2ptZHvGAEuDDwY925wnRVsvLdc7Rl5zSJE5n45yGOXFaeihwRearfnxKPi8vufhoDC0LOKRYIVSTRAPJJEq3dqTgtDJXkq4UKyJToMfCtF0x/https%3A%2F%2Fwww.debian.org%2F>" SUPPORT_URL="https://www.debian.org/support<https://secure-web.cisco.com/1Ck4yTAQD273G1LHpBEXmkgIu4kavtz3379Nq-aX3Btcot2VBSq9lKAZN5X2tWi6yT-GMfbuJYp3U621mE6Vk4vf9MJQ3SNMkBn19QcLWgLVnsQvOupmnVd31EWpgYIsutKyJ6Gro4S47Ac8b0GQLFE8VZNrWlbs5YKY3hcPuctYD13LVyjL5zlhP5zvg8wbUX73GBYMHp94c2qKPUfNnsUH1VxXudeB4ZtlnPl66TcO8gJAIYdQmVlw0tIOTdr5VVJ6EbdrRRHKJ7AJvsJWMIn8SAHHusAf7Kc8kCGc0obyScRr-5hU-ImcSro2vFhUe/https%3A%2F%2Fwww.debian.org%2Fsupport>" BUG_REPORT_URL="https://bugs.debian.org/<https://secure-web.cisco.com/1Kma-rbfMkalSm7CWHd2_J7JxCIrLBWgbW7-6BMP4GaNfBdSR7Y1Cb-8OMMnZjLZbNWg60FnjHKbHrwdsX_noY5nvLC_EaTbhL883q3UtHb3j3i0-Ebdd1vMvHvgn-xWeuYwRTmIZjWV8elT_g_phphQB7EMdaPCS6lHSGC1wkpOKYg9V_cNlEHHMsrtuRMXD-kkjj9h_ziH6Nu-25zVvTIfUbr0ODIxrtThOkJ4wopWByCBXPKInKG_-1gBNV8q_oPar_NJDH3N6D8Zil9T2dPkgZrSUkvP42xw2XnvQ_fYD1x9_57P8KPtPpCik1Sdc/https%3A%2F%2Fbugs.debian.org%2F>" No LSB modules are available. Distributor ID: Debian Description: Debian GNU/Linux 11 (bullseye) Release: 11 Codename: bullseye Icon name: computer-desktop Chassis: desktop Machine ID: 053ebf23707f49c8ad4e0684f4cf19d3 Boot ID: d0e6294d3b944286bef10e76c21e6401 Operating System: Debian GNU/Linux 11 (bullseye) Kernel: Linux 5.10.0-18-amd64 Architecture: x86-64 Any suggestions would be greatly appreciated. -- Thomas DeRamus (He/Him/His) Data Analyst Massachusetts General Hospital Brigham Alzheimer’s Clinical & Translational Research Unit 149 13th Street Charlestown, MA 02129 Phone: 205-834-5066 Email: tdera...@partners.org<mailto:tdera...@partners.org><mailto:tdera...@partners.org<mailto:tdera...@partners.org>>, tpdera...@gmail.com<mailto:tpdera...@gmail.com><mailto:tpdera...@gmail.com<mailto:tpdera...@gmail.com>> [https://ci3.googleusercontent.com/mail-sig/AIorK4we2sU30P2HyfDQF5hpEjYTt-9FTBK7cAVsP7EenrZ0nsKCf48fuYMtElj6Szn_2fpSPWr66eQ][https://ci3.googleusercontent.com/mail-sig/AIorK4yyY0DlImU0UONJrHTbPc5T3lJj8Kmu8SbDKJJ3XjcX6CgvVsvSueYKwficYFz4zXt6fZV8YIY<https://secure-web.cisco.com/1p4CoOEfGfFgyADIKdZuTyxCAxu1iAL4dl_xBOs0kNACeZBlKqxOU1hdF4bmbj_pz-ZhvRQyetDE_Br_g_kCZNaTln8L0uSdieFdFVX6W2p81hni5pzndOTrQPuenlZezHkcZDaXhr566__QDE7iZutarcm6bGONz_PdpQMeFGoqy9DtDbvmJLEf8x4C_4H1cFd1xYpj4sXozZ1vhTpbayy4IHfK9qLHsnzyPWDEtccLFMO2DjOQu1Kguw_rBvVlQlef7u_XufAbLPWJ5c6FUUC8K1vLFg1xB4-YJctdg4he3hvvz_e8kue-C0qn7a_Iw/https%3A%2F%2Fci3.googleusercontent.com%2Fmail-sig%2FAIorK4we2sU30P2HyfDQF5hpEjYTt-9FTBK7cAVsP7EenrZ0nsKCf48fuYMtElj6Szn_2fpSPWr66eQ%255D%255Bhttps%3A%2F%2Fci3.googleusercontent.com%2Fmail-sig%2FAIorK4yyY0DlImU0UONJrHTbPc5T3lJj8Kmu8SbDKJJ3XjcX6CgvVsvSueYKwficYFz4zXt6fZV8YIY>] “If knowledge can create problems, it is not through ignorance that we can solve them.” —Issac Asimov The information in this e-mail is intended only for the person to whom it is addressed. If you believe this e-mail was sent to you in error and the e-mail contains patient information, please contact the Mass General Brigham Compliance HelpLine at https://www.massgeneralbrigham.org/complianceline<https://secure-web.cisco.com/1nLR3j5jLEy6A5yyLXTgiSEPglT5zpYBKY0XNFwmEgtrG9PV13Lfd015mRGeFTVsVSFliuG7VXmI-O3Fw0DpFvFKA0GR9lj4jtJ3BQL8MblvkUTO-F6rUjk5jgtDZtSmOeRU6wInl00IqoCwYR3FSQD9_XBe3uQzm6mAviB2SqvIE4u4BFQbNnhQSh-z6E6eGVGAIuG64OGuTByh1AYgf2O3xP3PcXwQiUa6M504snWIaHca1eJ7HB03HRx91JqKUnwhXhOMXidWvZ8djcHw59faZKzc68R55-8RkBTvcDJ1Ln_5ouA6Z0VgC2oHGWPkW/https%3A%2F%2Fwww.massgeneralbrigham.org%2Fcomplianceline> <https://www.massgeneralbrigham.org/complianceline<https://secure-web.cisco.com/1nLR3j5jLEy6A5yyLXTgiSEPglT5zpYBKY0XNFwmEgtrG9PV13Lfd015mRGeFTVsVSFliuG7VXmI-O3Fw0DpFvFKA0GR9lj4jtJ3BQL8MblvkUTO-F6rUjk5jgtDZtSmOeRU6wInl00IqoCwYR3FSQD9_XBe3uQzm6mAviB2SqvIE4u4BFQbNnhQSh-z6E6eGVGAIuG64OGuTByh1AYgf2O3xP3PcXwQiUa6M504snWIaHca1eJ7HB03HRx91JqKUnwhXhOMXidWvZ8djcHw59faZKzc68R55-8RkBTvcDJ1Ln_5ouA6Z0VgC2oHGWPkW/https%3A%2F%2Fwww.massgeneralbrigham.org%2Fcomplianceline>> . Please note that this e-mail is not secure (encrypted). If you do not wish to continue communication over unencrypted e-mail, please notify the sender of this message immediately. Continuing to send or respond to e-mail after receiving this message means you understand and accept this risk and wish to continue to communicate over unencrypted e-mail. ______________________________________________ R-help@r-project.org<mailto:R-help@r-project.org> mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help<https://secure-web.cisco.com/169zqX-6K1A_ERyLXXpVduo-L_bdP58wU0tW6976VBgCPzCLprqR3-s8od_SWvcEnJWclWMXLbDu8eHwBNlKcyrCpUuuMxKuz3xGbpVjkqDLLHD7C45e3EKTiUq5JIEFLQ8XvR-IsvhuHJjqfmHKzaeEzciR71ac-nMIG0ZJHmMuQyVv3Yuety2ReZ0Xa5vdHW5zcq0HUKH16bzLXeqSWc6rRkE-i-y9WK_OX959hz6gDnIMSu8RNuwe0TndTxMQyHcG-fK9Ph0sG_12ZjFdw3qL2yzWsEI3sXZ6mFU2v7npLKE1634X53nJ-FnqFhRQ6/https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-help> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html<http://secure-web.cisco.com/1q--2yVlI7EkpVlbb7cUxHfn03FTKwZsan7XtRgYCrbX98T0WAQn7WNr8Hqt2ZVUuBgmAA5i5zh5M8nfzRoj6CFtWvdyzMvxMdS6eXYA3rlGcZ8LsHDHoUQbsEE9ui-hVicew7e5Ocd1Ed5OSDVeGDecKIEjUUZwTmS3DDa_ilt8dNQOGZaBGOdbnjk0bNYRDT_FiB077ShZRBmKmTX8838zgtFk-hVIGcfvbywCmhpGibHuQ7zyZhd-5wgj9b-c8d2XA6A8yno_2I9DNltls3YoqZcjFoczsuvvTx24DjXJGmaz8wRm3CEg1fuTONU0a/http%3A%2F%2Fwww.R-project.org%2Fposting-guide.html> and provide commented, minimal, self-contained, reproducible code. The information in this e-mail is intended only for the person to whom it is addressed. If you believe this e-mail was sent to you in error and the e-mail contains patient information, please contact the Mass General Brigham Compliance HelpLine at https://www.massgeneralbrigham.org/complianceline <https://www.massgeneralbrigham.org/complianceline> . Please note that this e-mail is not secure (encrypted). If you do not wish to continue communication over unencrypted e-mail, please notify the sender of this message immediately. Continuing to send or respond to e-mail after receiving this message means you understand and accept this risk and wish to continue to communicate over unencrypted e-mail. [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org 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.