У меня проблемы с тем, что iTerm2 не вставляет новые строки, когда я копирую определенное количество строк и вставляю его в любой интерпретатор. Например, если я скопирую приведенный ниже код и вставлю его в интерпретатор R, все строки будут свернуты и произойдет ошибка (аналогично, если я вставлю строку в интерпретатор Python).
Это работает, однако, если я вставлю строку в bash; это также работает, если я вставляю его в интерпретатор в Terminal.app (и pbcopy/pbpaste правильно регистрирует переводы строки), поэтому это должно быть проблемой с iTerm2. Любая помощь приветствуется.
gibbs <- function(theta.0, data, burnin=.25*niters, niters=1e3) {
# Sample using Gibbs.
#
# Args:
# theta.0: Starting value.
# data: list with x and y
# burnin: Number of initial samples to throw away.
# niters: Number of desired iterations.
#
# Returns:
# Matrix with (niters-burnin) columns being samples under the target
# distribution.
library(mvtnorm)
# Declare constants.
n <- length(data$y)
X <- cbind(rep(1, length(data$x)), data$x)
# Initialize output.
theta.out <- matrix(rep(NA, (n+2)*(burnin+niters)), nrow=n+2)
# Initialize proposals.
beta.props <- t(rmvnorm(burnin+niters, mean=rep(0,2), sigma=solve(t(X)%*%X)))
z.props <- matrix(abs(rnorm(n*(burnin+niters), mean=0, sd=1)), nrow=n)
z.props[data$y == 0, ] <- -z.props[data$y == 0, ]
# Initialize parameters.
beta.old <- theta.0[1:2]
z.old <- theta.0[3:(n+2)]
for (i in 1:(burnin+niters)) {
beta.new <- solve(t(X)%*%X)%*%t(X)%*%z.old + beta.props[, i]
z.new <- X%*%beta.new + z.props[, i]
# Save the draw.
theta.out[, i] <- c(z.new, beta.new)
beta.old <- beta.new
z.old <- z.new
}
# Remove the burn-in samples.
theta.out[, -(1:burnin)]
}