Skip to content
Prev 385392 / 398506 Next

Export R outputs to SAS dataset

Hi Daniel

Thanks,please find the code and output

 #--------R libraries---------
      library(tidyverse)
      library(MF)


MFSubj(lesion ~ group, calflung)
HLBoot(lesion ~ group,  calflung, compare = c("con", "vac"), b = 100,
          B = 100, alpha = 0.05, hpd = TRUE, bca = FALSE,
          return.boot = FALSE, trace.it = FALSE, seed = NULL)


10000 bootstrap samples
95% confidence intervals
Comparing vac to con


Mitigated Fraction

                observed median  lower  upper
Equal Tailed        0.44 0.4464 0.1360 0.7024
Highest Density     0.44 0.4464 0.1456 0.7088


Hodges-Lehmann

                observed   median     lower    upper
Equal Tailed    -0.07335 -0.07125 -0.170425 -0.01480
Highest Density -0.07335 -0.07125 -0.156350 -0.00975


Quartile Differences (quartiles of vac - quartiles of con)

     observed    median    lower     upper
Q25 -0.041500 -0.041300 -0.10340 -0.000905
Q50 -0.112525 -0.111175 -0.28115  0.019350
Q75 -0.168000 -0.170425 -0.38650  0.030000


Quartiles of con
    observed   median   lower   upper
Q25 0.054000 0.054000 0.01525 0.11275
Q50 0.139275 0.139275 0.06140 0.31000
Q75 0.315000 0.315000 0.17300 0.45250


Quartiles of vac
    observed  median   lower    upper
Q25  0.01250 0.01250 0.00125 0.026000
Q50  0.02675 0.02675 0.01665 0.144575
Q75  0.14700 0.14700 0.02810 0.292000


Best
Jose

On Sun, Aug 23, 2020 at 2:44 AM Daniel Nordlund <djnordlund at gmail.com>
wrote: