# ================================================================
# TUGAS BIOSTATISTIKA INTERMEDIATE
# NAMA : DIMAS WISNU WARDHANA
# NIM  : 2611018042
# ANALISIS PENGUKURAN BERULANG BERAT BADAN
# SUMBER: Ji et al. (2026) - Bofanglutide Phase 2b
#
# Artikel:
# https://doi.org/10.1038/s41392-026-02586-8
# PubMed:
# https://pubmed.ncbi.nlm.nih.gov/41760612/
# ClinicalTrials.gov: NCT06256562
#
# PENTING:
# Dataset individu yang dianalisis di sini adalah DATA SIMULASI untuk
# pembelajaran. Struktur kelompok dan pola longitudinal dikalibrasi
# terhadap artikel. Nilai week-30 mengikuti LSM publikasi; titik antara
# merupakan pendekatan didaktik terhadap trajectory Figure 2a.
# Dataset ini BUKAN raw data pasien asli.
# ================================================================

FILE_DATA <- "Bofanglutide_2026_Simulated_Long.csv"
FOLDER_HASIL <- "hasil_analisis_bofanglutide_2026"
KELOMPOK_FOKUS <- "Bofanglutide 24 mg QW"
INSTALL_PAKET_OTOMATIS <- TRUE
TAMPILKAN_GRAFIK_RSTUDIO <- TRUE
BUKA_HTML_OTOMATIS <- TRUE

