Advertisement
Guest User

Untitled

a guest
Jun 15th, 2019
66
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.85 KB | None | 0 0
  1. 1.00
  2. 0.49 1.00
  3. 0.53 0.57 1.00
  4. 0.49 0.46 0.48 1.00
  5. 0.51 0.53 0.57 0.57 1.00
  6. 0.33 0.30 0.31 0.24 0.38 1.00
  7. 0.32 0.21 0.23 0.22 0.32 0.43 1.00
  8. 0.20 0.16 0.14 0.12 0.17 0.27 0.33 1.00
  9. 0.19 0.08 0.07 0.19 0.23 0.24 0.26 0.25 1.00
  10. 0.30 0.27 0.24 0.21 0.32 0.34 0.54 0.46 0.28 1.00
  11. 0.37 0.35 0.37 0.29 0.36 0.37 0.32 0.29 0.30 0.35 1.00
  12. 0.21 0.20 0.18 0.16 0.27 0.40 0.58 0.45 0.27 0.59 0.31 1.00
  13.  
  14. data1 <- na.omit(as.vector(t(read.table('triangle-data.txt', fill = TRUE))))
  15. pt <- 12
  16. R <- matrix(0, nrow = pt , ncol = pt)
  17. for(i in 1:pt){
  18. R[i, 1:i] <- data1[(i*(i-1)/2 + 1): (i*(i+1)/2)]
  19. }
  20. R <- R + t(R) - diag(rep(1, pt))
  21. R
  22.  
  23. > dput(R)
  24. structure(c(1, 0.49, 0.53, 0.49, 0.51, 0.33, 0.32, 0.2, 0.19,
  25. 0.3, 0.37, 0.21, 0.49, 1, 0.57, 0.46, 0.53, 0.3, 0.21, 0.16,
  26. 0.08, 0.27, 0.35, 0.2, 0.53, 0.57, 1, 0.48, 0.57, 0.31, 0.23,
  27. 0.14, 0.07, 0.24, 0.37, 0.18, 0.49, 0.46, 0.48, 1, 0.57, 0.24,
  28. 0.22, 0.12, 0.19, 0.21, 0.29, 0.16, 0.51, 0.53, 0.57, 0.57, 1,
  29. 0.38, 0.32, 0.17, 0.23, 0.32, 0.36, 0.27, 0.33, 0.3, 0.31, 0.24,
  30. 0.38, 1, 0.43, 0.27, 0.24, 0.34, 0.37, 0.4, 0.32, 0.21, 0.23,
  31. 0.22, 0.32, 0.43, 1, 0.33, 0.26, 0.54, 0.32, 0.58, 0.2, 0.16,
  32. 0.14, 0.12, 0.17, 0.27, 0.33, 1, 0.25, 0.46, 0.29, 0.45, 0.19,
  33. 0.08, 0.07, 0.19, 0.23, 0.24, 0.26, 0.25, 1, 0.28, 0.3, 0.27,
  34. 0.3, 0.27, 0.24, 0.21, 0.32, 0.34, 0.54, 0.46, 0.28, 1, 0.35,
  35. 0.59, 0.37, 0.35, 0.37, 0.29, 0.36, 0.37, 0.32, 0.29, 0.3, 0.35,
  36. 1, 0.31, 0.21, 0.2, 0.18, 0.16, 0.27, 0.4, 0.58, 0.45, 0.27,
  37. 0.59, 0.31, 1), .Dim = c(12L, 12L))
  38.  
  39. a <- sapply(sapply(lapply(readLines("triangle.txt"),
  40. function(x) strsplit(x, " ")), "[", 1),
  41. function(x) na.omit(as.numeric(x)))
  42.  
  43. A <- do.call("rbind", a)
  44.  
  45. A[upper.tri(A)] <- 0
  46. A <- A + t(A) - diag(nrow(A))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement