User-defined function isn't running

General
,
1

I'm trying to run a script obtained from an online supplement to a methods paper. The only bits I have to edit are a few seed specifications at the top. Then I try to run the entire script, and nothing else--and nothing happens. It's like submitting inert text. I just get a cursor at the end. As far as I can tell, this is using basic R functions and doesn't rely on any particular packages (and none were recommended with the script). I'm VERY new to R, and really can't see what I am doing wrong. Will paste the script below--any suggestions? I am keeping my specification comments intact in case something about them is creating a problem, but you can otherwise ignore those comments.

I wanted to try to use debug() but that requires a function name and this is a UDF. Suggestions on that are also most welcome.

function () {
      #USER SPECIFICATION PORTION
      alpha<- 0.05 #DESIGNATED ALPHA
      power<- 0.9 #NOMINAL POWER
      beta1<- 1.26 #TREATMENT MEANS - from Luke paper Table 1 Pos affect act-ext & sham
      beta2<- 1.13
      sigsq<- 1.4377 #ERROR VARIANCE - hand-calculated pooled error term from Pos affect BP
      rn21<- 1 #GROUP SIZE RATIO
      tausq1<- 0.46 #COVARIATE VARIANCES - from Table 1 trait extraversion
      tausq2<- 0.46
      #END OF SPECIFICATION

      betad<-beta1-beta2
      sigma<-sqrt(sigsq)
      del<-betad/sigma
      numint<-50
      l<-numint+1
      dd<-1e-6
      coevec<-c(1,rep(c(4,2),numint/2-1),4,1)
      bl<-dd
      bu<-1-dd
      intb<-(bu-bl)/numint
      bvec<-bl+intb*(0:numint)
      kbpowerf<-function(){
        df<-n1+n2-4
        tcrit<-qt(1-alpha/2,df)
        dfk1<-n1-1
        dfk2<-n2-1
        dfk<-n1+n2-2
        dfb1<-dfk1/2
        dfb2<-dfk2/2
        wbpdf<-(intb/3)*coevec*dbeta(bvec,dfb1,dfb2)
        cl<-dd
        cu<-qchisq(1-dd,dfk)
        intc<-(cu-cl)/numint
        cvec<-cl+intc*(0:numint)
        wcpdf<-(intc/3)*coevec*dchisq(cvec,dfk)
        quan<-rep(0,l)
        for (i in seq(l)) {
          b1<-bvec[i]
          b2<-1-b1
          deltakbvec<-del/sqrt((1/(b1*tausq1)+1/(b2*tausq2))/cvec)
          quan[i]<-sum(wcpdf*(pt(-tcrit,df,deltakbvec)+
                                pt(tcrit,df,deltakbvec,lower.tail=FALSE)))
        }
        kbpower<-sum(wbpdf*quan)
      }
      n1<-9
      loop<-0
      kbpower<-0
      while(kbpower<power & loop<1000){
        n1<-n1+1
        n2<-n1*rn21
        loop<-loop+1
        kbpower<-kbpowerf()
      }
      kbn1<-n1
      kbn2<-n2
      print("kbn1,kbn2,kbpower")
      print(c(kbn1,kbn2,kbpower),digits=4)
}

System Information:

  • RStudio Edition: desktop
  • RStudio Version: 1.2.5019
  • OS Version: Mojave
  • R Version: 3.6.1 GUI 1.70 El Capitan build (7684)
  • sessionInfo(): R version 3.6.1 (2019-07-05)
    Platform: x86_64-apple-darwin15.6.0 (64-bit)
    Running under: macOS Mojave 10.14.6

Matrix products: default
BLAS: /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/3.6/Resources/lib/libRlapack.dylib

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

attached base packages:
[1] stats graphics grDevices utils datasets methods base

loaded via a namespace (and not attached):
[1] compiler_3.6.1 tools_3.6.1 packrat_0.5.0

Referred here from support.rstudio.com

2

You appear just to have defined a function but not have assigned it to anything. You also need to call the function. Try this:

my_func <- function() {
  # your code
}

my_func()
3

Eureka. Thank you for not only helping, but being gentle in pointing out what is clearly an obvious oversight. Resolved!

4

No problem. It would be worth your while to become familiar how to adapt your bare function into using parameters (in particular section 19.5):

5

This topic was automatically closed 7 days after the last reply. New replies are no longer allowed.