jalankan_analisis <- function() {

  waktu_mulai <- Sys.time()

  paket <- c(
    "dplyr","tidyr","ggplot2","afex","emmeans","rstatix","car",
    "effectsize","lme4","lmerTest","pbkrtest","performance",
    "ggpubr","knitr","htmltools","base64enc","tibble"
  )

  belum <- paket[!vapply(paket, requireNamespace, logical(1), quietly = TRUE)]
  if (length(belum)) {
    if (!INSTALL_PAKET_OTOMATIS) stop("Paket belum tersedia: ", paste(belum, collapse=", "))
    install.packages(belum, repos="https://cloud.r-project.org")
  }

  suppressPackageStartupMessages({
    library(dplyr); library(tidyr); library(ggplot2)
    library(afex); library(emmeans); library(rstatix)
    library(car); library(effectsize); library(lme4); library(lmerTest)
    library(performance); library(ggpubr); library(knitr)
    library(htmltools); library(base64enc); library(tibble)
  })

  opsi_lama <- options(contrasts = c("contr.sum","contr.poly"))
  on.exit(options(opsi_lama), add = TRUE)

  # ------------ lokasi file ------------
  sumber <- unlist(lapply(sys.frames(), function(fr) {
    if (!is.null(fr$ofile)) as.character(fr$ofile) else character(0)
  }))
  if (requireNamespace("rstudioapi", quietly=TRUE) && rstudioapi::isAvailable()) {
    aktif <- tryCatch(rstudioapi::getSourceEditorContext()$path, error=function(e) "")
    if (nzchar(aktif)) sumber <- c(sumber, aktif)
  }
  kandidat <- unique(c(FILE_DATA, file.path(dirname(sumber), FILE_DATA)))
  kandidat <- kandidat[file.exists(kandidat)]
  if (!length(kandidat)) stop("File data tidak ditemukan: ", FILE_DATA)
  path_data <- normalizePath(kandidat[1], winslash="/", mustWork=TRUE)

  folder <- file.path(dirname(path_data), FOLDER_HASIL)
  dir.create(folder, recursive=TRUE, showWarnings=FALSE)
  folder <- normalizePath(folder, winslash="/", mustWork=TRUE)

  # ------------ helper HTML/output ------------
  isi <- character(0)
  status <- data.frame()

  esc <- function(x) as.character(htmltools::htmlEscape(paste(x, collapse="\n")))
  tambah <- function(x) isi <<- c(isi, x)
  narasi <- function(x) tambah(paste0("<p>", esc(x), "</p>"))

  output <- function(x) {
    z <- capture.output(print(x))
    cat(paste(z, collapse="\n"), "\n")
    tambah(paste0("<pre>", esc(z), "</pre>"))
    invisible(x)
  }

  tabel <- function(x, nama) {
    x <- as.data.frame(x)
    write.csv(x, file.path(folder, paste0(nama, ".csv")),
              row.names=FALSE, na="", fileEncoding="UTF-8")
    if (nrow(x)) {
      tambah(paste0("<div class='table'>",
                    as.character(knitr::kable(x, format="html", digits=4)),
                    "</div>"))
    }
    print(x)
    invisible(x)
  }

  gambar <- function(p, nama, lebar=10, tinggi=6) {
    path <- file.path(folder, paste0(nama, ".png"))
    ggsave(path, p, width=lebar, height=tinggi, dpi=180)
    if (TAMPILKAN_GRAFIK_RSTUDIO) try(print(p), silent=TRUE)
    tambah(paste0("<figure><img src='",
                  base64enc::dataURI(file=path, mime="image/png"),
                  "'><figcaption>", esc(nama), "</figcaption></figure>"))
    message("Grafik tersimpan: ", path)
  }

  bagian <- function(judul, f) {
    message("\n>>> ", judul)
    tambah(paste0("<h2>", esc(judul), "</h2>"))
    warns <- character(0)
    ok <- tryCatch(
      withCallingHandlers({ f(); TRUE },
        warning=function(w) {
          warns <<- c(warns, conditionMessage(w))
          invokeRestart("muffleWarning")
        }),
      error=function(e) {
        narasi(paste("BAGIAN GAGAL:", conditionMessage(e)))
        message("Bagian gagal: ", conditionMessage(e))
        FALSE
      })
    if (length(warns)) narasi(paste("Peringatan:", paste(unique(warns), collapse="; ")))
    status <<- rbind(status, data.frame(
      bagian=judul,
      status=if(ok) "Selesai" else "Gagal",
      peringatan=paste(unique(warns), collapse="; ")
    ))
  }

  # ------------ data ------------
  dat <- read.csv(path_data, stringsAsFactors=FALSE, check.names=FALSE)

  level_kelompok <- c(
    "Placebo",
    "Bofanglutide 12 mg Q2W",
    "Bofanglutide 18 mg Q2W",
    "Bofanglutide 24 mg Q2W",
    "Bofanglutide 48 mg Q2W",
    "Bofanglutide 24 mg QW"
  )
  level_waktu <- c("W0","W4","W8","W12","W16","W20","W24","W30")
  minggu_nilai <- c(0,4,8,12,16,20,24,30)

  wajib <- c("id","treatment","week","time","baseline_weight_kg","pct_change","weight_kg")
  if (!all(wajib %in% names(dat))) stop("Kolom wajib tidak lengkap.")

  dat <- dat |>
    mutate(
      id=factor(id),
      treatment=factor(treatment, levels=level_kelompok),
      time=factor(time, levels=level_waktu),
      week=as.numeric(week),
      weight_kg=as.numeric(weight_kg),
      pct_change=as.numeric(pct_change)
    ) |>
    arrange(id, week)

  if (anyNA(dat$treatment) || anyNA(dat$time)) stop("Ada label kelompok/waktu tidak dikenali.")
  if (anyDuplicated(dat[c("id","time")])) stop("Duplikasi id-waktu terdeteksi.")

  dat_obs <- dat |> filter(!is.na(weight_kg))
  wide <- dat |> select(id,treatment,time,weight_kg) |>
    pivot_wider(names_from=time, values_from=weight_kg)

  lengkap_id <- wide |> filter(if_all(all_of(level_waktu), ~ !is.na(.x))) |> pull(id)
  cc <- dat |> filter(id %in% lengkap_id)
  d1 <- cc |> filter(treatment == KELOMPOK_FOKUS)

  # ================================================================
  bagian("1. Sumber, desain, dan integritas data", function() {
    narasi("Sumber ilmiah: Ji et al. (2026), Signal Transduction and Targeted Therapy, DOI 10.1038/s41392-026-02586-8, trial NCT06256562.")
    narasi("Analisis ini memakai data individu sintetis yang dibuat untuk pembelajaran. Enam kelompok mengikuti struktur uji fase 2b; delapan titik pengukuran adalah W0, W4, W8, W12, W16, W20, W24, dan W30.")
    narasi("Outcome utama untuk analisis repeated measures adalah berat badan (kg). Persentase perubahan juga dipertahankan sebagai variabel pendukung.")
    tabel(dat |> count(treatment, time, name="n"), "01_jumlah_per_sel")
    miss <- dat |> group_by(treatment,time) |>
      summarise(n=n(), n_tersedia=sum(!is.na(weight_kg)),
                n_hilang=sum(is.na(weight_kg)), .groups="drop")
    tabel(miss, "01_missing")
  })

  bagian("2. Statistik deskriptif", function() {
    desk <- dat |> group_by(treatment,time,week) |>
      summarise(
        n=sum(!is.na(weight_kg)),
        mean=mean(weight_kg,na.rm=TRUE),
        sd=sd(weight_kg,na.rm=TRUE),
        median=median(weight_kg,na.rm=TRUE),
        min=min(weight_kg,na.rm=TRUE),
        max=max(weight_kg,na.rm=TRUE),
        mean_pct=mean(pct_change,na.rm=TRUE),
        sd_pct=sd(pct_change,na.rm=TRUE),
        .groups="drop"
      )
    tabel(desk, "02_deskriptif")
  })

  bagian("2b. Validasi trajectory simulasi", function() {
    val <- dat |> group_by(treatment,time,week) |>
      summarise(mean_pct=mean(pct_change), sd_pct=sd(pct_change), .groups="drop")
    tabel(val, "02_validasi_trajectory")
    narasi("W30 dapat dibandingkan dengan LSM publikasi. Titik waktu antara merupakan trajectory simulasi yang mengikuti pola Figure 2a, bukan raw data asli.")
  })

  bagian("3. G1 Profil rerata berat badan dan 95% CI", function() {
    prof <- dat |> group_by(treatment,week) |>
      summarise(n=n(), mean=mean(weight_kg), sd=sd(weight_kg),
                se=sd/sqrt(n), lo=mean-qt(.975,n-1)*se,
                hi=mean+qt(.975,n-1)*se, .groups="drop")
    p <- ggplot(prof, aes(week,mean,group=treatment,colour=treatment)) +
      geom_line(linewidth=1) + geom_point(size=2.2) +
      geom_errorbar(aes(ymin=lo,ymax=hi), width=.4, alpha=.65) +
      scale_x_continuous(breaks=minggu_nilai) +
      theme_bw(base_size=12) + theme(legend.position="bottom") +
      labs(title="G1 Profil rerata berat badan sepanjang 30 minggu",
           x="Minggu", y="Berat badan (kg)", colour="Kelompok")
    gambar(p,"G1_profil_rerata_berat_badan",11,6.5)
  })

  bagian("4. G2 Lintasan individu", function() {
    p <- ggplot(dat, aes(week,weight_kg,group=id)) +
      geom_line(alpha=.15) +
      stat_summary(aes(group=treatment), fun=mean, geom="line",
                   linewidth=1.1, colour="black") +
      facet_wrap(~treatment, ncol=2, scales="free_y") +
      scale_x_continuous(breaks=minggu_nilai) +
      theme_bw(base_size=10) +
      labs(title="G2 Lintasan individu dan rerata tiap kelompok",
           x="Minggu", y="Berat badan (kg)")
    gambar(p,"G2_lintasan_individu",12,9)
  })

  bagian("5. G3 Korelasi antarwaktu", function() {
    ww <- wide |> select(all_of(level_waktu))
    cm <- cor(ww, use="pairwise.complete.obs")
    corr <- as.data.frame(as.table(cm))
    names(corr) <- c("waktu1","waktu2","r")
    p <- ggplot(corr, aes(waktu1,waktu2,fill=r)) +
      geom_tile() + geom_text(aes(label=sprintf("%.2f",r)), size=3) +
      theme_bw(base_size=11) +
      labs(title="G3 Korelasi berat badan antarwaktu",x=NULL,y=NULL,fill="r")
    gambar(p,"G3_korelasi_antarwaktu",8,7)
    tabel(cbind(waktu=rownames(cm), as.data.frame(cm)), "05_korelasi")
  })

  bagian("6. G4 Distribusi berat badan kelompok × waktu", function() {
    p <- ggplot(dat,aes(time,weight_kg,fill=treatment)) +
      geom_boxplot(outlier.alpha=.35) +
      facet_wrap(~treatment,ncol=2) +
      theme_bw(base_size=10) + theme(legend.position="none") +
      labs(title="G4 Distribusi berat badan menurut kelompok dan waktu",
           x="Waktu",y="Berat badan (kg)")
    gambar(p,"G4_boxplot_kelompok_waktu",12,8)
  })

  bagian("7. Normalitas dan outlier", function() {
    sh <- dat |> group_by(treatment,time) |>
      shapiro_test(weight_kg)
    tabel(sh,"07_shapiro_per_sel")

    out <- dat |> group_by(treatment,time) |> identify_outliers(weight_kg)
    tabel(out,"07_outlier")

    p <- ggplot(d1,aes(sample=weight_kg)) +
      stat_qq() + stat_qq_line() +
      facet_wrap(~time,scales="free") + theme_bw(base_size=10) +
      labs(title=paste("G5 Q-Q plot:",KELOMPOK_FOKUS))
    gambar(p,"G5_qq_kelompok_fokus",12,7)
  })

  bagian("8. Homogenitas varians antar kelompok", function() {
    lev <- dat |> group_by(time) |> levene_test(weight_kg ~ treatment)
    tabel(lev,"08_levene_per_waktu")
  })

  # ================================================================
  aov1 <- NULL
  bagian("9. RM ANOVA kelompok fokus", function() {
    aov1 <<- afex::aov_ez(
      id="id", dv="weight_kg", data=d1, within="time",
      anova_table=list(es=c("ges","pes"), correction="GG")
    )
    output(aov1); output(summary(aov1))
    tabel(as.data.frame(aov1$anova_table) |> rownames_to_column("efek"),
          "09_rm_anova_fokus")
    output(effectsize::eta_squared(aov1, partial=TRUE))

    e1 <- emmeans(aov1, ~ time)
    tabel(summary(pairs(e1, adjust="bonferroni"), infer=c(TRUE,TRUE)),
          "09_pairwise_bonferroni")

    fr <- d1 |> friedman_test(weight_kg ~ time | id)
    tabel(fr,"09_friedman")
    pw <- d1 |> pairwise_wilcox_test(weight_kg ~ time, paired=TRUE,
                                     p.adjust.method="bonferroni")
    tabel(pw,"09_wilcoxon_pairwise")
  })

  # ================================================================
  aov2 <- NULL
  bagian("10. Mixed Design ANOVA kelompok × waktu", function() {
    aov2 <<- afex::aov_ez(
      id="id", dv="weight_kg", data=cc,
      between="treatment", within="time",
      anova_table=list(es=c("ges","pes"), correction="GG")
    )
    output(aov2); output(summary(aov2))
    tabel(as.data.frame(aov2$anova_table) |> rownames_to_column("efek"),
          "10_mixed_anova")
    output(effectsize::eta_squared(aov2, partial=TRUE))
  })

  bagian("10b. G6 Visualisasi interaksi Mixed Design ANOVA", function() {
    if (is.null(aov2)) stop("Mixed ANOVA belum tersedia")
    p <- afex::afex_plot(
      aov2, x="time", trace="treatment",
      error="within", data_geom=geom_point,
      data_alpha=.12, legend_title="Kelompok"
    ) + theme_bw(base_size=11) +
      theme(legend.position="bottom") +
      labs(title="G6 Interaksi kelompok × waktu",
           x="Waktu",y="Berat badan (kg)")
    gambar(p,"G6_interaksi_mixed_anova",11,7)
  })

  bagian("10c. Simple effects dan post-hoc", function() {
    if (is.null(aov2)) stop("Mixed ANOVA belum tersedia")
    em <- emmeans(aov2, ~ treatment * time)

    tabel(summary(emmeans(aov2, ~ treatment | time)),
          "10_simple_kelompok_per_waktu")
    tabel(summary(pairs(emmeans(aov2, ~ treatment | time), adjust="tukey"),
                  infer=c(TRUE,TRUE)),
          "10_pairwise_kelompok_per_waktu")

    tabel(summary(emmeans(aov2, ~ time | treatment)),
          "10_simple_waktu_per_kelompok")
    tabel(summary(pairs(emmeans(aov2, ~ time | treatment), adjust="holm"),
                  infer=c(TRUE,TRUE)),
          "10_pairwise_waktu_per_kelompok")

    # Perbandingan perubahan W30-W0 antarkelompok
    ctime <- list("W30-W0"=c(-1,0,0,0,0,0,0,1))
    diffchg <- contrast(em,
      interaction=list(time=ctime, treatment="pairwise"),
      adjust="holm")
    tabel(summary(diffchg, infer=c(TRUE,TRUE)),
          "10_kontras_perubahan_W30_W0")
  })

  # ================================================================
  lmm <- NULL
  bagian("11. Linear Mixed Model", function() {
    # Waktu fixed effect kategorik; week numerik hanya untuk random slope.
    ctrl <- lmerControl(optimizer="bobyqa", optCtrl=list(maxfun=2e5))
    m0 <- lmer(weight_kg ~ treatment*time + (1|id),
               data=dat_obs, REML=TRUE, control=ctrl)
    m1 <- suppressWarnings(
      lmer(weight_kg ~ treatment*time + (1+week|id),
           data=dat_obs, REML=TRUE, control=ctrl)
    )

    if (!performance::check_singularity(m1)) {
      lmm <<- m1
      narasi("Model random intercept + random slope dipilih karena tidak singular.")
    } else {
      lmm <<- m0
      narasi("Model random slope singular; analisis menggunakan random intercept.")
    }

    output(summary(lmm))
    tab <- anova(lmm, type=3, ddf="Kenward-Roger")
    tabel(as.data.frame(tab) |> rownames_to_column("efek"),
          "11_lmm_type3")

    comp <- data.frame(
      model=c("Random intercept","Random intercept + slope"),
      AIC=c(AIC(m0),AIC(m1)),
      BIC=c(BIC(m0),BIC(m1)),
      singular=c(performance::check_singularity(m0),
                 performance::check_singularity(m1))
    )
    tabel(comp,"11_perbandingan_model")

    em <- emmeans(lmm, ~ treatment | time, lmer.df="kenward-roger")
    tabel(summary(em,infer=c(TRUE,TRUE)),"11_lmm_emmeans")
    tabel(summary(pairs(em,adjust="tukey"),infer=c(TRUE,TRUE)),
          "11_lmm_pairwise_kelompok")
  })

  bagian("12. Diagnostik LMM", function() {
    if (is.null(lmm)) stop("LMM belum tersedia")
    diag <- data.frame(
      fitted=fitted(lmm), residual=resid(lmm),
      treatment=dat_obs$treatment, time=dat_obs$time
    )

    p1 <- ggplot(diag,aes(sample=residual)) +
      stat_qq() + stat_qq_line() + theme_bw() +
      labs(title="G7 Q-Q residual LMM")
    gambar(p1,"G7_qq_residual_lmm",8,6)

    p2 <- ggplot(diag,aes(fitted,residual,colour=treatment)) +
      geom_point(alpha=.30) + geom_hline(yintercept=0,linetype=2) +
      theme_bw() + theme(legend.position="bottom") +
      labs(title="G8 Residual terhadap prediksi",
           x="Prediksi",y="Residual")
    gambar(p2,"G8_residual_vs_prediksi",9,6)

    p3 <- ggplot(diag,aes(time,residual,fill=treatment)) +
      geom_boxplot(outlier.alpha=.30) +
      facet_wrap(~treatment,ncol=2) + theme_bw() +
      theme(legend.position="none") +
      labs(title="G9 Residual menurut kelompok dan waktu")
    gambar(p3,"G9_residual_kelompok_waktu",12,8)

    re <- ranef(lmm)$id
    if ("(Intercept)" %in% names(re)) {
      p4 <- ggplot(data.frame(nilai=re[["(Intercept)"]]),aes(sample=nilai)) +
        stat_qq()+stat_qq_line()+theme_bw()+
        labs(title="G10 Q-Q random intercept")
      gambar(p4,"G10_qq_random_intercept",8,6)
    }
    if ("week" %in% names(re)) {
      p5 <- ggplot(data.frame(nilai=re[["week"]]),aes(sample=nilai)) +
        stat_qq()+stat_qq_line()+theme_bw()+
        labs(title="G11 Q-Q random slope minggu")
      gambar(p5,"G11_qq_random_slope",8,6)
    }

    output(performance::check_singularity(lmm))
    output(performance::check_convergence(lmm))
  })

  bagian("12b. G12 Diagnostik LMM tiga panel", function() {
    if (is.null(lmm)) stop("LMM belum tersedia")
    path <- file.path(folder,"G12_diagnostik_lmm_tiga_panel.png")
    png(path,width=2400,height=800,res=180)
    old <- par(mfrow=c(1,3),mar=c(4,4,3,1))
    qqnorm(resid(lmm),main="Q-Q residual"); qqline(resid(lmm))
    qqnorm(ranef(lmm)$id[,1],main="Q-Q random intercept")
    qqline(ranef(lmm)$id[,1])
    plot(fitted(lmm),resid(lmm),pch=16,col=adjustcolor("black",.25),
         xlab="Nilai prediksi",ylab="Residual",main="Residual vs prediksi")
    abline(h=0,lty=2)
    par(old); dev.off()
    tambah(paste0("<figure><img src='",
                  base64enc::dataURI(file=path,mime="image/png"),
                  "'><figcaption>G12 diagnostik LMM tiga panel</figcaption></figure>"))
  })

  bagian("13. G13 Rerata prediksi LMM dan 95% CI", function() {
    if (is.null(lmm)) stop("LMM belum tersedia")
    pred <- as.data.frame(summary(
      emmeans(lmm, ~ treatment | time, lmer.df="kenward-roger")
    ))
    pred$week <- minggu_nilai[as.integer(pred$time)]

    p <- ggplot(pred,aes(week,emmean,colour=treatment,group=treatment)) +
      geom_line(linewidth=1) + geom_point(size=2.2) +
      geom_errorbar(aes(ymin=lower.CL,ymax=upper.CL),width=.4) +
      scale_x_continuous(breaks=minggu_nilai) +
      theme_bw(base_size=11) + theme(legend.position="bottom") +
      labs(title="G13 Estimated marginal means LMM dan 95% CI",
           x="Minggu",y="Berat badan prediksi (kg)",colour="Kelompok")
    gambar(p,"G13_lmm_estimated_marginal_means",11,7)
  })

  bagian("14. Perubahan W30 dari baseline", function() {
    ch <- wide |>
      mutate(change_W30=W30-W0,
             pct_W30=(W30-W0)/W0*100)
    sumch <- ch |> group_by(treatment) |>
      summarise(n=n(), mean_change=mean(change_W30),
                sd_change=sd(change_W30),
                mean_pct=mean(pct_W30), sd_pct=sd(pct_W30),
                .groups="drop")
    tabel(sumch,"14_perubahan_W30")
  })

  bagian("15. Kesimpulan analitik otomatis", function() {
    narasi("Interpretasi utama harus memprioritaskan interaksi kelompok × waktu. Efek waktu saja hanya menyatakan bahwa berat badan berubah sepanjang follow-up, sedangkan interaksi menilai apakah pola perubahan berbeda menurut regimen.")
    narasi("Karena data individu merupakan simulasi, nilai p, effect size, residual, dan hasil post-hoc adalah hasil dataset sintetis. Nilai tersebut tidak boleh ditulis sebagai hasil asli uji klinis Ji et al. (2026).")
    narasi("Hasil sumber yang dapat disebut sebagai hasil artikel adalah struktur enam kelompok, jadwal kunjungan, dan endpoint week-30 yang dilaporkan pada publikasi.")
  })

  # ------------ status + session ------------
  write.csv(status,file.path(folder,"00_status_analisis.csv"),row.names=FALSE)
  capture.output(sessionInfo(), file=file.path(folder,"99_sessionInfo.txt"))

  # ------------ HTML ------------
  css <- "
  <style>
  body{font-family:Arial,sans-serif;max-width:1150px;margin:30px auto;line-height:1.5;color:#222}
  h1{font-size:24px} h2{font-size:18px;border-bottom:1px solid #ddd;padding-bottom:6px}
  table{border-collapse:collapse;margin:12px 0;font-size:12px}
  th,td{border:1px solid #ddd;padding:5px 7px}
  th{background:#f2f2f2}
  pre{background:#f7f7f7;padding:12px;white-space:pre-wrap;font-size:11px}
  figure{margin:24px 0} img{max-width:100%;height:auto}
  figcaption{font-size:12px;color:#555}
  </style>"
  header <- paste0(
    "<html><head><meta charset='utf-8'>",css,"</head><body>",
    "<h1>Analisis Repeated Measures Bofanglutide 2026</h1>",
    "<p><b>Catatan:</b> data individu adalah simulasi untuk pembelajaran, bukan raw data trial.</p>"
  )
  footer <- paste0(
    "<hr><p>Waktu proses: ", round(as.numeric(difftime(Sys.time(),waktu_mulai,units="secs")),1),
    " detik.</p></body></html>"
  )
  html <- paste(c(header,isi,footer),collapse="\n")
  html_path <- file.path(folder,"Laporan_Analisis_Bofanglutide_2026.html")
  writeLines(html,html_path,useBytes=TRUE)

  message("\nSELESAI. Folder hasil: ",folder)
  message("Laporan HTML: ",html_path)
  if (BUKA_HTML_OTOMATIS && interactive()) browseURL(html_path)
  invisible(list(folder=folder,status=status))
}

jalankan_analisis()
## Registered S3 method overwritten by 'lme4':
##   method           from
##   na.action.merMod car
## 
## >>> 1. Sumber, desain, dan integritas data
##                 treatment time  n
## 1                 Placebo   W0 66
## 2                 Placebo   W4 66
## 3                 Placebo   W8 66
## 4                 Placebo  W12 66
## 5                 Placebo  W16 66
## 6                 Placebo  W20 66
## 7                 Placebo  W24 66
## 8                 Placebo  W30 66
## 9  Bofanglutide 12 mg Q2W   W0 52
## 10 Bofanglutide 12 mg Q2W   W4 52
## 11 Bofanglutide 12 mg Q2W   W8 52
## 12 Bofanglutide 12 mg Q2W  W12 52
## 13 Bofanglutide 12 mg Q2W  W16 52
## 14 Bofanglutide 12 mg Q2W  W20 52
## 15 Bofanglutide 12 mg Q2W  W24 52
## 16 Bofanglutide 12 mg Q2W  W30 52
## 17 Bofanglutide 18 mg Q2W   W0 53
## 18 Bofanglutide 18 mg Q2W   W4 53
## 19 Bofanglutide 18 mg Q2W   W8 53
## 20 Bofanglutide 18 mg Q2W  W12 53
## 21 Bofanglutide 18 mg Q2W  W16 53
## 22 Bofanglutide 18 mg Q2W  W20 53
## 23 Bofanglutide 18 mg Q2W  W24 53
## 24 Bofanglutide 18 mg Q2W  W30 53
## 25 Bofanglutide 24 mg Q2W   W0 52
## 26 Bofanglutide 24 mg Q2W   W4 52
## 27 Bofanglutide 24 mg Q2W   W8 52
## 28 Bofanglutide 24 mg Q2W  W12 52
## 29 Bofanglutide 24 mg Q2W  W16 52
## 30 Bofanglutide 24 mg Q2W  W20 52
## 31 Bofanglutide 24 mg Q2W  W24 52
## 32 Bofanglutide 24 mg Q2W  W30 52
## 33 Bofanglutide 48 mg Q2W   W0 64
## 34 Bofanglutide 48 mg Q2W   W4 64
## 35 Bofanglutide 48 mg Q2W   W8 64
## 36 Bofanglutide 48 mg Q2W  W12 64
## 37 Bofanglutide 48 mg Q2W  W16 64
## 38 Bofanglutide 48 mg Q2W  W20 64
## 39 Bofanglutide 48 mg Q2W  W24 64
## 40 Bofanglutide 48 mg Q2W  W30 64
## 41  Bofanglutide 24 mg QW   W0 53
## 42  Bofanglutide 24 mg QW   W4 53
## 43  Bofanglutide 24 mg QW   W8 53
## 44  Bofanglutide 24 mg QW  W12 53
## 45  Bofanglutide 24 mg QW  W16 53
## 46  Bofanglutide 24 mg QW  W20 53
## 47  Bofanglutide 24 mg QW  W24 53
## 48  Bofanglutide 24 mg QW  W30 53
##                 treatment time  n n_tersedia n_hilang
## 1                 Placebo   W0 66         66        0
## 2                 Placebo   W4 66         66        0
## 3                 Placebo   W8 66         66        0
## 4                 Placebo  W12 66         66        0
## 5                 Placebo  W16 66         66        0
## 6                 Placebo  W20 66         66        0
## 7                 Placebo  W24 66         66        0
## 8                 Placebo  W30 66         66        0
## 9  Bofanglutide 12 mg Q2W   W0 52         52        0
## 10 Bofanglutide 12 mg Q2W   W4 52         52        0
## 11 Bofanglutide 12 mg Q2W   W8 52         52        0
## 12 Bofanglutide 12 mg Q2W  W12 52         52        0
## 13 Bofanglutide 12 mg Q2W  W16 52         52        0
## 14 Bofanglutide 12 mg Q2W  W20 52         52        0
## 15 Bofanglutide 12 mg Q2W  W24 52         52        0
## 16 Bofanglutide 12 mg Q2W  W30 52         52        0
## 17 Bofanglutide 18 mg Q2W   W0 53         53        0
## 18 Bofanglutide 18 mg Q2W   W4 53         53        0
## 19 Bofanglutide 18 mg Q2W   W8 53         53        0
## 20 Bofanglutide 18 mg Q2W  W12 53         53        0
## 21 Bofanglutide 18 mg Q2W  W16 53         53        0
## 22 Bofanglutide 18 mg Q2W  W20 53         53        0
## 23 Bofanglutide 18 mg Q2W  W24 53         53        0
## 24 Bofanglutide 18 mg Q2W  W30 53         53        0
## 25 Bofanglutide 24 mg Q2W   W0 52         52        0
## 26 Bofanglutide 24 mg Q2W   W4 52         52        0
## 27 Bofanglutide 24 mg Q2W   W8 52         52        0
## 28 Bofanglutide 24 mg Q2W  W12 52         52        0
## 29 Bofanglutide 24 mg Q2W  W16 52         52        0
## 30 Bofanglutide 24 mg Q2W  W20 52         52        0
## 31 Bofanglutide 24 mg Q2W  W24 52         52        0
## 32 Bofanglutide 24 mg Q2W  W30 52         52        0
## 33 Bofanglutide 48 mg Q2W   W0 64         64        0
## 34 Bofanglutide 48 mg Q2W   W4 64         64        0
## 35 Bofanglutide 48 mg Q2W   W8 64         64        0
## 36 Bofanglutide 48 mg Q2W  W12 64         64        0
## 37 Bofanglutide 48 mg Q2W  W16 64         64        0
## 38 Bofanglutide 48 mg Q2W  W20 64         64        0
## 39 Bofanglutide 48 mg Q2W  W24 64         64        0
## 40 Bofanglutide 48 mg Q2W  W30 64         64        0
## 41  Bofanglutide 24 mg QW   W0 53         53        0
## 42  Bofanglutide 24 mg QW   W4 53         53        0
## 43  Bofanglutide 24 mg QW   W8 53         53        0
## 44  Bofanglutide 24 mg QW  W12 53         53        0
## 45  Bofanglutide 24 mg QW  W16 53         53        0
## 46  Bofanglutide 24 mg QW  W20 53         53        0
## 47  Bofanglutide 24 mg QW  W24 53         53        0
## 48  Bofanglutide 24 mg QW  W30 53         53        0
## 
## >>> 2. Statistik deskriptif
##                 treatment time week  n     mean       sd  median    min     max
## 1                 Placebo   W0    0 66 95.60005 16.29755 95.6205 62.539 133.605
## 2                 Placebo   W4    4 66 95.03509 16.36649 95.6135 61.254 133.484
## 3                 Placebo   W8    8 66 94.43308 16.17452 94.8290 60.908 130.749
## 4                 Placebo  W12   12 66 94.15509 16.22767 95.4275 59.949 132.721
## 5                 Placebo  W16   16 66 94.34009 16.71254 94.9100 59.063 133.542
## 6                 Placebo  W20   20 66 94.32868 16.84561 93.4720 59.421 131.553
## 7                 Placebo  W24   24 66 94.53573 16.40236 97.2815 62.396 131.173
## 8                 Placebo  W30   30 66 94.64341 17.42706 95.6750 58.628 133.564
## 9  Bofanglutide 12 mg Q2W   W0    0 52 95.59994 15.94797 95.3110 65.450 123.970
## 10 Bofanglutide 12 mg Q2W   W4    4 52 93.35688 15.54533 93.7500 64.368 123.437
## 11 Bofanglutide 12 mg Q2W   W8    8 52 91.07038 15.17271 91.3295 63.591 117.707
## 12 Bofanglutide 12 mg Q2W  W12   12 52 89.74763 14.54793 90.5920 62.334 115.173
## 13 Bofanglutide 12 mg Q2W  W16   16 52 88.41281 15.16858 88.6935 60.191 115.035
## 14 Bofanglutide 12 mg Q2W  W20   20 52 87.35044 14.80497 87.3985 59.574 112.552
## 15 Bofanglutide 12 mg Q2W  W24   24 52 86.50338 14.28792 88.6465 57.271 113.543
## 16 Bofanglutide 12 mg Q2W  W30   30 52 86.24948 15.22747 87.3320 62.005 116.218
## 17 Bofanglutide 18 mg Q2W   W0    0 53 95.60002 18.08652 97.0410 62.974 136.241
## 18 Bofanglutide 18 mg Q2W   W4    4 53 92.94379 17.87625 91.4280 59.545 132.802
## 19 Bofanglutide 18 mg Q2W   W8    8 53 90.54409 17.37334 89.9380 58.318 128.533
## 20 Bofanglutide 18 mg Q2W  W12   12 53 89.16662 17.01255 88.4810 57.022 126.634
## 21 Bofanglutide 18 mg Q2W  W16   16 53 87.41823 17.10679 86.6940 54.147 119.885
## 22 Bofanglutide 18 mg Q2W  W20   20 53 85.94349 16.26878 83.6240 55.434 117.612
## 23 Bofanglutide 18 mg Q2W  W24   24 53 85.01032 16.55063 84.7390 52.897 117.131
## 24 Bofanglutide 18 mg Q2W  W30   30 53 83.63896 16.64148 81.6670 52.323 115.748
## 25 Bofanglutide 24 mg Q2W   W0    0 52 95.60002 14.08287 95.8860 62.058 123.561
## 26 Bofanglutide 24 mg Q2W   W4    4 52 92.63042 13.13257 93.7135 59.324 116.758
## 27 Bofanglutide 24 mg Q2W   W8    8 52 90.05506 12.74527 91.8540 56.766 111.924
## 28 Bofanglutide 24 mg Q2W  W12   12 52 88.57467 12.44859 90.3005 56.856 112.784
## 29 Bofanglutide 24 mg Q2W  W16   16 52 86.69915 12.55222 87.8240 56.341 110.556
## 30 Bofanglutide 24 mg Q2W  W20   20 52 85.04369 12.23788 86.7020 55.607 110.152
## 31 Bofanglutide 24 mg Q2W  W24   24 52 83.98846 12.05581 86.5465 52.913 108.800
## 32 Bofanglutide 24 mg Q2W  W30   30 52 82.36304 11.90554 84.3125 53.488 106.416
## 33 Bofanglutide 48 mg Q2W   W0    0 64 95.60002 17.78114 95.1815 64.191 148.418
## 34 Bofanglutide 48 mg Q2W   W4    4 64 92.43313 16.66829 92.9055 61.027 140.193
## 35 Bofanglutide 48 mg Q2W   W8    8 64 89.76919 16.41872 88.4350 61.842 135.021
## 36 Bofanglutide 48 mg Q2W  W12   12 64 87.99405 15.74837 87.2085 59.525 127.346
## 37 Bofanglutide 48 mg Q2W  W16   16 64 85.88577 15.57402 84.6650 59.692 129.047
## 38 Bofanglutide 48 mg Q2W  W20   20 64 83.87587 15.74840 82.6980 57.761 130.197
## 39 Bofanglutide 48 mg Q2W  W24   24 64 81.88464 14.53236 81.4650 56.931 121.564
## 40 Bofanglutide 48 mg Q2W  W30   30 64 80.02989 14.67692 77.8790 54.017 122.937
## 41  Bofanglutide 24 mg QW   W0    0 53 95.59994 17.50856 97.5550 65.288 143.257
## 42  Bofanglutide 24 mg QW   W4    4 53 92.00428 17.13206 92.2530 61.725 138.404
## 43  Bofanglutide 24 mg QW   W8    8 53 89.15998 16.67567 90.1930 60.045 132.797
## 44  Bofanglutide 24 mg QW  W12   12 53 87.06358 16.55994 86.6860 58.586 132.162
## 45  Bofanglutide 24 mg QW  W16   16 53 85.06030 16.19944 84.9020 56.647 124.003
## 46  Bofanglutide 24 mg QW  W20   20 53 83.11740 16.12978 83.2390 54.380 122.067
## 47  Bofanglutide 24 mg QW  W24   24 53 81.31128 15.46738 82.2850 54.269 114.079
## 48  Bofanglutide 24 mg QW  W30   30 53 79.75681 15.51428 80.4560 51.666 119.028
##       mean_pct   sd_pct
## 1    0.0000000 0.000000
## 2   -0.5999697 2.070872
## 3   -1.2000303 2.220168
## 4   -1.4999545 3.005713
## 5   -1.4000606 3.038797
## 6   -1.3999848 3.632361
## 7   -1.1000000 3.589256
## 8   -1.1499242 4.584249
## 9    0.0000000 0.000000
## 10  -2.3000000 2.082521
## 11  -4.7000192 2.294067
## 12  -5.9999808 2.906479
## 13  -7.5000385 3.558521
## 14  -8.5999808 3.119501
## 15  -9.4000385 3.818683
## 16  -9.7499615 5.024100
## 17   0.0000000 0.000000
## 18  -2.8000000 2.135692
## 19  -5.2999623 2.800287
## 20  -6.7000000 3.100216
## 21  -8.6000566 3.709293
## 22  -9.9998868 4.010521
## 23 -11.1000189 4.130342
## 24 -12.5500000 5.019624
## 25   0.0000000 0.000000
## 26  -3.0000192 2.107573
## 27  -5.7000192 2.296999
## 28  -7.2000192 3.309627
## 29  -9.2000577 3.550893
## 30 -10.8999808 4.009778
## 31 -12.0000000 4.250161
## 32 -13.6600192 4.984187
## 33   0.0000000 0.000000
## 34  -3.1999531 2.148943
## 35  -5.9999531 2.621524
## 36  -7.8000156 2.940092
## 37 -10.0000000 3.434752
## 38 -12.2000000 3.539502
## 39 -14.0999844 4.523457
## 40 -16.0899688 4.808888
## 41   0.0000000 0.000000
## 42  -3.8000000 1.886321
## 43  -6.7999434 2.325056
## 44  -9.0000000 3.177839
## 45 -11.1000000 3.215738
## 46 -13.1999811 3.329883
## 47 -15.0000189 4.127251
## 48 -16.6900189 4.129449
## 
## >>> 2b. Validasi trajectory simulasi
##                 treatment time week    mean_pct   sd_pct
## 1                 Placebo   W0    0   0.0000000 0.000000
## 2                 Placebo   W4    4  -0.5999697 2.070872
## 3                 Placebo   W8    8  -1.2000303 2.220168
## 4                 Placebo  W12   12  -1.4999545 3.005713
## 5                 Placebo  W16   16  -1.4000606 3.038797
## 6                 Placebo  W20   20  -1.3999848 3.632361
## 7                 Placebo  W24   24  -1.1000000 3.589256
## 8                 Placebo  W30   30  -1.1499242 4.584249
## 9  Bofanglutide 12 mg Q2W   W0    0   0.0000000 0.000000
## 10 Bofanglutide 12 mg Q2W   W4    4  -2.3000000 2.082521
## 11 Bofanglutide 12 mg Q2W   W8    8  -4.7000192 2.294067
## 12 Bofanglutide 12 mg Q2W  W12   12  -5.9999808 2.906479
## 13 Bofanglutide 12 mg Q2W  W16   16  -7.5000385 3.558521
## 14 Bofanglutide 12 mg Q2W  W20   20  -8.5999808 3.119501
## 15 Bofanglutide 12 mg Q2W  W24   24  -9.4000385 3.818683
## 16 Bofanglutide 12 mg Q2W  W30   30  -9.7499615 5.024100
## 17 Bofanglutide 18 mg Q2W   W0    0   0.0000000 0.000000
## 18 Bofanglutide 18 mg Q2W   W4    4  -2.8000000 2.135692
## 19 Bofanglutide 18 mg Q2W   W8    8  -5.2999623 2.800287
## 20 Bofanglutide 18 mg Q2W  W12   12  -6.7000000 3.100216
## 21 Bofanglutide 18 mg Q2W  W16   16  -8.6000566 3.709293
## 22 Bofanglutide 18 mg Q2W  W20   20  -9.9998868 4.010521
## 23 Bofanglutide 18 mg Q2W  W24   24 -11.1000189 4.130342
## 24 Bofanglutide 18 mg Q2W  W30   30 -12.5500000 5.019624
## 25 Bofanglutide 24 mg Q2W   W0    0   0.0000000 0.000000
## 26 Bofanglutide 24 mg Q2W   W4    4  -3.0000192 2.107573
## 27 Bofanglutide 24 mg Q2W   W8    8  -5.7000192 2.296999
## 28 Bofanglutide 24 mg Q2W  W12   12  -7.2000192 3.309627
## 29 Bofanglutide 24 mg Q2W  W16   16  -9.2000577 3.550893
## 30 Bofanglutide 24 mg Q2W  W20   20 -10.8999808 4.009778
## 31 Bofanglutide 24 mg Q2W  W24   24 -12.0000000 4.250161
## 32 Bofanglutide 24 mg Q2W  W30   30 -13.6600192 4.984187
## 33 Bofanglutide 48 mg Q2W   W0    0   0.0000000 0.000000
## 34 Bofanglutide 48 mg Q2W   W4    4  -3.1999531 2.148943
## 35 Bofanglutide 48 mg Q2W   W8    8  -5.9999531 2.621524
## 36 Bofanglutide 48 mg Q2W  W12   12  -7.8000156 2.940092
## 37 Bofanglutide 48 mg Q2W  W16   16 -10.0000000 3.434752
## 38 Bofanglutide 48 mg Q2W  W20   20 -12.2000000 3.539502
## 39 Bofanglutide 48 mg Q2W  W24   24 -14.0999844 4.523457
## 40 Bofanglutide 48 mg Q2W  W30   30 -16.0899688 4.808888
## 41  Bofanglutide 24 mg QW   W0    0   0.0000000 0.000000
## 42  Bofanglutide 24 mg QW   W4    4  -3.8000000 1.886321
## 43  Bofanglutide 24 mg QW   W8    8  -6.7999434 2.325056
## 44  Bofanglutide 24 mg QW  W12   12  -9.0000000 3.177839
## 45  Bofanglutide 24 mg QW  W16   16 -11.1000000 3.215738
## 46  Bofanglutide 24 mg QW  W20   20 -13.1999811 3.329883
## 47  Bofanglutide 24 mg QW  W24   24 -15.0000189 4.127251
## 48  Bofanglutide 24 mg QW  W30   30 -16.6900189 4.129449
## 
## >>> 3. G1 Profil rerata berat badan dan 95% CI

## Grafik tersimpan: C:/Users/masdi/OneDrive/Desktop/2611018042_DIMAS_Biostat/hasil_analisis_bofanglutide_2026/G1_profil_rerata_berat_badan.png
## 
## >>> 4. G2 Lintasan individu

## Grafik tersimpan: C:/Users/masdi/OneDrive/Desktop/2611018042_DIMAS_Biostat/hasil_analisis_bofanglutide_2026/G2_lintasan_individu.png
## 
## >>> 5. G3 Korelasi antarwaktu

## Grafik tersimpan: C:/Users/masdi/OneDrive/Desktop/2611018042_DIMAS_Biostat/hasil_analisis_bofanglutide_2026/G3_korelasi_antarwaktu.png
##     waktu        W0        W4        W8       W12       W16       W20       W24
## W0     W0 1.0000000 0.9903542 0.9821573 0.9697463 0.9562284 0.9424946 0.9201176
## W4     W4 0.9903542 1.0000000 0.9921575 0.9858528 0.9776368 0.9674997 0.9512399
## W8     W8 0.9821573 0.9921575 1.0000000 0.9896071 0.9833011 0.9766764 0.9633456
## W12   W12 0.9697463 0.9858528 0.9896071 1.0000000 0.9855963 0.9800257 0.9698493
## W16   W16 0.9562284 0.9776368 0.9833011 0.9855963 1.0000000 0.9841911 0.9763677
## W20   W20 0.9424946 0.9674997 0.9766764 0.9800257 0.9841911 1.0000000 0.9774602
## W24   W24 0.9201176 0.9512399 0.9633456 0.9698493 0.9763677 0.9774602 1.0000000
## W30   W30 0.8991787 0.9351122 0.9494696 0.9588368 0.9668664 0.9719486 0.9753610
##           W30
## W0  0.8991787
## W4  0.9351122
## W8  0.9494696
## W12 0.9588368
## W16 0.9668664
## W20 0.9719486
## W24 0.9753610
## W30 1.0000000
## 
## >>> 6. G4 Distribusi berat badan kelompok × waktu

## Grafik tersimpan: C:/Users/masdi/OneDrive/Desktop/2611018042_DIMAS_Biostat/hasil_analisis_bofanglutide_2026/G4_boxplot_kelompok_waktu.png
## 
## >>> 7. Normalitas dan outlier
##                 treatment time  variable statistic         p
## 1                 Placebo   W0 weight_kg 0.9891159 0.8343321
## 2                 Placebo   W4 weight_kg 0.9887882 0.8178338
## 3                 Placebo   W8 weight_kg 0.9830754 0.5058327
## 4                 Placebo  W12 weight_kg 0.9861431 0.6729420
## 5                 Placebo  W16 weight_kg 0.9872351 0.7343276
## 6                 Placebo  W20 weight_kg 0.9839519 0.5516676
## 7                 Placebo  W24 weight_kg 0.9798946 0.3599343
## 8                 Placebo  W30 weight_kg 0.9842498 0.5676853
## 9  Bofanglutide 12 mg Q2W   W0 weight_kg 0.9706806 0.2257653
## 10 Bofanglutide 12 mg Q2W   W4 weight_kg 0.9769322 0.4050672
## 11 Bofanglutide 12 mg Q2W   W8 weight_kg 0.9682441 0.1778724
## 12 Bofanglutide 12 mg Q2W  W12 weight_kg 0.9682746 0.1784066
## 13 Bofanglutide 12 mg Q2W  W16 weight_kg 0.9695293 0.2018041
## 14 Bofanglutide 12 mg Q2W  W20 weight_kg 0.9692818 0.1969705
## 15 Bofanglutide 12 mg Q2W  W24 weight_kg 0.9794560 0.5029588
## 16 Bofanglutide 12 mg Q2W  W30 weight_kg 0.9656844 0.1380298
## 17 Bofanglutide 18 mg Q2W   W0 weight_kg 0.9780059 0.4322013
## 18 Bofanglutide 18 mg Q2W   W4 weight_kg 0.9831321 0.6543469
## 19 Bofanglutide 18 mg Q2W   W8 weight_kg 0.9771818 0.4014922
## 20 Bofanglutide 18 mg Q2W  W12 weight_kg 0.9818972 0.5971166
## 21 Bofanglutide 18 mg Q2W  W16 weight_kg 0.9792392 0.4811453
## 22 Bofanglutide 18 mg Q2W  W20 weight_kg 0.9744847 0.3125742
## 23 Bofanglutide 18 mg Q2W  W24 weight_kg 0.9755236 0.3447229
## 24 Bofanglutide 18 mg Q2W  W30 weight_kg 0.9714072 0.2320061
## 25 Bofanglutide 24 mg Q2W   W0 weight_kg 0.9836122 0.6881380
## 26 Bofanglutide 24 mg Q2W   W4 weight_kg 0.9774602 0.4243574
## 27 Bofanglutide 24 mg Q2W   W8 weight_kg 0.9707061 0.2263251
## 28 Bofanglutide 24 mg Q2W  W12 weight_kg 0.9698832 0.2089048
## 29 Bofanglutide 24 mg Q2W  W16 weight_kg 0.9785358 0.4656361
## 30 Bofanglutide 24 mg Q2W  W20 weight_kg 0.9804774 0.5463422
## 31 Bofanglutide 24 mg Q2W  W24 weight_kg 0.9793879 0.5001405
## 32 Bofanglutide 24 mg Q2W  W30 weight_kg 0.9844169 0.7251495
## 33 Bofanglutide 48 mg Q2W   W0 weight_kg 0.9724266 0.1616567
## 34 Bofanglutide 48 mg Q2W   W4 weight_kg 0.9795433 0.3662587
## 35 Bofanglutide 48 mg Q2W   W8 weight_kg 0.9707997 0.1331340
## 36 Bofanglutide 48 mg Q2W  W12 weight_kg 0.9761387 0.2498887
## 37 Bofanglutide 48 mg Q2W  W16 weight_kg 0.9730048 0.1731430
## 38 Bofanglutide 48 mg Q2W  W20 weight_kg 0.9685054 0.1010945
## 39 Bofanglutide 48 mg Q2W  W24 weight_kg 0.9768633 0.2715536
## 40 Bofanglutide 48 mg Q2W  W30 weight_kg 0.9763559 0.2562138
## 41  Bofanglutide 24 mg QW   W0 weight_kg 0.9702029 0.2059548
## 42  Bofanglutide 24 mg QW   W4 weight_kg 0.9762499 0.3687577
## 43  Bofanglutide 24 mg QW   W8 weight_kg 0.9781945 0.4394592
## 44  Bofanglutide 24 mg QW  W12 weight_kg 0.9715876 0.2361595
## 45  Bofanglutide 24 mg QW  W16 weight_kg 0.9765059 0.3775373
## 46  Bofanglutide 24 mg QW  W20 weight_kg 0.9753182 0.3381589
## 47  Bofanglutide 24 mg QW  W24 weight_kg 0.9764161 0.3744384
## 48  Bofanglutide 24 mg QW  W30 weight_kg 0.9771408 0.4000045
##                 treatment time     id week baseline_weight_kg pct_change
## 1                 Placebo  W12 BFG340   12            133.605     -0.662
## 2                 Placebo  W16 BFG340   16            133.605     -0.047
## 3  Bofanglutide 24 mg Q2W   W0 BFG138    0             62.058      0.000
## 4  Bofanglutide 24 mg Q2W   W0 BFG143    0             62.058      0.000
## 5  Bofanglutide 48 mg Q2W   W0 BFG176    0            148.418      0.000
## 6  Bofanglutide 48 mg Q2W   W4 BFG176    4            148.418     -5.542
## 7  Bofanglutide 48 mg Q2W   W8 BFG176    8            148.418     -9.027
## 8  Bofanglutide 48 mg Q2W  W12 BFG176   12            148.418    -14.198
## 9  Bofanglutide 48 mg Q2W  W16 BFG176   16            148.418    -13.052
## 10 Bofanglutide 48 mg Q2W  W20 BFG176   20            148.418    -12.277
## 11 Bofanglutide 48 mg Q2W  W24 BFG176   24            148.418    -18.094
## 12 Bofanglutide 48 mg Q2W  W30 BFG176   30            148.418    -17.169
## 13  Bofanglutide 24 mg QW   W0 BFG234    0            143.257      0.000
## 14  Bofanglutide 24 mg QW   W4 BFG234    4            143.257     -3.388
## 15  Bofanglutide 24 mg QW  W12 BFG234   12            143.257     -7.745
## 16  Bofanglutide 24 mg QW  W30 BFG234   30            143.257    -16.913
##    weight_kg is.outlier is.extreme
## 1    132.721       TRUE      FALSE
## 2    133.542       TRUE      FALSE
## 3     62.058       TRUE      FALSE
## 4     62.058       TRUE      FALSE
## 5    148.418       TRUE      FALSE
## 6    140.193       TRUE      FALSE
## 7    135.021       TRUE      FALSE
## 8    127.346       TRUE      FALSE
## 9    129.047       TRUE      FALSE
## 10   130.197       TRUE      FALSE
## 11   121.564       TRUE      FALSE
## 12   122.937       TRUE      FALSE
## 13   143.257       TRUE      FALSE
## 14   138.404       TRUE      FALSE
## 15   132.162       TRUE      FALSE
## 16   119.028       TRUE      FALSE

## Grafik tersimpan: C:/Users/masdi/OneDrive/Desktop/2611018042_DIMAS_Biostat/hasil_analisis_bofanglutide_2026/G5_qq_kelompok_fokus.png
## 
## >>> 8. Homogenitas varians antar kelompok
##   time df1 df2 statistic         p
## 1   W0   5 334 0.9061904 0.4771255
## 2   W4   5 334 1.2126824 0.3027384
## 3   W8   5 334 1.0361465 0.3961973
## 4  W12   5 334 1.2143288 0.3019560
## 5  W16   5 334 1.1058572 0.3570071
## 6  W20   5 334 1.1294070 0.3444446
## 7  W24   5 334 1.3260995 0.2525502
## 8  W30   5 334 1.3957018 0.2253290
## 
## >>> 9. RM ANOVA kelompok fokus
## Anova Table (Type 3 tests)
## 
## Response: weight_kg
##   Effect           df  MSE          F  ges  pes p.value
## 1   time 4.05, 210.76 8.82 305.90 *** .089 .855   <.001
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '+' 0.1 ' ' 1
## 
## Sphericity correction method: GG  
## 
## Univariate Type III Repeated-Measures ANOVA Assuming Sphericity
## 
##              Sum Sq num Df Error SS den Df F value    Pr(>F)    
## (Intercept) 3182325      1   110196     52  1501.7 < 2.2e-16 ***
## time          10932      7     1858    364   305.9 < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## 
## Mauchly Tests for Sphericity
## 
##      Test statistic    p-value
## time       0.085055 5.8987e-14
## 
## 
## Greenhouse-Geisser and Huynh-Feldt Corrections
##  for Departure from Sphericity
## 
##      GG eps Pr(>F[GG])    
## time  0.579  < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
##         HF eps   Pr(>F[HF])
## time 0.6340656 5.167047e-95 
##   efek   num Df   den Df      MSE        F       pes        ges       Pr(>F)
## 1 time 4.053031 210.7576 8.817394 305.8987 0.8547075 0.08888757 5.303935e-87
## # Effect Size for ANOVA (Type III)
## 
## Parameter | Eta2 (partial) |       95% CI
## -----------------------------------------
## time      |           0.85 | [0.84, 1.00]
## 
## - One-sided CIs: upper bound fixed at [1.00]. 
##  contrast   estimate        SE df  lower.CL  upper.CL t.ratio p.value
##  W0 - W4    3.595660 0.2612705 52  2.735239  4.456081  13.762 <0.0001
##  W0 - W8    6.439962 0.3421168 52  5.313297  7.566628  18.824 <0.0001
##  W0 - W12   8.536358 0.4640678 52  7.008081 10.064636  18.395 <0.0001
##  W0 - W16  10.539642 0.4814928 52  8.953980 12.125303  21.890 <0.0001
##  W0 - W20  12.482547 0.4762227 52 10.914242 14.050853  26.212 <0.0001
##  W0 - W24  14.288660 0.6628854 52 12.105634 16.471687  21.555 <0.0001
##  W0 - W30  15.843132 0.6235394 52 13.789680 17.896584  25.408 <0.0001
##  W4 - W8    2.844302 0.2472814 52  2.029950  3.658654  11.502 <0.0001
##  W4 - W12   4.940698 0.3272174 52  3.863100  6.018297  15.099 <0.0001
##  W4 - W16   6.943981 0.3482164 52  5.797228  8.090734  19.942 <0.0001
##  W4 - W20   8.886887 0.3510435 52  7.730824 10.042950  25.316 <0.0001
##  W4 - W24  10.693000 0.5379472 52  8.921422 12.464578  19.877 <0.0001
##  W4 - W30  12.247472 0.4915110 52 10.628818 13.866125  24.918 <0.0001
##  W8 - W12   2.096396 0.3110107 52  1.072170  3.120623   6.741 <0.0001
##  W8 - W16   4.099679 0.3715564 52  2.876063  5.323296  11.034 <0.0001
##  W8 - W20   6.042585 0.3402483 52  4.922073  7.163097  17.759 <0.0001
##  W8 - W24   7.848698 0.4901047 52  6.234676  9.462720  16.014 <0.0001
##  W8 - W30   9.403170 0.4214876 52  8.015119 10.791221  22.309 <0.0001
##  W12 - W16  2.003283 0.3958360 52  0.699708  3.306858   5.061  0.0002
##  W12 - W20  3.946189 0.3757691 52  2.708699  5.183679  10.502 <0.0001
##  W12 - W24  5.752302 0.5163315 52  4.051909  7.452695  11.141 <0.0001
##  W12 - W30  7.306774 0.4327228 52  5.881723  8.731824  16.886 <0.0001
##  W16 - W20  1.942906 0.3474007 52  0.798839  3.086972   5.593 <0.0001
##  W16 - W24  3.749019 0.4819767 52  2.161764  5.336274   7.778 <0.0001
##  W16 - W30  5.303491 0.4658178 52  3.769450  6.837531  11.385 <0.0001
##  W20 - W24  1.806113 0.4609174 52  0.288211  3.324015   3.919  0.0073
##  W20 - W30  3.360585 0.4304118 52  1.943144  4.778025   7.808 <0.0001
##  W24 - W30  1.554472 0.5246831 52 -0.173425  3.282368   2.963  0.1285
## 
## Confidence level used: 0.95 
## Conf-level adjustment: bonferroni method for 28 estimates 
## P value adjustment: bonferroni method for 28 tests
## Bagian gagal: not an unreplicated complete block design
## 
## >>> 10. Mixed Design ANOVA kelompok × waktu
## Anova Table (Type 3 tests)
## 
## Response: weight_kg
##           Effect             df     MSE          F  ges  pes p.value
## 1      treatment         5, 334 1947.84     2.17 + .031 .031    .057
## 2           time  3.85, 1286.79    9.56 939.90 *** .050 .738   <.001
## 3 treatment:time 19.26, 1286.79    9.56  43.53 *** .012 .395   <.001
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '+' 0.1 ' ' 1
## 
## Sphericity correction method: GG  
## 
## Univariate Type III Repeated-Measures ANOVA Assuming Sphericity
## 
##                  Sum Sq num Df Error SS den Df    F value  Pr(>F)    
## (Intercept)    21416509      1   650580    334 10994.9794 < 2e-16 ***
## treatment         21116      5   650580    334     2.1681 0.05731 .  
## time              34622      7    12303   2338   939.9018 < 2e-16 ***
## treatment:time     8018     35    12303   2338    43.5308 < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## 
## Mauchly Tests for Sphericity
## 
##                Test statistic     p-value
## time                  0.09846 1.9874e-144
## treatment:time        0.09846 1.9874e-144
## 
## 
## Greenhouse-Geisser and Huynh-Feldt Corrections
##  for Departure from Sphericity
## 
##                 GG eps Pr(>F[GG])    
## time           0.55038  < 2.2e-16 ***
## treatment:time 0.55038  < 2.2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
##                   HF eps    Pr(>F[HF])
## time           0.5576031  0.000000e+00
## treatment:time 0.5576031 5.930385e-127 
##             efek    num Df   den Df         MSE          F        pes
## 1      treatment  5.000000  334.000 1947.844347   2.168099 0.03143627
## 2           time  3.852652 1286.786    9.561255 939.901792 0.73781338
## 3 treatment:time 19.263258 1286.786    9.561255  43.530810 0.39454809
##          ges        Pr(>F)
## 1 0.03087081  5.731106e-02
## 2 0.04963744  0.000000e+00
## 3 0.01195042 2.405936e-125
## # Effect Size for ANOVA (Type III)
## 
## Parameter      | Eta2 (partial) |       95% CI
## ----------------------------------------------
## treatment      |           0.03 | [0.00, 1.00]
## time           |           0.74 | [0.73, 1.00]
## treatment:time |           0.39 | [0.37, 1.00]
## 
## - One-sided CIs: upper bound fixed at [1.00].
## 
## >>> 10b. G6 Visualisasi interaksi Mixed Design ANOVA

## Grafik tersimpan: C:/Users/masdi/OneDrive/Desktop/2611018042_DIMAS_Biostat/hasil_analisis_bofanglutide_2026/G6_interaksi_mixed_anova.png
## 
## >>> 10c. Simple effects dan post-hoc
## time = W0:
##  treatment                emmean       SE  df lower.CL  upper.CL
##  Placebo                95.60005 2.056410 334 91.55490  99.64519
##  Bofanglutide 12 mg Q2W 95.59994 2.316754 334 91.04267 100.15721
##  Bofanglutide 18 mg Q2W 95.60002 2.294794 334 91.08595 100.11409
##  Bofanglutide 24 mg Q2W 95.60002 2.316754 334 91.04275 100.15729
##  Bofanglutide 48 mg Q2W 95.60002 2.088294 334 91.49215  99.70788
##  Bofanglutide 24 mg QW  95.59994 2.294794 334 91.08587 100.11401
## 
## time = W4:
##  treatment                emmean       SE  df lower.CL  upper.CL
##  Placebo                95.03509 1.997029 334 91.10675  98.96343
##  Bofanglutide 12 mg Q2W 93.35688 2.249855 334 88.93121  97.78256
##  Bofanglutide 18 mg Q2W 92.94379 2.228529 334 88.56007  97.32751
##  Bofanglutide 24 mg Q2W 92.63042 2.249855 334 88.20475  97.05610
##  Bofanglutide 48 mg Q2W 92.43312 2.027992 334 88.44388  96.42237
##  Bofanglutide 24 mg QW  92.00428 2.228529 334 87.62056  96.38800
## 
## time = W8:
##  treatment                emmean       SE  df lower.CL  upper.CL
##  Placebo                94.43308 1.954027 334 90.58933  98.27683
##  Bofanglutide 12 mg Q2W 91.07038 2.201409 334 86.74001  95.40076
##  Bofanglutide 18 mg Q2W 90.54409 2.180543 334 86.25477  94.83342
##  Bofanglutide 24 mg Q2W 90.05506 2.201409 334 85.72468  94.38543
##  Bofanglutide 48 mg Q2W 89.76919 1.984324 334 85.86584  93.67253
##  Bofanglutide 24 mg QW  89.15998 2.180543 334 84.87065  93.44931
## 
## time = W12:
##  treatment                emmean       SE  df lower.CL  upper.CL
##  Placebo                94.15509 1.913918 334 90.39024  97.91994
##  Bofanglutide 12 mg Q2W 89.74763 2.156223 334 85.50615  93.98912
##  Bofanglutide 18 mg Q2W 89.16662 2.135784 334 84.96534  93.36791
##  Bofanglutide 24 mg Q2W 88.57467 2.156223 334 84.33319  92.81616
##  Bofanglutide 48 mg Q2W 87.99405 1.943593 334 84.17082  91.81727
##  Bofanglutide 24 mg QW  87.06358 2.135784 334 82.86230  91.26487
## 
## time = W16:
##  treatment                emmean       SE  df lower.CL  upper.CL
##  Placebo                94.34009 1.929497 334 90.54459  98.13559
##  Bofanglutide 12 mg Q2W 88.41281 2.173774 334 84.13679  92.68882
##  Bofanglutide 18 mg Q2W 87.41823 2.153169 334 83.18274  91.65371
##  Bofanglutide 24 mg Q2W 86.69915 2.173774 334 82.42314  90.97517
##  Bofanglutide 48 mg Q2W 85.88577 1.959414 334 82.03142  89.74011
##  Bofanglutide 24 mg QW  85.06030 2.153169 334 80.82482  89.29578
## 
## time = W20:
##  treatment                emmean       SE  df lower.CL  upper.CL
##  Placebo                94.32868 1.907146 334 90.57715  98.08021
##  Bofanglutide 12 mg Q2W 87.35044 2.148593 334 83.12396  91.57692
##  Bofanglutide 18 mg Q2W 85.94349 2.128227 334 81.75707  90.12991
##  Bofanglutide 24 mg Q2W 85.04369 2.148593 334 80.81721  89.27017
##  Bofanglutide 48 mg Q2W 83.87588 1.936716 334 80.06618  87.68557
##  Bofanglutide 24 mg QW  83.11740 2.128227 334 78.93098  87.30381
## 
## time = W24:
##  treatment                emmean       SE  df lower.CL  upper.CL
##  Placebo                94.53573 1.848223 334 90.90010  98.17135
##  Bofanglutide 12 mg Q2W 86.50338 2.082210 334 82.40749  90.59928
##  Bofanglutide 18 mg Q2W 85.01032 2.062473 334 80.95325  89.06740
##  Bofanglutide 24 mg Q2W 83.98846 2.082210 334 79.89256  88.08436
##  Bofanglutide 48 mg Q2W 81.88464 1.876879 334 78.19265  85.57663
##  Bofanglutide 24 mg QW  81.31128 2.062473 334 77.25421  85.36836
## 
## time = W30:
##  treatment                emmean       SE  df lower.CL  upper.CL
##  Placebo                94.64341 1.896465 334 90.91289  98.37393
##  Bofanglutide 12 mg Q2W 86.24948 2.136560 334 82.04667  90.45229
##  Bofanglutide 18 mg Q2W 83.63896 2.116308 334 79.47599  87.80193
##  Bofanglutide 24 mg Q2W 82.36304 2.136560 334 78.16023  86.56585
##  Bofanglutide 48 mg Q2W 80.02989 1.925869 334 76.24153  83.81825
##  Bofanglutide 24 mg QW  79.75681 2.116308 334 75.59384  83.91978
## 
## Confidence level used: 0.95 
## time = W0:
##  contrast                                         estimate       SE  df
##  Placebo - Bofanglutide 12 mg Q2W                 0.000103 3.097769 334
##  Placebo - Bofanglutide 18 mg Q2W                 0.000027 3.081380 334
##  Placebo - Bofanglutide 24 mg Q2W                 0.000026 3.097769 334
##  Placebo - Bofanglutide 48 mg Q2W                 0.000030 2.930835 334
##  Placebo - Bofanglutide 24 mg QW                  0.000102 3.081380 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 18 mg Q2W -0.000077 3.260894 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg Q2W -0.000077 3.276386 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 48 mg Q2W -0.000073 3.119026 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg QW  -0.000001 3.260894 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg Q2W  0.000000 3.260894 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 48 mg Q2W  0.000003 3.102749 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg QW   0.000075 3.245329 334
##  Bofanglutide 24 mg Q2W - Bofanglutide 48 mg Q2W  0.000004 3.119026 334
##  Bofanglutide 24 mg Q2W - Bofanglutide 24 mg QW   0.000076 3.260894 334
##  Bofanglutide 48 mg Q2W - Bofanglutide 24 mg QW   0.000072 3.102749 334
##   lower.CL  upper.CL t.ratio p.value
##  -8.879165  8.879371   0.000  1.0000
##  -8.832264  8.832318   0.000  1.0000
##  -8.879242  8.879294   0.000  1.0000
##  -8.400748  8.400808   0.000  1.0000
##  -8.832189  8.832393   0.000  1.0000
##  -9.346917  9.346764   0.000  1.0000
##  -9.391321  9.391167   0.000  1.0000
##  -8.940271  8.940124   0.000  1.0000
##  -9.346842  9.346839   0.000  1.0000
##  -9.346841  9.346840   0.000  1.0000
##  -8.893540  8.893546   0.000  1.0000
##  -9.302150  9.302301   0.000  1.0000
##  -8.940194  8.940201   0.000  1.0000
##  -9.346765  9.346916   0.000  1.0000
##  -8.893471  8.893615   0.000  1.0000
## 
## time = W4:
##  contrast                                         estimate       SE  df
##  Placebo - Bofanglutide 12 mg Q2W                 1.678206 3.008317 334
##  Placebo - Bofanglutide 18 mg Q2W                 2.091298 2.992401 334
##  Placebo - Bofanglutide 24 mg Q2W                 2.404668 3.008317 334
##  Placebo - Bofanglutide 48 mg Q2W                 2.601966 2.846203 334
##  Placebo - Bofanglutide 24 mg QW                  3.030808 2.992401 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 18 mg Q2W  0.413092 3.166732 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg Q2W  0.726462 3.181776 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 48 mg Q2W  0.923760 3.028960 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg QW   1.352602 3.166732 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg Q2W  0.313369 3.166732 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 48 mg Q2W  0.510667 3.013153 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg QW   0.939509 3.151616 334
##  Bofanglutide 24 mg Q2W - Bofanglutide 48 mg Q2W  0.197298 3.028960 334
##  Bofanglutide 24 mg Q2W - Bofanglutide 24 mg QW   0.626140 3.166732 334
##  Bofanglutide 48 mg Q2W - Bofanglutide 24 mg QW   0.428842 3.013153 334
##   lower.CL  upper.CL t.ratio p.value
##  -6.944661 10.301074   0.558  0.9936
##  -6.485949 10.668546   0.699  0.9820
##  -6.218200 11.027535   0.799  0.9675
##  -5.556229 10.760161   0.914  0.9427
##  -5.546439 11.608055   1.013  0.9135
##  -8.663846  9.490031   0.130  1.0000
##  -8.393598  9.846521   0.228  0.9999
##  -7.758278  9.605798   0.305  0.9996
##  -7.724337 10.429540   0.427  0.9982
##  -8.763569  9.390308   0.099  1.0000
##  -8.126063  9.147398   0.169  1.0000
##  -8.094102  9.973121   0.298  0.9997
##  -8.484740  8.879336   0.065  1.0000
##  -8.450798  9.703078   0.198  1.0000
##  -8.207888  9.065572   0.142  1.0000
## 
## time = W8:
##  contrast                                         estimate       SE  df
##  Placebo - Bofanglutide 12 mg Q2W                 3.362691 2.943539 334
##  Placebo - Bofanglutide 18 mg Q2W                 3.888981 2.927966 334
##  Placebo - Bofanglutide 24 mg Q2W                 4.378018 2.943539 334
##  Placebo - Bofanglutide 48 mg Q2W                 4.663888 2.784917 334
##  Placebo - Bofanglutide 24 mg QW                  5.273095 2.927966 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 18 mg Q2W  0.526290 3.098543 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg Q2W  1.015327 3.113263 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 48 mg Q2W  1.301197 2.963738 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg QW   1.910403 3.098543 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg Q2W  0.489037 3.098543 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 48 mg Q2W  0.774907 2.948272 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg QW   1.384113 3.083753 334
##  Bofanglutide 24 mg Q2W - Bofanglutide 48 mg Q2W  0.285870 2.963738 334
##  Bofanglutide 24 mg Q2W - Bofanglutide 24 mg QW   0.895077 3.098543 334
##  Bofanglutide 48 mg Q2W - Bofanglutide 24 mg QW   0.609206 2.948272 334
##   lower.CL  upper.CL t.ratio p.value
##  -5.074502 11.799884   1.142  0.8633
##  -4.503574 12.281536   1.328  0.7692
##  -4.059175 12.815211   1.487  0.6726
##  -3.318638 12.646414   1.675  0.5495
##  -3.119460 13.665650   1.801  0.4665
##  -8.355196  9.407777   0.170  1.0000
##  -7.908352  9.939006   0.326  0.9995
##  -7.193892  9.796286   0.439  0.9979
##  -6.971083 10.791890   0.617  0.9898
##  -8.392450  9.370523   0.158  1.0000
##  -7.675850  9.225664   0.263  0.9998
##  -7.454979 10.223206   0.449  0.9977
##  -8.209219  8.780960   0.096  1.0000
##  -7.986410  9.776563   0.289  0.9997
##  -7.841551  9.059964   0.207  0.9999
## 
## time = W12:
##  contrast                                         estimate       SE  df
##  Placebo - Bofanglutide 12 mg Q2W                 4.407456 2.883119 334
##  Placebo - Bofanglutide 18 mg Q2W                 4.988468 2.867866 334
##  Placebo - Bofanglutide 24 mg Q2W                 5.580418 2.883119 334
##  Placebo - Bofanglutide 48 mg Q2W                 6.161044 2.727752 334
##  Placebo - Bofanglutide 24 mg QW                  7.091506 2.867866 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 18 mg Q2W  0.581012 3.034941 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg Q2W  1.172962 3.049359 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 48 mg Q2W  1.753588 2.902903 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg QW   2.684050 3.034941 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg Q2W  0.591950 3.034941 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 48 mg Q2W  1.172576 2.887754 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg QW   2.103038 3.020455 334
##  Bofanglutide 24 mg Q2W - Bofanglutide 48 mg Q2W  0.580626 2.902903 334
##  Bofanglutide 24 mg Q2W - Bofanglutide 24 mg QW   1.511088 3.034941 334
##  Bofanglutide 48 mg Q2W - Bofanglutide 24 mg QW   0.930462 2.887754 334
##   lower.CL  upper.CL t.ratio p.value
##  -3.856552 12.671464   1.529  0.6460
##  -3.231818 13.208754   1.739  0.5067
##  -2.683590 13.844426   1.936  0.3823
##  -1.657629 13.979718   2.259  0.2142
##  -1.128780 15.311792   2.473  0.1352
##  -8.118170  9.280194   0.191  1.0000
##  -7.567547  9.913470   0.385  0.9989
##  -6.567128 10.074304   0.604  0.9907
##  -6.015132 11.383231   0.884  0.9501
##  -8.107232  9.291131   0.195  1.0000
##  -7.104718  9.449869   0.406  0.9986
##  -6.554620 10.760696   0.696  0.9823
##  -7.740090  8.901342   0.200  1.0000
##  -7.188093 10.210270   0.498  0.9962
##  -7.346832  9.207756   0.322  0.9995
## 
## time = W16:
##  contrast                                         estimate       SE  df
##  Placebo - Bofanglutide 12 mg Q2W                 5.927283 2.906588 334
##  Placebo - Bofanglutide 18 mg Q2W                 6.921864 2.891210 334
##  Placebo - Bofanglutide 24 mg Q2W                 7.640937 2.906588 334
##  Placebo - Bofanglutide 48 mg Q2W                 8.454325 2.749956 334
##  Placebo - Bofanglutide 24 mg QW                  9.279789 2.891210 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 18 mg Q2W  0.994581 3.059646 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg Q2W  1.713654 3.074181 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 48 mg Q2W  2.527042 2.926533 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg QW   3.352506 3.059646 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg Q2W  0.719073 3.059646 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 48 mg Q2W  1.532461 2.911261 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg QW   2.357925 3.045041 334
##  Bofanglutide 24 mg Q2W - Bofanglutide 48 mg Q2W  0.813388 2.926533 334
##  Bofanglutide 24 mg Q2W - Bofanglutide 24 mg QW   1.638852 3.059646 334
##  Bofanglutide 48 mg Q2W - Bofanglutide 24 mg QW   0.825464 2.911261 334
##   lower.CL  upper.CL t.ratio p.value
##  -2.403994 14.258560   2.039  0.3224
##  -1.365335 15.209064   2.394  0.1612
##  -0.690340 15.972214   2.629  0.0932
##   0.572008 16.336643   3.074  0.0275
##   0.992590 17.566989   3.210  0.0181
##  -7.775412  9.764575   0.325  0.9995
##  -7.098002 10.525310   0.557  0.9936
##  -5.861405 10.915489   0.863  0.9548
##  -5.417487 12.122499   1.096  0.8829
##  -8.050921  9.489066   0.235  0.9999
##  -6.812210  9.877132   0.526  0.9951
##  -6.370207 11.086056   0.774  0.9717
##  -7.575059  9.201835   0.278  0.9998
##  -7.131141 10.408845   0.536  0.9947
##  -7.519207  9.170135   0.284  0.9998
## 
## time = W20:
##  contrast                                         estimate       SE  df
##  Placebo - Bofanglutide 12 mg Q2W                 6.978240 2.872918 334
##  Placebo - Bofanglutide 18 mg Q2W                 8.385191 2.857719 334
##  Placebo - Bofanglutide 24 mg Q2W                 9.284990 2.872918 334
##  Placebo - Bofanglutide 48 mg Q2W                10.452807 2.718101 334
##  Placebo - Bofanglutide 24 mg QW                 11.211286 2.857719 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 18 mg Q2W  1.406952 3.024203 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg Q2W  2.306750 3.038570 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 48 mg Q2W  3.474567 2.892632 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg QW   4.233046 3.024203 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg Q2W  0.899798 3.024203 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 48 mg Q2W  2.067616 2.877537 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg QW   2.826094 3.009768 334
##  Bofanglutide 24 mg Q2W - Bofanglutide 48 mg Q2W  1.167817 2.892632 334
##  Bofanglutide 24 mg Q2W - Bofanglutide 24 mg QW   1.926296 3.024203 334
##  Bofanglutide 48 mg Q2W - Bofanglutide 24 mg QW   0.758479 2.877537 334
##   lower.CL  upper.CL t.ratio p.value
##  -1.256529 15.213008   2.429  0.1493
##   0.193990 16.576393   2.934  0.0413
##   1.050221 17.519758   3.232  0.0169
##   2.661797 18.243817   3.846  0.0020
##   3.020084 19.402487   3.923  0.0015
##  -7.261451 10.075354   0.465  0.9973
##  -6.402833 11.016333   0.759  0.9740
##  -4.816708 11.765843   1.201  0.8361
##  -4.435356 12.901449   1.400  0.7272
##  -7.768604  9.568201   0.298  0.9997
##  -6.180392 10.315623   0.719  0.9796
##  -5.800931 11.453120   0.939  0.9361
##  -7.123458  9.459093   0.404  0.9986
##  -6.742106 10.594699   0.637  0.9881
##  -7.489528  9.006486   0.264  0.9998
## 
## time = W24:
##  contrast                                         estimate       SE  df
##  Placebo - Bofanglutide 12 mg Q2W                 8.032343 2.784156 334
##  Placebo - Bofanglutide 18 mg Q2W                 9.525407 2.769427 334
##  Placebo - Bofanglutide 24 mg Q2W                10.547266 2.784156 334
##  Placebo - Bofanglutide 48 mg Q2W                12.651087 2.634123 334
##  Placebo - Bofanglutide 24 mg QW                 13.224444 2.769427 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 18 mg Q2W  1.493064 2.930767 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg Q2W  2.514923 2.944690 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 48 mg Q2W  4.618744 2.803261 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg QW   5.192102 2.930767 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg Q2W  1.021859 2.930767 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 48 mg Q2W  3.125680 2.788632 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg QW   3.699038 2.916778 334
##  Bofanglutide 24 mg Q2W - Bofanglutide 48 mg Q2W  2.103821 2.803261 334
##  Bofanglutide 24 mg Q2W - Bofanglutide 24 mg QW   2.677179 2.930767 334
##  Bofanglutide 48 mg Q2W - Bofanglutide 24 mg QW   0.573358 2.788632 334
##   lower.CL  upper.CL t.ratio p.value
##   0.051996 16.012689   2.885  0.0475
##   1.587281 17.463532   3.439  0.0085
##   2.566919 18.527612   3.788  0.0025
##   5.100788 20.201385   4.803 <0.0001
##   5.286319 21.162570   4.775 <0.0001
##  -6.907519  9.893647   0.509  0.9958
##  -5.925568 10.955414   0.854  0.9569
##  -3.416364 12.653852   1.648  0.5675
##  -3.208481 13.592685   1.772  0.4856
##  -7.378724  9.422442   0.349  0.9993
##  -4.867496 11.118857   1.121  0.8725
##  -4.661447 12.059522   1.268  0.8022
##  -5.931287 10.138929   0.750  0.9753
##  -5.723405 11.077762   0.913  0.9429
##  -7.419819  8.566534   0.206  0.9999
## 
## time = W30:
##  contrast                                         estimate       SE  df
##  Placebo - Bofanglutide 12 mg Q2W                 8.393928 2.856828 334
##  Placebo - Bofanglutide 18 mg Q2W                11.004447 2.841714 334
##  Placebo - Bofanglutide 24 mg Q2W                12.280371 2.856828 334
##  Placebo - Bofanglutide 48 mg Q2W                14.613518 2.702878 334
##  Placebo - Bofanglutide 24 mg QW                 14.886598 2.841714 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 18 mg Q2W  2.610519 3.007265 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg Q2W  3.886442 3.021552 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 48 mg Q2W  6.219590 2.876432 334
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg QW   6.492669 3.007265 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg Q2W  1.275924 3.007265 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 48 mg Q2W  3.609072 2.861421 334
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg QW   3.882151 2.992911 334
##  Bofanglutide 24 mg Q2W - Bofanglutide 48 mg Q2W  2.333148 2.876432 334
##  Bofanglutide 24 mg Q2W - Bofanglutide 24 mg QW   2.606227 3.007265 334
##  Bofanglutide 48 mg Q2W - Bofanglutide 24 mg QW   0.273079 2.861421 334
##   lower.CL  upper.CL t.ratio p.value
##   0.205280 16.582576   2.938  0.0409
##   2.859122 19.149772   3.872  0.0018
##   4.091723 20.469019   4.299  0.0003
##   6.866144 22.360893   5.407 <0.0001
##   6.741273 23.031923   5.239 <0.0001
##  -6.009335 11.230372   0.868  0.9538
##  -4.774361 12.547245   1.286  0.7925
##  -2.025249 14.464429   2.162  0.2584
##  -2.127184 15.112523   2.159  0.2600
##  -7.343930  9.895777   0.424  0.9982
##  -4.592741 11.810884   1.261  0.8058
##  -4.696557 12.460859   1.297  0.7865
##  -5.911691 10.577987   0.811  0.9654
##  -6.013626 11.226080   0.867  0.9541
##  -7.928733  8.474892   0.095  1.0000
## 
## Confidence level used: 0.95 
## Conf-level adjustment: tukey method for comparing a family of 6 estimates 
## P value adjustment: tukey method for comparing a family of 6 estimates 
## treatment = Placebo:
##  time   emmean       SE  df lower.CL  upper.CL
##  W0   95.60005 2.056410 334 91.55490  99.64519
##  W4   95.03509 1.997029 334 91.10675  98.96343
##  W8   94.43308 1.954027 334 90.58933  98.27683
##  W12  94.15509 1.913918 334 90.39024  97.91994
##  W16  94.34009 1.929497 334 90.54459  98.13559
##  W20  94.32868 1.907146 334 90.57715  98.08021
##  W24  94.53573 1.848223 334 90.90010  98.17135
##  W30  94.64341 1.896465 334 90.91289  98.37393
## 
## treatment = Bofanglutide 12 mg Q2W:
##  time   emmean       SE  df lower.CL  upper.CL
##  W0   95.59994 2.316754 334 91.04267 100.15721
##  W4   93.35688 2.249855 334 88.93121  97.78256
##  W8   91.07038 2.201409 334 86.74001  95.40076
##  W12  89.74763 2.156223 334 85.50615  93.98912
##  W16  88.41281 2.173774 334 84.13679  92.68882
##  W20  87.35044 2.148593 334 83.12396  91.57692
##  W24  86.50338 2.082210 334 82.40749  90.59928
##  W30  86.24948 2.136560 334 82.04667  90.45229
## 
## treatment = Bofanglutide 18 mg Q2W:
##  time   emmean       SE  df lower.CL  upper.CL
##  W0   95.60002 2.294794 334 91.08595 100.11409
##  W4   92.94379 2.228529 334 88.56007  97.32751
##  W8   90.54409 2.180543 334 86.25477  94.83342
##  W12  89.16662 2.135784 334 84.96534  93.36791
##  W16  87.41823 2.153169 334 83.18274  91.65371
##  W20  85.94349 2.128227 334 81.75707  90.12991
##  W24  85.01032 2.062473 334 80.95325  89.06740
##  W30  83.63896 2.116308 334 79.47599  87.80193
## 
## treatment = Bofanglutide 24 mg Q2W:
##  time   emmean       SE  df lower.CL  upper.CL
##  W0   95.60002 2.316754 334 91.04275 100.15729
##  W4   92.63042 2.249855 334 88.20475  97.05610
##  W8   90.05506 2.201409 334 85.72468  94.38543
##  W12  88.57467 2.156223 334 84.33319  92.81616
##  W16  86.69915 2.173774 334 82.42314  90.97517
##  W20  85.04369 2.148593 334 80.81721  89.27017
##  W24  83.98846 2.082210 334 79.89256  88.08436
##  W30  82.36304 2.136560 334 78.16023  86.56585
## 
## treatment = Bofanglutide 48 mg Q2W:
##  time   emmean       SE  df lower.CL  upper.CL
##  W0   95.60002 2.088294 334 91.49215  99.70788
##  W4   92.43312 2.027992 334 88.44388  96.42237
##  W8   89.76919 1.984324 334 85.86584  93.67253
##  W12  87.99405 1.943593 334 84.17082  91.81727
##  W16  85.88577 1.959414 334 82.03142  89.74011
##  W20  83.87588 1.936716 334 80.06618  87.68557
##  W24  81.88464 1.876879 334 78.19265  85.57663
##  W30  80.02989 1.925869 334 76.24153  83.81825
## 
## treatment = Bofanglutide 24 mg QW:
##  time   emmean       SE  df lower.CL  upper.CL
##  W0   95.59994 2.294794 334 91.08587 100.11401
##  W4   92.00428 2.228529 334 87.62056  96.38800
##  W8   89.15998 2.180543 334 84.87065  93.44931
##  W12  87.06358 2.135784 334 82.86230  91.26487
##  W16  85.06030 2.153169 334 80.82482  89.29578
##  W20  83.11740 2.128227 334 78.93098  87.30381
##  W24  81.31128 2.062473 334 77.25421  85.36836
##  W30  79.75681 2.116308 334 75.59384  83.91978
## 
## Confidence level used: 0.95 
## treatment = Placebo:
##  contrast   estimate        SE  df  lower.CL  upper.CL t.ratio p.value
##  W0 - W4    0.564955 0.2575529 334 -0.246099  1.376008   2.194  0.6370
##  W0 - W8    1.166970 0.3193125 334  0.161431  2.172508   3.655  0.0084
##  W0 - W12   1.444955 0.4089557 334  0.157123  2.732786   3.533  0.0126
##  W0 - W16   1.259955 0.4507161 334 -0.159384  2.679293   2.795  0.1371
##  W0 - W20   1.271364 0.4862115 334 -0.259752  2.802480   2.615  0.2240
##  W0 - W24   1.064318 0.5645604 334 -0.713524  2.842160   1.885  1.0000
##  W0 - W30   0.956636 0.6292416 334 -1.024892  2.938165   1.520  1.0000
##  W4 - W8    0.602015 0.2306387 334 -0.124283  1.328313   2.610  0.2240
##  W4 - W12   0.880000 0.2911476 334 -0.036845  1.796845   3.023  0.0702
##  W4 - W16   0.695000 0.3210209 334 -0.315918  1.705918   2.165  0.6531
##  W4 - W20   0.706409 0.3567087 334 -0.416893  1.829711   1.980  0.9697
##  W4 - W24   0.499364 0.4237984 334 -0.835209  1.833936   1.178  1.0000
##  W4 - W30   0.391682 0.4832275 334 -1.130037  1.913401   0.811  1.0000
##  W8 - W12   0.277985 0.2733867 334 -0.582930  1.138900   1.017  1.0000
##  W8 - W16   0.092985 0.3146368 334 -0.897830  1.083799   0.296  1.0000
##  W8 - W20   0.104394 0.3338975 334 -0.947074  1.155862   0.313  1.0000
##  W8 - W24  -0.102652 0.3938832 334 -1.343019  1.137716  -0.261  1.0000
##  W8 - W30  -0.210333 0.4540221 334 -1.640082  1.219416  -0.463  1.0000
##  W12 - W16 -0.185000 0.3168306 334 -1.182723  0.812723  -0.584  1.0000
##  W12 - W20 -0.173591 0.3413907 334 -1.248656  0.901474  -0.508  1.0000
##  W12 - W24 -0.380636 0.3829282 334 -1.586506  0.825233  -0.994  1.0000
##  W12 - W30 -0.488318 0.4333379 334 -1.852931  0.876295  -1.127  1.0000
##  W16 - W20  0.011409 0.3386714 334 -1.055092  1.077910   0.034  1.0000
##  W16 - W24 -0.195636 0.3855890 334 -1.409885  1.018612  -0.507  1.0000
##  W16 - W30 -0.303318 0.4365030 334 -1.677898  1.071262  -0.695  1.0000
##  W20 - W24 -0.207045 0.4049823 334 -1.482365  1.068274  -0.511  1.0000
##  W20 - W30 -0.314727 0.4360342 334 -1.687831  1.058377  -0.722  1.0000
##  W24 - W30 -0.107682 0.4335492 334 -1.472960  1.257597  -0.248  1.0000
## 
## treatment = Bofanglutide 12 mg Q2W:
##  contrast   estimate        SE  df  lower.CL  upper.CL t.ratio p.value
##  W0 - W4    2.243058 0.2901595 334  1.329324  3.156791   7.730 <0.0001
##  W0 - W8    4.529558 0.3597379 334  3.396716  5.662399  12.591 <0.0001
##  W0 - W12   5.852308 0.4607301 334  4.401435  7.303181  12.702 <0.0001
##  W0 - W16   7.187135 0.5077774 334  5.588106  8.786163  14.154 <0.0001
##  W0 - W20   8.249500 0.5477666 334  6.524543  9.974457  15.060 <0.0001
##  W0 - W24   9.096558 0.6360345 334  7.093638 11.099477  14.302 <0.0001
##  W0 - W30   9.350462 0.7089045 334  7.118069 11.582854  13.190 <0.0001
##  W4 - W8    2.286500 0.2598379 334  1.468251  3.104749   8.800 <0.0001
##  W4 - W12   3.609250 0.3280073 334  2.576331  4.642169  11.004 <0.0001
##  W4 - W16   4.944077 0.3616626 334  3.805175  6.082979  13.670 <0.0001
##  W4 - W20   6.006442 0.4018685 334  4.740929  7.271956  14.946 <0.0001
##  W4 - W24   6.853500 0.4774519 334  5.349969  8.357031  14.354 <0.0001
##  W4 - W30   7.107404 0.5444048 334  5.393033  8.821775  13.055 <0.0001
##  W8 - W12   1.322750 0.3079978 334  0.352842  2.292658   4.295  0.0001
##  W8 - W16   2.657577 0.3544703 334  1.541324  3.773830   7.497 <0.0001
##  W8 - W20   3.719942 0.3761694 334  2.535357  4.904527   9.889 <0.0001
##  W8 - W24   4.567000 0.4437494 334  3.169600  5.964400  10.292 <0.0001
##  W8 - W30   4.820904 0.5115019 334  3.210147  6.431661   9.425 <0.0001
##  W12 - W16  1.334827 0.3569418 334  0.210791  2.458863   3.740  0.0011
##  W12 - W20  2.397192 0.3846113 334  1.186023  3.608362   6.233 <0.0001
##  W12 - W24  3.244250 0.4314074 334  1.885716  4.602784   7.520 <0.0001
##  W12 - W30  3.498154 0.4881991 334  1.960779  5.035529   7.165 <0.0001
##  W16 - W20  1.062365 0.3815477 334 -0.139156  2.263887   2.784  0.0227
##  W16 - W24  1.909423 0.4344051 334  0.541449  3.277397   4.395  0.0001
##  W16 - W30  2.163327 0.4917649 334  0.614723  3.711931   4.399  0.0001
##  W20 - W24  0.847058 0.4562536 334 -0.589719  2.283834   1.857  0.1285
##  W20 - W30  1.100962 0.4912367 334 -0.445979  2.647902   2.241  0.0770
##  W24 - W30  0.253904 0.4884371 334 -1.284221  1.792028   0.520  0.6035
## 
## treatment = Bofanglutide 18 mg Q2W:
##  contrast   estimate        SE  df  lower.CL  upper.CL t.ratio p.value
##  W0 - W4    2.656226 0.2874091 334  1.751154  3.561299   9.242 <0.0001
##  W0 - W8    5.055925 0.3563280 334  3.933821  6.178028  14.189 <0.0001
##  W0 - W12   6.433396 0.4563629 334  4.996276  7.870517  14.097 <0.0001
##  W0 - W16   8.181792 0.5029643 334  6.597921  9.765664  16.267 <0.0001
##  W0 - W20   9.656528 0.5425744 334  7.947922 11.365135  17.798 <0.0001
##  W0 - W24  10.589698 0.6300056 334  8.605764 12.573632  16.809 <0.0001
##  W0 - W30  11.961057 0.7021849 334  9.749824 14.172289  17.034 <0.0001
##  W4 - W8    2.399698 0.2573749 334  1.589206  3.210191   9.324 <0.0001
##  W4 - W12   3.777170 0.3248982 334  2.754041  4.800298  11.626 <0.0001
##  W4 - W16   5.525566 0.3582344 334  4.397459  6.653673  15.424 <0.0001
##  W4 - W20   7.000302 0.3980593 334  5.746784  8.253820  17.586 <0.0001
##  W4 - W24   7.933472 0.4729262 334  6.444192  9.422751  16.775 <0.0001
##  W4 - W30   9.304830 0.5392444 334  7.606710 11.002951  17.255 <0.0001
##  W8 - W12   1.377472 0.3050784 334  0.416757  2.338186   4.515 <0.0001
##  W8 - W16   3.125868 0.3511103 334  2.020196  4.231540   8.903 <0.0001
##  W8 - W20   4.600604 0.3726037 334  3.427247  5.773960  12.347 <0.0001
##  W8 - W24   5.533774 0.4395431 334  4.149620  6.917927  12.590 <0.0001
##  W8 - W30   6.905132 0.5066534 334  5.309643  8.500621  13.629 <0.0001
##  W12 - W16  1.748396 0.3535584 334  0.635015  2.861778   4.945 <0.0001
##  W12 - W20  3.223132 0.3809656 334  2.023443  4.422821   8.460 <0.0001
##  W12 - W24  4.156302 0.4273182 334  2.810645  5.501958   9.726 <0.0001
##  W12 - W30  5.527660 0.4835715 334  4.004858  7.050463  11.431 <0.0001
##  W16 - W20  1.474736 0.3779310 334  0.284603  2.664869   3.902  0.0003
##  W16 - W24  2.407906 0.4302874 334  1.052899  3.762913   5.596 <0.0001
##  W16 - W30  3.779264 0.4871035 334  2.245339  5.313189   7.759 <0.0001
##  W20 - W24  0.933170 0.4519289 334 -0.489988  2.356327   2.065  0.0397
##  W20 - W30  2.304528 0.4865804 334  0.772251  3.836806   4.736 <0.0001
##  W24 - W30  1.371358 0.4838073 334 -0.152186  2.894903   2.835  0.0097
## 
## treatment = Bofanglutide 24 mg Q2W:
##  contrast   estimate        SE  df  lower.CL  upper.CL t.ratio p.value
##  W0 - W4    2.969596 0.2901595 334  2.055862  3.883330  10.234 <0.0001
##  W0 - W8    5.544962 0.3597379 334  4.412120  6.677803  15.414 <0.0001
##  W0 - W12   7.025346 0.4607301 334  5.574473  8.476219  15.248 <0.0001
##  W0 - W16   8.900865 0.5077774 334  7.301837 10.499894  17.529 <0.0001
##  W0 - W20  10.556327 0.5477666 334  8.831370 12.281284  19.272 <0.0001
##  W0 - W24  11.611558 0.6360345 334  9.608638 13.614477  18.256 <0.0001
##  W0 - W30  13.236981 0.7089045 334 11.004588 15.469373  18.672 <0.0001
##  W4 - W8    2.575365 0.2598379 334  1.757117  3.393614   9.911 <0.0001
##  W4 - W12   4.055750 0.3280073 334  3.022831  5.088669  12.365 <0.0001
##  W4 - W16   5.931269 0.3616626 334  4.792367  7.070171  16.400 <0.0001
##  W4 - W20   7.586731 0.4018685 334  6.321217  8.852245  18.879 <0.0001
##  W4 - W24   8.641962 0.4774519 334  7.138430 10.145493  18.100 <0.0001
##  W4 - W30  10.267385 0.5444048 334  8.553014 11.981755  18.860 <0.0001
##  W8 - W12   1.480385 0.3079978 334  0.510477  2.450293   4.806 <0.0001
##  W8 - W16   3.355904 0.3544703 334  2.239651  4.472157   9.467 <0.0001
##  W8 - W20   5.011365 0.3761694 334  3.826780  6.195950  13.322 <0.0001
##  W8 - W24   6.066596 0.4437494 334  4.669197  7.463996  13.671 <0.0001
##  W8 - W30   7.692019 0.5115019 334  6.081262  9.302777  15.038 <0.0001
##  W12 - W16  1.875519 0.3569418 334  0.751483  2.999555   5.254 <0.0001
##  W12 - W20  3.530981 0.3846113 334  2.319811  4.742150   9.181 <0.0001
##  W12 - W24  4.586212 0.4314074 334  3.227678  5.944745  10.631 <0.0001
##  W12 - W30  6.211635 0.4881991 334  4.674260  7.749010  12.724 <0.0001
##  W16 - W20  1.655462 0.3815477 334  0.453940  2.856983   4.339 <0.0001
##  W16 - W24  2.710692 0.4344051 334  1.342719  4.078666   6.240 <0.0001
##  W16 - W30  4.336115 0.4917649 334  2.787511  5.884719   8.817 <0.0001
##  W20 - W24  1.055231 0.4562536 334 -0.381546  2.492007   2.313  0.0213
##  W20 - W30  2.680654 0.4912367 334  1.133713  4.227595   5.457 <0.0001
##  W24 - W30  1.625423 0.4884371 334  0.087298  3.163548   3.328  0.0019
## 
## treatment = Bofanglutide 48 mg Q2W:
##  contrast   estimate        SE  df  lower.CL  upper.CL t.ratio p.value
##  W0 - W4    3.166891 0.2615462 334  2.343262  3.990519  12.108 <0.0001
##  W0 - W8    5.830828 0.3242634 334  4.809699  6.851957  17.982 <0.0001
##  W0 - W12   7.605969 0.4152965 334  6.298170  8.913768  18.315 <0.0001
##  W0 - W16   9.714250 0.4577044 334  8.272905 11.155595  21.224 <0.0001
##  W0 - W20  11.724141 0.4937501 334 10.169285 13.278996  23.745 <0.0001
##  W0 - W24  13.715375 0.5733138 334 11.909968 15.520782  23.923 <0.0001
##  W0 - W30  15.570125 0.6389979 334 13.557873 17.582377  24.366 <0.0001
##  W4 - W8    2.663938 0.2342147 334  1.926378  3.401497  11.374 <0.0001
##  W4 - W12   4.439078 0.2956618 334  3.508017  5.370139  15.014 <0.0001
##  W4 - W16   6.547359 0.3259982 334  5.520767  7.573952  20.084 <0.0001
##  W4 - W20   8.557250 0.3622394 334  7.416531  9.697969  23.623 <0.0001
##  W4 - W24  10.548484 0.4303693 334  9.193220 11.903749  24.510 <0.0001
##  W4 - W30  12.403234 0.4907198 334 10.857921 13.948547  25.276 <0.0001
##  W8 - W12   1.775141 0.2776255 334  0.900877  2.649404   6.394 <0.0001
##  W8 - W16   3.883422 0.3195152 334  2.877245  4.889599  12.154 <0.0001
##  W8 - W20   5.893312 0.3390745 334  4.825542  6.961083  17.381 <0.0001
##  W8 - W24   7.884547 0.3999903 334  6.624948  9.144146  19.712 <0.0001
##  W8 - W30   9.739297 0.4610616 334  8.287380 11.191214  21.124 <0.0001
##  W12 - W16  2.108281 0.3217430 334  1.095089  3.121474   6.553 <0.0001
##  W12 - W20  4.118172 0.3466839 334  3.026439  5.209905  11.879 <0.0001
##  W12 - W24  6.109406 0.3888654 334  4.884840  7.333972  15.711 <0.0001
##  W12 - W30  7.964156 0.4400567 334  6.578385  9.349927  18.098 <0.0001
##  W16 - W20  2.009891 0.3439224 334  0.926854  3.092928   5.844 <0.0001
##  W16 - W24  4.001125 0.3915675 334  2.768050  5.234200  10.218 <0.0001
##  W16 - W30  5.855875 0.4432709 334  4.459982  7.251768  13.211 <0.0001
##  W20 - W24  1.991234 0.4112615 334  0.696142  3.286327   4.842 <0.0001
##  W20 - W30  3.845984 0.4427948 334  2.451591  5.240378   8.686 <0.0001
##  W24 - W30  1.854750 0.4402713 334  0.468303  3.241197   4.213 <0.0001
## 
## treatment = Bofanglutide 24 mg QW:
##  contrast   estimate        SE  df  lower.CL  upper.CL t.ratio p.value
##  W0 - W4    3.595660 0.2874091 334  2.690588  4.500733  12.511 <0.0001
##  W0 - W8    6.439962 0.3563280 334  5.317859  7.562065  18.073 <0.0001
##  W0 - W12   8.536358 0.4563629 334  7.099238  9.973479  18.705 <0.0001
##  W0 - W16  10.539642 0.5029643 334  8.955770 12.123513  20.955 <0.0001
##  W0 - W20  12.482547 0.5425744 334 10.773940 14.191154  23.006 <0.0001
##  W0 - W24  14.288660 0.6300056 334 12.304726 16.272595  22.680 <0.0001
##  W0 - W30  15.843132 0.7021849 334 13.631900 18.054364  22.563 <0.0001
##  W4 - W8    2.844302 0.2573749 334  2.033809  3.654794  11.051 <0.0001
##  W4 - W12   4.940698 0.3248982 334  3.917570  5.963826  15.207 <0.0001
##  W4 - W16   6.943981 0.3582344 334  5.815875  8.072088  19.384 <0.0001
##  W4 - W20   8.886887 0.3980593 334  7.633369 10.140405  22.326 <0.0001
##  W4 - W24  10.693000 0.4729262 334  9.203721 12.182279  22.610 <0.0001
##  W4 - W30  12.247472 0.5392444 334 10.549351 13.945592  22.712 <0.0001
##  W8 - W12   2.096396 0.3050784 334  1.135682  3.057111   6.872 <0.0001
##  W8 - W16   4.099679 0.3511103 334  2.994007  5.205351  11.676 <0.0001
##  W8 - W20   6.042585 0.3726037 334  4.869228  7.215941  16.217 <0.0001
##  W8 - W24   7.848698 0.4395431 334  6.464544  9.232852  17.856 <0.0001
##  W8 - W30   9.403170 0.5066534 334  7.807681 10.998659  18.559 <0.0001
##  W12 - W16  2.003283 0.3535584 334  0.889901  3.116665   5.666 <0.0001
##  W12 - W20  3.946189 0.3809656 334  2.746500  5.145877  10.358 <0.0001
##  W12 - W24  5.752302 0.4273182 334  4.406645  7.097958  13.461 <0.0001
##  W12 - W30  7.306774 0.4835715 334  5.783971  8.829576  15.110 <0.0001
##  W16 - W20  1.942906 0.3779310 334  0.752773  3.133038   5.141 <0.0001
##  W16 - W24  3.749019 0.4302874 334  2.394012  5.104026   8.713 <0.0001
##  W16 - W30  5.303491 0.4871035 334  3.769566  6.837415  10.888 <0.0001
##  W20 - W24  1.806113 0.4519289 334  0.382956  3.229271   3.996  0.0002
##  W20 - W30  3.360585 0.4865804 334  1.828307  4.892862   6.907 <0.0001
##  W24 - W30  1.554472 0.4838073 334  0.030927  3.078017   3.213  0.0014
## 
## Confidence level used: 0.95 
## Conf-level adjustment: bonferroni method for 28 estimates 
## P value adjustment: holm method for 28 tests 
##  time_custom treatment_pairwise                               estimate
##  W30-W0      Placebo - Bofanglutide 12 mg Q2W                 8.393825
##  W30-W0      Placebo - Bofanglutide 18 mg Q2W                11.004420
##  W30-W0      Placebo - Bofanglutide 24 mg Q2W                12.280344
##  W30-W0      Placebo - Bofanglutide 48 mg Q2W                14.613489
##  W30-W0      Placebo - Bofanglutide 24 mg QW                 14.886496
##  W30-W0      Bofanglutide 12 mg Q2W - Bofanglutide 18 mg Q2W  2.610595
##  W30-W0      Bofanglutide 12 mg Q2W - Bofanglutide 24 mg Q2W  3.886519
##  W30-W0      Bofanglutide 12 mg Q2W - Bofanglutide 48 mg Q2W  6.219663
##  W30-W0      Bofanglutide 12 mg Q2W - Bofanglutide 24 mg QW   6.492671
##  W30-W0      Bofanglutide 18 mg Q2W - Bofanglutide 24 mg Q2W  1.275924
##  W30-W0      Bofanglutide 18 mg Q2W - Bofanglutide 48 mg Q2W  3.609068
##  W30-W0      Bofanglutide 18 mg Q2W - Bofanglutide 24 mg QW   3.882075
##  W30-W0      Bofanglutide 24 mg Q2W - Bofanglutide 48 mg Q2W  2.333144
##  W30-W0      Bofanglutide 24 mg Q2W - Bofanglutide 24 mg QW   2.606151
##  W30-W0      Bofanglutide 48 mg Q2W - Bofanglutide 24 mg QW   0.273007
##         SE  df  lower.CL  upper.CL t.ratio p.value
##  0.9478875 334  5.591428 11.196222   8.855 <0.0001
##  0.9428726 334  8.216850 13.791991  11.671 <0.0001
##  0.9478875 334  9.477948 15.082741  12.955 <0.0001
##  0.8968073 334 11.962109 17.264869  16.295 <0.0001
##  0.9428726 334 12.098925 17.674066  15.788 <0.0001
##  0.9978022 334 -0.339373  5.560563   2.616  0.0465
##  1.0025424 334  0.922537  6.850501   3.877  0.0009
##  0.9543919 334  3.398037  9.041290   6.517 <0.0001
##  0.9978022 334  3.542703  9.442638   6.507 <0.0001
##  0.9978022 334 -1.674044  4.225892   1.279  0.4038
##  0.9494114 334  0.802166  6.415970   3.801  0.0010
##  0.9930394 334  0.946189  6.817962   3.909  0.0009
##  0.9543919 334 -0.488483  5.154771   2.445  0.0465
##  0.9978022 334 -0.343817  5.556119   2.612  0.0465
##  0.9494114 334 -2.533895  3.079909   0.288  0.7739
## 
## Confidence level used: 0.95 
## Conf-level adjustment: bonferroni method for 15 estimates 
## P value adjustment: holm method for 15 tests
## 
## >>> 11. Linear Mixed Model
## 
## Correlation matrix not shown by default, as p = 48 > 12.
## Use print(x, correlation=TRUE)  or
##     vcov(x)        if you need it
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: weight_kg ~ treatment * time + (1 + week | id)
##    Data: dat_obs
## Control: ctrl
## 
## REML criterion at convergence: 13760.8
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -3.4310 -0.5439  0.0134  0.5609  3.0670 
## 
## Random effects:
##  Groups   Name        Variance  Std.Dev. Corr  
##  id       (Intercept) 268.21717 16.3773        
##           week          0.01671  0.1293  -0.47 
##  Residual               3.51648  1.8752        
## Number of obs: 2720, groups:  id, 340
## 
## Fixed effects:
##                    Estimate Std. Error         df t value Pr(>|t|)    
## (Intercept)        89.19017    0.85059  334.00010 104.857  < 2e-16 ***
## treatment1          5.44373    1.78407  334.00007   3.051 0.002461 ** 
## treatment2          0.59620    1.96088  334.00007   0.304 0.761283    
## treatment3         -0.40698    1.94580  334.00007  -0.209 0.834451    
## treatment4         -1.07086    1.96088  334.00007  -0.546 0.585354    
## treatment5         -2.00611    1.80548  334.00007  -1.111 0.267316    
## time1               6.40982    0.13865  705.27407  46.230  < 2e-16 ***
## time2               3.87709    0.11983 1255.57618  32.356  < 2e-16 ***
## time3               1.64846    0.10527 2214.95625  15.659  < 2e-16 ***
## time4               0.26010    0.09692 2133.25922   2.684 0.007341 ** 
## time5              -1.22078    0.09641 2085.37111 -12.662  < 2e-16 ***
## time6              -2.58024    0.10385 2285.31299 -24.846  < 2e-16 ***
## time7              -3.65120    0.11774 1364.61670 -31.011  < 2e-16 ***
## treatment1:time1   -5.44368    0.29081  705.27406 -18.719  < 2e-16 ***
## treatment2:time1   -0.59625    0.31963  705.27406  -1.865 0.062536 .  
## treatment3:time1    0.40700    0.31717  705.27406   1.283 0.199837    
## treatment4:time1    1.07088    0.31963  705.27406   3.350 0.000850 ***
## treatment5:time1    2.00612    0.29430  705.27406   6.817 2.00e-11 ***
## treatment1:time2   -3.47590    0.25133 1255.57617 -13.830  < 2e-16 ***
## treatment2:time2   -0.30658    0.27624 1255.57617  -1.110 0.267290    
## treatment3:time2    0.28351    0.27412 1255.57617   1.034 0.301214    
## treatment4:time2    0.63402    0.27624 1255.57617   2.295 0.021889 *  
## treatment5:time2    1.37196    0.25435 1255.57617   5.394 8.23e-08 ***
## treatment1:time3   -1.84928    0.22080 2214.95625  -8.375  < 2e-16 ***
## treatment2:time3   -0.36444    0.24269 2214.95625  -1.502 0.133317    
## treatment3:time3    0.11245    0.24082 2214.95625   0.467 0.640593    
## treatment4:time3    0.28729    0.24269 2214.95625   1.184 0.236628    
## treatment5:time3    0.93666    0.22345 2214.95625   4.192 2.88e-05 ***
## treatment1:time4   -0.73891    0.20329 2133.25922  -3.635 0.000285 ***
## treatment2:time4   -0.29884    0.22344 2133.25922  -1.337 0.181226    
## treatment3:time4    0.12333    0.22172 2133.25922   0.556 0.578110    
## treatment4:time4    0.19526    0.22344 2133.25922   0.874 0.382294    
## treatment5:time4    0.54988    0.20573 2133.25922   2.673 0.007580 ** 
## treatment1:time5    0.92697    0.20222 2085.37111   4.584 4.83e-06 ***
## treatment2:time5   -0.15278    0.22226 2085.37111  -0.687 0.491909    
## treatment3:time5   -0.14418    0.22055 2085.37111  -0.654 0.513352    
## treatment4:time5   -0.19938    0.22226 2085.37111  -0.897 0.369792    
## treatment5:time5   -0.07752    0.20464 2085.37111  -0.379 0.704873    
## treatment1:time6    2.27502    0.21782 2285.31298  10.445  < 2e-16 ***
## treatment2:time6    0.14432    0.23940 2285.31298   0.603 0.546689    
## treatment3:time6   -0.25946    0.23756 2285.31298  -1.092 0.274877    
## treatment4:time6   -0.49538    0.23940 2285.31298  -2.069 0.038637 *  
## treatment5:time6   -0.72795    0.22043 2285.31299  -3.302 0.000973 ***
## treatment1:time7    3.55303    0.24695 1364.61669  14.388  < 2e-16 ***
## treatment2:time7    0.36822    0.27142 1364.61669   1.357 0.175127    
## treatment3:time7   -0.12167    0.26934 1364.61669  -0.452 0.651540    
## treatment4:time7   -0.47965    0.27142 1364.61669  -1.767 0.077425 .  
## treatment5:time7   -1.64822    0.24991 1364.61669  -6.595 6.06e-11 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
##             efek     Sum Sq     Mean Sq NumDF    DenDF   F value       Pr(>F)
## 1      treatment   38.12036    7.624072     5  334.000   2.16810 5.731099e-02
## 2           time 8723.55097 1246.221567     7 1909.994 354.15954 0.000000e+00
## 3 treatment:time 1886.22787   53.892225    35 2199.682  15.31334 3.142811e-80
##                      model      AIC      BIC singular
## 1         Random intercept 14342.25 14637.67    FALSE
## 2 Random intercept + slope 13864.80 14172.04    FALSE
## time = W0:
##  treatment                emmean       SE     df lower.CL  upper.CL t.ratio
##  Placebo                95.60005 2.029083 339.22 91.60888  99.59121  47.115
##  Bofanglutide 12 mg Q2W 95.59994 2.285967 339.22 91.10349 100.09640  41.820
##  Bofanglutide 18 mg Q2W 95.60002 2.264299 339.22 91.14618 100.05385  42.221
##  Bofanglutide 24 mg Q2W 95.60002 2.285967 339.22 91.10356 100.09648  41.820
##  Bofanglutide 48 mg Q2W 95.60002 2.060543 339.22 91.54696  99.65307  46.396
##  Bofanglutide 24 mg QW  95.59994 2.264299 339.22 91.14611 100.05378  42.221
##  p.value
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
## 
## time = W4:
##  treatment                emmean       SE     df lower.CL  upper.CL t.ratio
##  Placebo                95.03509 1.999931 340.60 91.10132  98.96886  47.519
##  Bofanglutide 12 mg Q2W 93.35688 2.253125 340.60 88.92509  97.78868  41.434
##  Bofanglutide 18 mg Q2W 92.94379 2.231768 340.60 88.55401  97.33357  41.646
##  Bofanglutide 24 mg Q2W 92.63042 2.253125 340.60 88.19863  97.06221  41.112
##  Bofanglutide 48 mg Q2W 92.43313 2.030939 340.60 88.43836  96.42789  45.513
##  Bofanglutide 24 mg QW  92.00428 2.231768 340.60 87.61450  96.39407  41.225
##  p.value
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
## 
## time = W8:
##  treatment                emmean       SE     df lower.CL  upper.CL t.ratio
##  Placebo                94.43308 1.972401 341.64 90.55350  98.31266  47.877
##  Bofanglutide 12 mg Q2W 91.07038 2.222110 341.64 86.69964  95.44112  40.984
##  Bofanglutide 18 mg Q2W 90.54409 2.201047 341.64 86.21478  94.87340  41.137
##  Bofanglutide 24 mg Q2W 90.05506 2.222110 341.64 85.68432  94.42580  40.527
##  Bofanglutide 48 mg Q2W 89.76919 2.002983 341.64 85.82946  93.70892  44.818
##  Bofanglutide 24 mg QW  89.15998 2.201047 341.64 84.83067  93.48929  40.508
##  p.value
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
## 
## time = W12:
##  treatment                emmean       SE     df lower.CL  upper.CL t.ratio
##  Placebo                94.15509 1.946565 342.30 90.32636  97.98383  48.370
##  Bofanglutide 12 mg Q2W 89.74763 2.193002 342.30 85.43418  94.06109  40.925
##  Bofanglutide 18 mg Q2W 89.16662 2.172215 342.30 84.89405  93.43919  41.049
##  Bofanglutide 24 mg Q2W 88.57467 2.193002 342.30 84.26122  92.88813  40.390
##  Bofanglutide 48 mg Q2W 87.99405 1.976746 342.30 84.10595  91.88214  44.515
##  Bofanglutide 24 mg QW  87.06358 2.172215 342.30 82.79101  91.33616  40.081
##  p.value
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
## 
## time = W16:
##  treatment                emmean       SE     df lower.CL  upper.CL t.ratio
##  Placebo                94.34009 1.922488 342.54 90.55872  98.12146  49.072
##  Bofanglutide 12 mg Q2W 88.41281 2.165878 342.54 84.15271  92.67290  40.821
##  Bofanglutide 18 mg Q2W 87.41823 2.145348 342.54 83.19851  91.63794  40.748
##  Bofanglutide 24 mg Q2W 86.69915 2.165878 342.54 82.43906  90.95925  40.030
##  Bofanglutide 48 mg Q2W 85.88577 1.952296 342.54 82.04577  89.72576  43.992
##  Bofanglutide 24 mg QW  85.06030 2.145348 342.54 80.84059  89.28002  39.649
##  p.value
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
## 
## time = W20:
##  treatment                emmean       SE     df lower.CL  upper.CL t.ratio
##  Placebo                94.32868 1.900238 342.32 90.59107  98.06630  49.640
##  Bofanglutide 12 mg Q2W 87.35044 2.140811 342.32 83.13964  91.56124  40.802
##  Bofanglutide 18 mg Q2W 85.94349 2.120519 342.32 81.77260  90.11438  40.529
##  Bofanglutide 24 mg Q2W 85.04369 2.140811 342.32 80.83289  89.25449  39.725
##  Bofanglutide 48 mg Q2W 83.87588 1.929701 342.32 80.08031  87.67144  43.466
##  Bofanglutide 24 mg QW  83.11740 2.120519 342.32 78.94651  87.28828  39.197
##  p.value
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
## 
## time = W24:
##  treatment                emmean       SE     df lower.CL  upper.CL t.ratio
##  Placebo                94.53573 1.879881 341.62 90.83813  98.23333  50.288
##  Bofanglutide 12 mg Q2W 86.50338 2.117877 341.62 82.33766  90.66910  40.844
##  Bofanglutide 18 mg Q2W 85.01032 2.097802 341.62 80.88409  89.13655  40.524
##  Bofanglutide 24 mg Q2W 83.98846 2.117877 341.62 79.82274  88.15418  39.657
##  Bofanglutide 48 mg Q2W 81.88464 1.909028 341.62 78.12971  85.63957  42.893
##  Bofanglutide 24 mg QW  81.31128 2.097802 341.62 77.18505  85.43752  38.760
##  p.value
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
## 
## time = W30:
##  treatment                emmean       SE     df lower.CL  upper.CL t.ratio
##  Placebo                94.64341 1.853027 339.62 90.99855  98.28826  51.075
##  Bofanglutide 12 mg Q2W 86.24948 2.087623 339.62 82.14318  90.35578  41.315
##  Bofanglutide 18 mg Q2W 83.63896 2.067835 339.62 79.57159  87.70634  40.448
##  Bofanglutide 24 mg Q2W 82.36304 2.087623 339.62 78.25674  86.46934  39.453
##  Bofanglutide 48 mg Q2W 80.02989 1.881758 339.62 76.32852  83.73126  42.529
##  Bofanglutide 24 mg QW  79.75681 2.067835 339.62 75.68943  83.82419  38.570
##  p.value
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
##  <0.0001
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## time = W0:
##  contrast                                         estimate       SE     df
##  Placebo - Bofanglutide 12 mg Q2W                 0.000103 3.056603 339.22
##  Placebo - Bofanglutide 18 mg Q2W                 0.000027 3.040432 339.22
##  Placebo - Bofanglutide 24 mg Q2W                 0.000026 3.056603 339.22
##  Placebo - Bofanglutide 48 mg Q2W                 0.000030 2.891888 339.22
##  Placebo - Bofanglutide 24 mg QW                  0.000102 3.040432 339.22
##  Bofanglutide 12 mg Q2W - Bofanglutide 18 mg Q2W -0.000077 3.217561 339.22
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg Q2W -0.000077 3.232846 339.22
##  Bofanglutide 12 mg Q2W - Bofanglutide 48 mg Q2W -0.000073 3.077578 339.22
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg QW  -0.000001 3.217561 339.22
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg Q2W  0.000000 3.217561 339.22
##  Bofanglutide 18 mg Q2W - Bofanglutide 48 mg Q2W  0.000003 3.061517 339.22
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg QW   0.000075 3.202202 339.22
##  Bofanglutide 24 mg Q2W - Bofanglutide 48 mg Q2W  0.000004 3.077578 339.22
##  Bofanglutide 24 mg Q2W - Bofanglutide 24 mg QW   0.000076 3.217561 339.22
##  Bofanglutide 48 mg Q2W - Bofanglutide 24 mg QW   0.000072 3.061517 339.22
##   lower.CL  upper.CL t.ratio p.value
##  -8.760384  8.760590   0.000  1.0000
##  -8.714112  8.714165   0.000  1.0000
##  -8.760461  8.760513   0.000  1.0000
##  -8.288368  8.288428   0.000  1.0000
##  -8.714037  8.714241   0.000  1.0000
##  -9.221882  9.221728   0.000  1.0000
##  -9.265691  9.265537   0.000  1.0000
##  -8.820675  8.820529   0.000  1.0000
##  -9.221806  9.221804   0.000  1.0000
##  -9.221805  9.221805   0.000  1.0000
##  -8.774568  8.774574   0.000  1.0000
##  -9.177711  9.177862   0.000  1.0000
##  -8.820598  8.820606   0.000  1.0000
##  -9.221729  9.221881   0.000  1.0000
##  -8.774499  8.774643   0.000  1.0000
## 
## time = W4:
##  contrast                                         estimate       SE     df
##  Placebo - Bofanglutide 12 mg Q2W                 1.678206 3.012689 340.60
##  Placebo - Bofanglutide 18 mg Q2W                 2.091298 2.996750 340.60
##  Placebo - Bofanglutide 24 mg Q2W                 2.404668 3.012689 340.60
##  Placebo - Bofanglutide 48 mg Q2W                 2.601966 2.850340 340.60
##  Placebo - Bofanglutide 24 mg QW                  3.030808 2.996750 340.60
##  Bofanglutide 12 mg Q2W - Bofanglutide 18 mg Q2W  0.413092 3.171334 340.60
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg Q2W  0.726462 3.186399 340.60
##  Bofanglutide 12 mg Q2W - Bofanglutide 48 mg Q2W  0.923760 3.033362 340.60
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg QW   1.352602 3.171334 340.60
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg Q2W  0.313369 3.171334 340.60
##  Bofanglutide 18 mg Q2W - Bofanglutide 48 mg Q2W  0.510667 3.017532 340.60
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg QW   0.939509 3.156196 340.60
##  Bofanglutide 24 mg Q2W - Bofanglutide 48 mg Q2W  0.197298 3.033362 340.60
##  Bofanglutide 24 mg Q2W - Bofanglutide 24 mg QW   0.626140 3.171334 340.60
##  Bofanglutide 48 mg Q2W - Bofanglutide 24 mg QW   0.428842 3.017532 340.60
##   lower.CL  upper.CL t.ratio p.value
##  -6.956217 10.312630   0.557  0.9936
##  -6.497444 10.680041   0.698  0.9821
##  -6.229756 11.039091   0.798  0.9677
##  -5.567162 10.771094   0.913  0.9431
##  -5.557934 11.619550   1.011  0.9140
##  -8.676011  9.502195   0.130  1.0000
##  -8.405820  9.858743   0.228  0.9999
##  -7.769914  9.617433   0.305  0.9996
##  -7.736501 10.441705   0.427  0.9982
##  -8.775734  9.402472   0.099  1.0000
##  -8.137637  9.158972   0.169  1.0000
##  -8.106209  9.985227   0.298  0.9997
##  -8.496375  8.890971   0.065  1.0000
##  -8.462963  9.715243   0.197  1.0000
##  -8.219463  9.077147   0.142  1.0000
## 
## time = W8:
##  contrast                                         estimate       SE     df
##  Placebo - Bofanglutide 12 mg Q2W                 3.362691 2.971219 341.64
##  Placebo - Bofanglutide 18 mg Q2W                 3.888981 2.955499 341.64
##  Placebo - Bofanglutide 24 mg Q2W                 4.378018 2.971219 341.64
##  Placebo - Bofanglutide 48 mg Q2W                 4.663888 2.811105 341.64
##  Placebo - Bofanglutide 24 mg QW                  5.273095 2.955499 341.64
##  Bofanglutide 12 mg Q2W - Bofanglutide 18 mg Q2W  0.526290 3.127680 341.64
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg Q2W  1.015327 3.142539 341.64
##  Bofanglutide 12 mg Q2W - Bofanglutide 48 mg Q2W  1.301197 2.991608 341.64
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg QW   1.910403 3.127680 341.64
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg Q2W  0.489037 3.127680 341.64
##  Bofanglutide 18 mg Q2W - Bofanglutide 48 mg Q2W  0.774907 2.975996 341.64
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg QW   1.384113 3.112751 341.64
##  Bofanglutide 24 mg Q2W - Bofanglutide 48 mg Q2W  0.285870 2.991608 341.64
##  Bofanglutide 24 mg Q2W - Bofanglutide 24 mg QW   0.895077 3.127680 341.64
##  Bofanglutide 48 mg Q2W - Bofanglutide 24 mg QW   0.609206 2.975996 341.64
##   lower.CL  upper.CL t.ratio p.value
##  -5.152732 11.878114   1.132  0.8679
##  -4.581390 12.359353   1.316  0.7762
##  -4.137405 12.893441   1.473  0.6815
##  -3.392652 12.720429   1.659  0.5599
##  -3.197277 13.743466   1.784  0.4774
##  -8.437546  9.490126   0.168  1.0000
##  -7.991093 10.021747   0.323  0.9995
##  -7.272659  9.875053   0.435  0.9980
##  -7.053433 10.874240   0.611  0.9902
##  -8.474800  9.452873   0.156  1.0000
##  -7.754206  9.304020   0.260  0.9998
##  -7.536936 10.305162   0.445  0.9978
##  -8.287986  8.859727   0.096  1.0000
##  -8.068760  9.858913   0.286  0.9997
##  -7.919907  9.138320   0.205  1.0000
## 
## time = W12:
##  contrast                                         estimate       SE     df
##  Placebo - Bofanglutide 12 mg Q2W                 4.407456 2.932298 342.30
##  Placebo - Bofanglutide 18 mg Q2W                 4.988468 2.916785 342.30
##  Placebo - Bofanglutide 24 mg Q2W                 5.580418 2.932298 342.30
##  Placebo - Bofanglutide 48 mg Q2W                 6.161044 2.774281 342.30
##  Placebo - Bofanglutide 24 mg QW                  7.091506 2.916785 342.30
##  Bofanglutide 12 mg Q2W - Bofanglutide 18 mg Q2W  0.581012 3.086710 342.30
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg Q2W  1.172962 3.101374 342.30
##  Bofanglutide 12 mg Q2W - Bofanglutide 48 mg Q2W  1.753588 2.952420 342.30
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg QW   2.684050 3.086710 342.30
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg Q2W  0.591950 3.086710 342.30
##  Bofanglutide 18 mg Q2W - Bofanglutide 48 mg Q2W  1.172576 2.937013 342.30
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg QW   2.103038 3.071976 342.30
##  Bofanglutide 24 mg Q2W - Bofanglutide 48 mg Q2W  0.580626 2.952420 342.30
##  Bofanglutide 24 mg Q2W - Bofanglutide 24 mg QW   1.511088 3.086710 342.30
##  Bofanglutide 48 mg Q2W - Bofanglutide 24 mg QW   0.930462 2.937013 342.30
##   lower.CL  upper.CL t.ratio p.value
##  -3.996329 12.811242   1.503  0.6625
##  -3.370856 13.347793   1.710  0.5259
##  -2.823368 13.984203   1.903  0.4020
##  -1.789875 14.111963   2.221  0.2309
##  -1.267818 15.450830   2.431  0.1484
##  -8.265308  9.427332   0.188  1.0000
##  -7.715384 10.061307   0.378  0.9990
##  -6.707865 10.215040   0.594  0.9914
##  -6.162270 11.530370   0.870  0.9535
##  -8.254370  9.438269   0.192  1.0000
##  -7.244720  9.589872   0.399  0.9987
##  -6.701056 10.907132   0.685  0.9836
##  -7.880827  9.042079   0.197  1.0000
##  -7.335232 10.357408   0.490  0.9965
##  -7.486834  9.347758   0.317  0.9996
## 
## time = W16:
##  contrast                                         estimate       SE     df
##  Placebo - Bofanglutide 12 mg Q2W                 5.927283 2.896029 342.54
##  Placebo - Bofanglutide 18 mg Q2W                 6.921864 2.880708 342.54
##  Placebo - Bofanglutide 24 mg Q2W                 7.640937 2.896029 342.54
##  Placebo - Bofanglutide 48 mg Q2W                 8.454325 2.739967 342.54
##  Placebo - Bofanglutide 24 mg QW                  9.279789 2.880708 342.54
##  Bofanglutide 12 mg Q2W - Bofanglutide 18 mg Q2W  0.994581 3.048531 342.54
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg Q2W  1.713654 3.063014 342.54
##  Bofanglutide 12 mg Q2W - Bofanglutide 48 mg Q2W  2.527042 2.915902 342.54
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg QW   3.352506 3.048531 342.54
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg Q2W  0.719073 3.048531 342.54
##  Bofanglutide 18 mg Q2W - Bofanglutide 48 mg Q2W  1.532461 2.900685 342.54
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg QW   2.357925 3.033980 342.54
##  Bofanglutide 24 mg Q2W - Bofanglutide 48 mg Q2W  0.813388 2.915902 342.54
##  Bofanglutide 24 mg Q2W - Bofanglutide 24 mg QW   1.638852 3.048531 342.54
##  Bofanglutide 48 mg Q2W - Bofanglutide 24 mg QW   0.825464 2.900685 342.54
##   lower.CL  upper.CL t.ratio p.value
##  -2.372524 14.227091   2.047  0.3182
##  -1.334032 15.177761   2.403  0.1581
##  -0.658871 15.940745   2.638  0.0910
##   0.601781 16.306869   3.086  0.0265
##   1.023892 17.535686   3.221  0.0174
##  -7.742285  9.731448   0.326  0.9995
##  -7.064718 10.492026   0.559  0.9935
##  -5.829719 10.883803   0.867  0.9541
##  -5.384361 12.089372   1.100  0.8813
##  -8.017794  9.455939   0.236  0.9999
##  -6.780690  9.845612   0.528  0.9950
##  -6.337239 11.053088   0.777  0.9712
##  -7.543373  9.170150   0.279  0.9998
##  -7.098015 10.375719   0.538  0.9946
##  -7.487687  9.138615   0.285  0.9997
## 
## time = W20:
##  contrast                                         estimate       SE     df
##  Placebo - Bofanglutide 12 mg Q2W                 6.978240 2.862513 342.32
##  Placebo - Bofanglutide 18 mg Q2W                 8.385191 2.847368 342.32
##  Placebo - Bofanglutide 24 mg Q2W                 9.284990 2.862513 342.32
##  Placebo - Bofanglutide 48 mg Q2W                10.452807 2.708256 342.32
##  Placebo - Bofanglutide 24 mg QW                 11.211286 2.847368 342.32
##  Bofanglutide 12 mg Q2W - Bofanglutide 18 mg Q2W  1.406952 3.013250 342.32
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg Q2W  2.306750 3.027564 342.32
##  Bofanglutide 12 mg Q2W - Bofanglutide 48 mg Q2W  3.474567 2.882155 342.32
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg QW   4.233046 3.013250 342.32
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg Q2W  0.899798 3.013250 342.32
##  Bofanglutide 18 mg Q2W - Bofanglutide 48 mg Q2W  2.067616 2.867115 342.32
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg QW   2.826094 2.998866 342.32
##  Bofanglutide 24 mg Q2W - Bofanglutide 48 mg Q2W  1.167817 2.882155 342.32
##  Bofanglutide 24 mg Q2W - Bofanglutide 24 mg QW   1.926296 3.013250 342.32
##  Bofanglutide 48 mg Q2W - Bofanglutide 24 mg QW   0.758479 2.867115 342.32
##   lower.CL  upper.CL t.ratio p.value
##  -1.225541 15.182020   2.438  0.1463
##   0.224814 16.545569   2.945  0.0401
##   1.081209 17.488770   3.244  0.0162
##   2.691115 18.214499   3.860  0.0019
##   3.050908 19.371663   3.937  0.0014
##  -7.228831 10.042735   0.467  0.9972
##  -6.370058 10.983558   0.762  0.9736
##  -4.785508 11.734643   1.206  0.8340
##  -4.402737 12.868829   1.405  0.7241
##  -7.735985  9.535581   0.299  0.9997
##  -6.149354 10.284585   0.721  0.9793
##  -5.768468 11.420656   0.942  0.9352
##  -7.092258  9.427893   0.405  0.9986
##  -6.709487 10.562079   0.639  0.9880
##  -7.458491  8.975449   0.265  0.9998
## 
## time = W24:
##  contrast                                         estimate       SE     df
##  Placebo - Bofanglutide 12 mg Q2W                 8.032343 2.831846 341.62
##  Placebo - Bofanglutide 18 mg Q2W                 9.525407 2.816864 341.62
##  Placebo - Bofanglutide 24 mg Q2W                10.547266 2.831846 341.62
##  Placebo - Bofanglutide 48 mg Q2W                12.651087 2.679243 341.62
##  Placebo - Bofanglutide 24 mg QW                 13.224444 2.816864 341.62
##  Bofanglutide 12 mg Q2W - Bofanglutide 18 mg Q2W  1.493064 2.980968 341.62
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg Q2W  2.514923 2.995130 341.62
##  Bofanglutide 12 mg Q2W - Bofanglutide 48 mg Q2W  4.618744 2.851279 341.62
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg QW   5.192102 2.980968 341.62
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg Q2W  1.021859 2.980968 341.62
##  Bofanglutide 18 mg Q2W - Bofanglutide 48 mg Q2W  3.125680 2.836399 341.62
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg QW   3.699038 2.966739 341.62
##  Bofanglutide 24 mg Q2W - Bofanglutide 48 mg Q2W  2.103821 2.851279 341.62
##  Bofanglutide 24 mg Q2W - Bofanglutide 24 mg QW   2.677179 2.980968 341.62
##  Bofanglutide 48 mg Q2W - Bofanglutide 24 mg QW   0.573358 2.836399 341.62
##   lower.CL  upper.CL t.ratio p.value
##  -0.083644 16.148330   2.836  0.0542
##   1.452358 17.598455   3.382  0.0104
##   2.431279 18.663253   3.725  0.0031
##   4.972457 20.329716   4.722 <0.0001
##   5.151396 21.297493   4.695 <0.0001
##  -7.050302 10.036430   0.501  0.9961
##  -6.069030 11.098876   0.840  0.9599
##  -3.552935 12.790423   1.620  0.5859
##  -3.351265 13.735468   1.742  0.5051
##  -7.521507  9.565226   0.343  0.9994
##  -5.003355 11.254715   1.102  0.8804
##  -4.803549 12.201624   1.247  0.8133
##  -6.067859 10.275500   0.738  0.9771
##  -5.866188 11.220545   0.898  0.9468
##  -7.555677  8.702393   0.202  1.0000
## 
## time = W30:
##  contrast                                         estimate       SE     df
##  Placebo - Bofanglutide 12 mg Q2W                 8.393928 2.791394 339.62
##  Placebo - Bofanglutide 18 mg Q2W                11.004447 2.776626 339.62
##  Placebo - Bofanglutide 24 mg Q2W                12.280371 2.791394 339.62
##  Placebo - Bofanglutide 48 mg Q2W                14.613518 2.640970 339.62
##  Placebo - Bofanglutide 24 mg QW                 14.886598 2.776626 339.62
##  Bofanglutide 12 mg Q2W - Bofanglutide 18 mg Q2W  2.610519 2.938386 339.62
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg Q2W  3.886442 2.952345 339.62
##  Bofanglutide 12 mg Q2W - Bofanglutide 48 mg Q2W  6.219590 2.810549 339.62
##  Bofanglutide 12 mg Q2W - Bofanglutide 24 mg QW   6.492669 2.938386 339.62
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg Q2W  1.275924 2.938386 339.62
##  Bofanglutide 18 mg Q2W - Bofanglutide 48 mg Q2W  3.609072 2.795882 339.62
##  Bofanglutide 18 mg Q2W - Bofanglutide 24 mg QW   3.882151 2.924360 339.62
##  Bofanglutide 24 mg Q2W - Bofanglutide 48 mg Q2W  2.333148 2.810549 339.62
##  Bofanglutide 24 mg Q2W - Bofanglutide 24 mg QW   2.606227 2.938386 339.62
##  Bofanglutide 48 mg Q2W - Bofanglutide 24 mg QW   0.273079 2.795882 339.62
##   lower.CL  upper.CL t.ratio p.value
##   0.393608 16.394249   3.007  0.0335
##   3.046453 18.962441   3.963  0.0013
##   4.280050 20.280691   4.399  0.0002
##   7.044323 22.182714   5.533 <0.0001
##   6.928604 22.844592   5.361 <0.0001
##  -5.811090 11.032127   0.888  0.9491
##  -4.575174 12.348059   1.316  0.7759
##  -1.835629 14.274809   2.213  0.2345
##  -1.928939 14.914278   2.210  0.2360
##  -7.145685  9.697533   0.434  0.9980
##  -4.404111 11.622254   1.291  0.7900
##  -4.499259 12.263561   1.328  0.7696
##  -5.722071 10.388367   0.830  0.9618
##  -5.815382 11.027836   0.887  0.9495
##  -7.740103  8.286262   0.098  1.0000
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## Conf-level adjustment: tukey method for comparing a family of 6 estimates 
## P value adjustment: tukey method for comparing a family of 6 estimates
## 
## >>> 12. Diagnostik LMM

## Grafik tersimpan: C:/Users/masdi/OneDrive/Desktop/2611018042_DIMAS_Biostat/hasil_analisis_bofanglutide_2026/G7_qq_residual_lmm.png

## Grafik tersimpan: C:/Users/masdi/OneDrive/Desktop/2611018042_DIMAS_Biostat/hasil_analisis_bofanglutide_2026/G8_residual_vs_prediksi.png

## Grafik tersimpan: C:/Users/masdi/OneDrive/Desktop/2611018042_DIMAS_Biostat/hasil_analisis_bofanglutide_2026/G9_residual_kelompok_waktu.png

## Grafik tersimpan: C:/Users/masdi/OneDrive/Desktop/2611018042_DIMAS_Biostat/hasil_analisis_bofanglutide_2026/G10_qq_random_intercept.png

## Grafik tersimpan: C:/Users/masdi/OneDrive/Desktop/2611018042_DIMAS_Biostat/hasil_analisis_bofanglutide_2026/G11_qq_random_slope.png
## [1] FALSE 
## [1] TRUE
## attr(,"gradient")
## [1] 2.655066e-06
## 
## >>> 12b. G12 Diagnostik LMM tiga panel
## 
## >>> 13. G13 Rerata prediksi LMM dan 95% CI

## Grafik tersimpan: C:/Users/masdi/OneDrive/Desktop/2611018042_DIMAS_Biostat/hasil_analisis_bofanglutide_2026/G13_lmm_estimated_marginal_means.png
## 
## >>> 14. Perubahan W30 dari baseline
##                treatment  n mean_change sd_change   mean_pct   sd_pct
## 1                Placebo 66  -0.9566364  4.325979  -1.150060 4.584279
## 2 Bofanglutide 12 mg Q2W 52  -9.3504615  5.055447  -9.750027 5.024006
## 3 Bofanglutide 18 mg Q2W 53 -11.9610566  5.187674 -12.549959 5.019726
## 4 Bofanglutide 24 mg Q2W 52 -13.2369808  5.604434 -13.659977 4.984154
## 5 Bofanglutide 48 mg Q2W 64 -15.5701250  5.822110 -16.090012 4.808863
## 6  Bofanglutide 24 mg QW 53 -15.8431321  4.539435 -16.689924 4.129423
## 
## >>> 15. Kesimpulan analitik otomatis
## 
## SELESAI. Folder hasil: C:/Users/masdi/OneDrive/Desktop/2611018042_DIMAS_Biostat/hasil_analisis_bofanglutide_2026
## Laporan HTML: C:/Users/masdi/OneDrive/Desktop/2611018042_DIMAS_Biostat/hasil_analisis_bofanglutide_2026/Laporan_Analisis_Bofanglutide_2026.html