Differences

This shows you the differences between two versions of the page.

--- sna_eg_stanford:lab06 [2019/12/03 06:18] – [Lab 06: Step 4 deleted] hkimscil
+++ sna_eg_stanford:lab06 [2019/12/11 01:40] (current) – [no pdf] hkimscil
@@ Line 39: / Line 39: @@
 library(NetCluster)
+# install.packages(c("sna", "NetCluster"))
 ###
 #2. LOADING AND FORMATTING DATA
@@ Line 55: / Line 56: @@
 m182_friend <- delete_edges(m182_full, E(m182_full)[E(m182_full)$friend_tie==0])
 summary(m182_friend)
+m182_friend[]
 m182_social <- delete_edges(m182_full, E(m182_full)[E(m182_full)$social_tie==0])
 summary(m182_social)
+m182_social[]
 m182_task <- delete_edges(m182_full, E(m182_full)[E(m182_full)$task_tie==0])
 summary(m182_task)
+m182_task[]
 # Look at the plots for each sub-graph
@@ Line 992: / Line 996: @@
 # Look at the plots for each sub-graph
-pdf("6.1_m182_studentnet_friend_social_task_plots.pdf", width = 10)
+# pdf("6.1_m182_studentnet_friend_social_task_plots.pdf", width = 10)
 par(mfrow = c(1,3))
@@ Line 1003: / Line 1007: @@
 task_layout <- layout.fruchterman.reingold(m182_task)
 plot(m182_task, layout=task_layout, main = "task", edge.arrow.size=.5)
-dev.off()
+# dev.off()
+# back to normal window pannel (no partition)
+par(mfrow = c(1,1))
 ###
@@ Line 1084: / Line 1091: @@
 # the "method" parameter.
-pdf("6.2_m182_studentnet_social_hclust.pdf")
+# pdf("6.2_m182_studentnet_social_hclust.pdf")
 m182_task_social_hclust <- hclust(m182_task_social_dist)
 plot(m182_task_social_hclust)
-dev.off()
+# dev.off()
 # cutree() allows us to use the output of hclust() to set
@@ Line 1124: / Line 1131: @@
 # with the observed correlation matrix.
-pdf("6.3_m182_studentnet_task_social_clustered_observed_corrs.pdf")
+# pdf("6.3_m182_studentnet_task_social_clustered_observed_corrs.pdf")
 clustered_observed_cors <-clustConfigurations(num_vertices,m182_task_social_hclust,m182_task_social_cors)
 clustered_observed_cors
 plot(clustered_observed_cors$correlations)
-dev.off()
+# dev.off()
 clustered_observed_cors$correlations
@@ Line 1278: / Line 1285: @@
 nS <- nScree(ev$values, ap$eigen$qevpea)
-pdf("6.6_m182_studentnet_task_social_pca_scree.pdf")
+# pdf("6.6_m182_studentnet_task_social_pca_scree.pdf")
 plotnScree(nS)
@@ Line 1285: / Line 1292: @@
 plotnScree(nS)
 abline(h=1)
-dev.off()
+# dev.off()
@@ Line 1322: / Line 1329: @@
 m182_task_factor_hclust <- hclust(m182_task_factor_dist)
-pdf("6.7_m182_studentnet_task_social_pca_hclust.pdf")
+# pdf("6.7_m182_studentnet_task_social_pca_hclust.pdf")
 plot(m182_task_factor_hclust)
-dev.off()
+# dev.off()
 # And compare to NetCluster based on correlations and triads:
-pdf("6.8_m182_task_cluster_by_correlation_PCA_Triads.pdf")
+# pdf("6.8_m182_task_cluster_by_correlation_PCA_Triads.pdf")
-par(mfrow = c(1,3))
+par(mfrow = c(1,2))
 plot(m182_task_social_hclust, main = "Correlation")
 plot(m182_task_factor_hclust, main = "PCA")
-plot(m182_task_triad_hclust, main = "Triads")
+# the below one that didn't work
-dev.off()
+# plot(m182_task_triad_hclust, main = "Triads")
+# dev.off()