如何将函数输出转为矩阵？R语言马尔可夫链循环报错问题

Hey there Maria! Let's work through this problem step by step—first, that "unexpected '}'" error almost always comes from mismatched curly braces or messy indentation that makes it hard to track where blocks start/end. Let's fix that and build a solution that fits your Markov chain needs.

First, Let's Break Down Your Requirements #

You need to:

• Run iterations 1 to 45 for a Markov chain
• Use adjacent results (n-th and (n+1)-th iterations) in calculations
• Only execute those calculations when a value N is an integer

Fixing the Brace Error & Building the Solution #

Let's start with a complete, working example tailored to your use case. I'll include comments to explain each part:

# 1. Initialize storage for your Markov chain states (45 iterations total)
markov_states <- numeric(45)
markov_states[1] <- 1  # Set your initial state here (adjust as needed)

# 2. Define your transition matrix (replace with your actual matrix)
transition_matrix <- matrix(
  c(0.6, 0.4, 0.3, 0.7),  # Example transition probabilities
  nrow = 2, byrow = TRUE
)

# 3. Run iterations 1 to 44 (so we can access n and n+1 up to 45)
for(n in 1:44) {
  # Check if N is an integer (use this if N is a separate variable; if N = n, skip to next line)
  # if(floor(N) == N) {
  
  # If you mean the iteration number n should be integer (which it already is, but just in case)
  if(floor(n) == n) {
    # Calculate the (n+1)-th state using the n-th state
    markov_states[n+1] <- sample(
      x = c(1, 2),  # Your possible states (adjust to match your chain)
      size = 1,
      prob = transition_matrix[markov_states[n], ]
    )
    
    # --------------------------
    # Add your custom calculation here using markov_states[n] and markov_states[n+1]
    cat("Iteration", n, "->", n+1, ":", markov_states[n], "→", markov_states[n+1], "\n")
    # --------------------------
  }
}

Key Notes to Avoid Errors & Meet Your Needs #

1. Fixing the Brace Mismatch:
   Always make sure every { has a corresponding }. Use consistent indentation (2-4 spaces per block) to easily spot missing braces—this is the most common cause of your "unexpected '}'" error.

2. Reliable Integer Check:

   • Use floor(N) == N instead of is.integer(N) because R treats numbers like 5.0 as numeric (not integer type), but you still want to count them as integers.
   • If N is a separate parameter (not the iteration number n), just replace the condition with if(floor(N) == N) { ... }.

3. Accessing Adjacent Iterations:
   We loop from 1 to 44 instead of 1 to 45, so we can always access both the n-th and (n+1)-th states without going out of bounds (since the 45th state is the last one).

Example of a Bad Brace Structure (What Caused Your Error) #

This is a common mistake that triggers the "unexpected '}'" message:

for(n in 1:44)
  if(floor(n) == n) {
    markov_states[n+1] <- sample(...)
  }
# Missing the closing } for the for loop!

内容的提问来源于stack exchange，提问作者Maria Denk