title: Getting summary statistics from model using permutation values in R
body: <p>I am trying to obtain the summary statistics (<code>summary()</code>) for the linear model (below), which uses 1000 permutations of the original dataset to create a 1000 random dataset (large matrix).</p>
  <pre><code>random_model &lt;- rep(NA,1000)
  
  for (i in c(1:1000)) {
    random_data &lt;- final_data
    random_data$weighted_degree &lt;- rowSums(node.perm_1000[i,,],na.rm=T)
    random_model[i] &lt;- coef(lm(weighted_degree ~ age + sex + age*sex, data=random_data))
  }
  </code></pre>
  <p>I am not simply trying to compare the models to get an overall p-value but I want to get a t-value for each of the variables in the model that uses the random permutations as well.</p>