Estou com problemas com o iTerm2 não colando novas linhas quando copio um determinado número de linhas e o colo em qualquer intérprete. Por exemplo, se eu copiar o código abaixo e colá-lo no interpretador R, todas as linhas serão recolhidas e serão excluídas (da mesma forma, se eu colar a string no interpretador Python).
No entanto, funciona se eu colar a string no bash; Ele também funciona se eu colá-lo em um intérprete no Terminal.app (e pbcopy/pbpaste
registra as novas linhas corretamente), por isso tem que ser um problema com iTerm2. Qualquer ajuda é apreciada.
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)]
}