R version 3.5.1 (2018-07-02) -- "Feather Spray" Copyright (C) 2018 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. [Workspace loaded from ~/.RData] > getwd() [1] "C:/Users/naranjo/Documents" > setwd("stat4620") > getwd() [1] "C:/Users/naranjo/Documents/stat4620" > > ophrs<-c(359,413,25,130,90,50,50,487,102,194,55,74,97) > mean(ophrs) [1] 163.5385 > dexp(ophrs,rate=1/150) [1] 0.0006088317 0.0004247675 0.0056432115 0.0028023359 0.0036587442 [6] 0.0047768754 0.0047768754 0.0002593578 0.0033774466 0.0018290316 [11] 0.0046202708 0.0040705847 0.0034919256 > prod(dexp(ophrs,rate=1/150)) [1] 3.592636e-35 > prod(dexp(ophrs,rate=1/160)) [1] 3.76489e-35 > prod(dexp(ophrs,rate=1/163)) [1] 3.776437e-35 > prod(dexp(ophrs,rate=1/164)) [1] 3.776509e-35 > prod(dexp(ophrs,rate=1/163.5385)) [1] 3.776704e-35 > prod(dexp(ophrs,rate=1/164)) [1] 3.776509e-35 > store<-rep(0,170) > for (i in 150:170){store[i]<-prod(dexp(ophrs,rate=1/i))} > cbind(seq(150,170),store[150:170]) [,1] [,2] [1,] 150 3.592636e-35 [2,] 151 3.619624e-35 [3,] 152 3.644391e-35 [4,] 153 3.666951e-35 [5,] 154 3.687319e-35 [6,] 155 3.705517e-35 [7,] 156 3.721568e-35 [8,] 157 3.735499e-35 [9,] 158 3.747341e-35 [10,] 159 3.757127e-35 [11,] 160 3.764890e-35 [12,] 161 3.770669e-35 [13,] 162 3.774504e-35 [14,] 163 3.776437e-35 [15,] 164 3.776509e-35 [16,] 165 3.774767e-35 [17,] 166 3.771256e-35 [18,] 167 3.766024e-35 [19,] 168 3.759119e-35 [20,] 169 3.750591e-35 [21,] 170 3.740489e-35