如何將R中的函數輸出直接保存到工作負載中而不顯示它？

Question

我寫一個函數中的R是這樣的：

MH <- function(Nsim,var,initial) { 
Total <- 0 
X=rep(runif(1),Nsim) 
X[1]=initial               
for (i in 2:Nsim){ 
    Y=X[i-1]+rnorm(1, mean=0, sd=sqrt(var)) 
    alpha= fx(Y)/fx(X[i-1]) 
    u <- runif(1) 
    X[i]=X[i-1]+(Y-X[i-1])*(u<alpha) 
    Total=Total+as.numeric(u<alpha) 
    } 
    accept=Total/Nsim 
    accept 
} 

我想我的功能輸出爲「接受」，同時節省了向量X到工作區的情況下直接顯示它，這樣我可以在以後做的事情與X 。 我知道這是一個基本問題，但有人可以幫助我嗎？

Answer 1

如何像這樣

MH <- function(Nsim,var,initial) { 
Total <- 0 
X=rep(runif(1),Nsim) 
X[1]=initial               
for (i in 2:Nsim){ 
    Y=X[i-1]+rnorm(1, mean=0, sd=sqrt(var)) 
    alpha= fx(Y)/fx(X[i-1]) 
    u <- runif(1) 
    X[i]=X[i-1]+(Y-X[i-1])*(u<alpha) 
    Total=Total+as.numeric(u<alpha) 
    } 
    accept=Total/Nsim 

    # put output into a list and give it class MH 
    output <- list(accept = accept, X = X) 
    class(output) <- "MH" 
    output 
} 

# Write a custom print function for objects of 
# class "MH" 
print.MH <- function(x, ...){ 
    print(x[["accept"]]) 
} 

out <- MH(10, 1, 1) 
out # only prints the acceptance rate 
out$accept # but can access accept 
out$X # and X individually