R version 2.11.0 (2010-04-22) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(56 + ,34 + ,5 + ,155 + ,124.9 + ,65.55 + ,44.07 + ,4 + ,-2 + ,-2 + ,-2 + ,5 + ,8 + ,11 + ,3 + ,70 + ,40 + ,16 + ,137 + ,550.95 + ,191.95 + ,113.6 + ,4 + ,4 + ,4 + ,-2 + ,8 + ,2 + ,11 + ,3 + ,47 + ,28 + ,15 + ,72 + ,431.64 + ,165.86 + ,112.46 + ,4 + ,4 + ,-2 + ,4 + ,4 + ,4 + ,7 + ,6 + ,12 + ,10 + ,0 + ,0 + ,0 + ,0 + ,0 + ,-2 + ,-2 + ,-2 + ,-2 + ,8 + ,0 + ,12 + ,8 + ,31 + ,20 + ,3 + ,75 + ,111.97 + ,62.78 + ,72.72 + ,0 + ,-2 + ,4 + ,4 + ,6 + ,10 + ,7 + ,6 + ,46 + ,41 + ,2 + ,160 + ,222.41 + ,94.36 + ,70.55 + ,4 + ,0 + ,-2 + ,4 + ,10 + ,8 + ,12 + ,12 + ,38 + ,37 + ,15 + ,165 + ,192.3 + ,94.45 + ,68.92 + ,4 + ,4 + ,4 + ,4 + ,6 + ,7 + ,4 + ,7 + ,33 + ,30 + ,9 + ,211 + ,238.98 + ,106.73 + ,71.66 + ,4 + ,4 + ,4 + ,4 + ,9 + ,7 + ,13 + ,13 + ,63 + ,36 + ,9 + ,83 + ,306.95 + ,95.4 + ,72.51 + ,4 + ,4 + ,-2 + ,-2 + ,10 + ,6 + ,10 + ,5 + ,28 + ,15 + ,17 + ,105 + ,377.49 + ,133.03 + ,93.05 + ,4 + ,4 + ,4 + ,4 + ,3 + ,0 + ,8 + ,8 + ,45 + ,31 + ,13 + ,138 + ,433.17 + ,139.95 + ,145.32 + ,4 + ,-2 + ,-2 + ,4 + ,10 + ,8 + ,11 + ,9 + ,24 + ,17 + ,8 + ,63 + ,142.63 + ,75.71 + ,74.55 + ,4 + ,4 + ,4 + ,4 + ,2 + ,6 + ,6 + ,5 + ,36 + ,25 + ,6 + ,139 + ,158.74 + ,73.83 + ,53.7 + ,4 + ,0 + ,4 + ,4 + ,6 + ,8 + ,8 + ,10 + ,30 + ,25 + ,38 + ,98 + ,199.55 + ,79.16 + ,65.84 + ,4 + ,0 + ,-2 + ,4 + ,5 + ,1 + ,6 + ,2 + ,58 + ,52 + ,8 + ,142 + ,224.65 + ,106.73 + ,80.73 + ,4 + ,0 + ,4 + ,4 + ,2 + ,6 + ,13 + ,11 + ,51 + ,34 + ,12 + ,143 + ,270.93 + ,120.9 + ,87.7 + ,4 + ,4 + ,4 + ,4 + ,7 + ,8 + ,14 + ,11 + ,45 + ,42 + ,21 + ,167 + ,241.92 + ,92.99 + ,61.24 + ,0 + ,4 + ,4 + ,4 + ,6 + ,6 + ,5 + ,10 + ,65 + ,40 + ,7 + ,134 + ,224.99 + ,93.1 + ,94.19 + ,4 + ,-2 + ,4 + ,4 + ,4 + ,4 + ,10 + ,10 + ,64 + ,49 + ,8 + ,135 + ,172.84 + ,76.53 + ,84.26 + ,4 + ,4 + ,4 + ,0 + ,4 + ,7 + ,6 + ,7 + ,24 + ,24 + ,3 + ,83 + ,99.61 + ,68.4 + ,49.84 + ,4 + ,0 + ,-2 + ,4 + ,4 + ,10 + ,7 + ,4 + ,51 + ,33 + ,32 + ,180 + ,395.78 + ,167.04 + ,124.37 + ,4 + ,4 + ,4 + ,4 + ,2 + ,2 + ,10 + ,9 + ,32 + ,26 + ,18 + ,122 + ,443.16 + ,133.9 + ,77.13 + ,4 + ,0 + ,4 + ,-2 + ,8 + ,6 + ,10 + ,6 + ,16 + ,10 + ,2 + ,79 + ,98 + ,61.52 + ,66.85 + ,4 + ,4 + ,4 + ,4 + ,5 + ,8 + ,7 + ,3 + ,38 + ,31 + ,12 + ,147 + ,196.75 + ,86.12 + ,57.65 + ,4 + ,-2 + ,-2 + ,4 + ,8 + ,4 + ,7 + ,8 + ,75 + ,52 + ,10 + ,178 + ,301.11 + ,125.24 + ,76.37 + ,4 + ,4 + ,4 + ,4 + ,7 + ,6 + ,11 + ,14 + ,56 + ,31 + ,15 + ,190 + ,281.28 + ,113.75 + ,85.74 + ,4 + ,4 + ,4 + ,4 + ,4 + ,4 + ,8 + ,0 + ,38 + ,34 + ,17 + ,152 + ,310.81 + ,114.51 + ,74.23 + ,4 + ,0 + ,4 + ,4 + ,2 + ,1 + ,6 + ,0 + ,46 + ,43 + ,18 + ,177 + ,370.86 + ,143.62 + ,74.06 + ,4 + ,4 + ,-2 + ,4 + ,8 + ,7 + ,12 + ,6 + ,32 + ,31 + ,14 + ,142 + ,446.36 + ,165.41 + ,90.25 + ,4 + ,-2 + ,-2 + ,4 + ,10 + ,4 + ,14 + ,8 + ,62 + ,60 + ,18 + ,127 + ,336.62 + ,115.19 + ,75.11 + ,4 + ,-2 + ,4 + ,4 + ,6 + ,7 + ,12 + ,7 + ,46 + ,36 + ,6 + ,125 + ,134.66 + ,75.43 + ,69.89 + ,4 + ,-2 + ,-2 + ,-2 + ,8 + ,2 + ,8 + ,5 + ,54 + ,53 + ,17 + ,180 + ,333.21 + ,126.88 + ,71.85 + ,4 + ,0 + ,-2 + ,4 + ,4 + ,10 + ,12 + ,9 + ,43 + ,33 + ,17 + ,108 + ,326.98 + ,113.95 + ,70.19 + ,4 + ,0 + ,4 + ,4 + ,10 + ,9 + ,12 + ,8 + ,74 + ,51 + ,26 + ,175 + ,520.95 + ,215.17 + ,202.68 + ,4 + ,4 + ,4 + ,4 + ,10 + ,9 + ,14 + ,12 + ,24 + ,24 + ,3 + ,137 + ,136.04 + ,72.16 + ,52.71 + ,4 + ,4 + ,4 + ,4 + ,4 + ,10 + ,4 + ,12 + ,42 + ,29 + ,15 + ,168 + ,176.99 + ,87.94 + ,53.83 + ,4 + ,0 + ,4 + ,4 + ,4 + ,2 + ,7 + ,7 + ,33 + ,31 + ,13 + ,178 + ,448.54 + ,173.34 + ,97.73 + ,4 + ,-2 + ,-2 + ,4 + ,6 + ,8 + ,12 + ,10 + ,44 + ,41 + ,13 + ,111 + ,171.14 + ,81.27 + ,62.99 + ,4 + ,0 + ,-2 + ,4 + ,5 + ,4 + ,3 + ,7 + ,53 + ,37 + ,8 + ,163 + ,459.07 + ,170.49 + ,132.84 + ,4 + ,4 + ,4 + ,4 + ,6 + ,8 + ,1 + ,4 + ,31 + ,31 + ,9 + ,159 + ,158.74 + ,80.87 + ,68.29 + ,4 + ,-2 + ,4 + ,4 + ,10 + ,2 + ,8 + ,0 + ,41 + ,17 + ,17 + ,160 + ,447.5 + ,168.71 + ,104.5 + ,4 + ,-2 + ,4 + ,4 + ,9 + ,7 + ,7 + ,5 + ,51 + ,37 + ,28 + ,140 + ,581.06 + ,200.45 + ,110.77 + ,4 + ,0 + ,4 + ,4 + ,6 + ,4 + ,8 + ,6 + ,35 + ,31 + ,17 + ,111 + ,386.25 + ,133.57 + ,87.19 + ,4 + ,0 + ,4 + ,4 + ,3 + ,2 + ,8 + ,10 + ,28 + ,25 + ,17 + ,140 + ,359.99 + ,126.07 + ,78 + ,4 + ,0 + ,0 + ,4 + ,6 + ,0 + ,12 + ,12 + ,41 + ,33 + ,8 + ,133 + ,280.26 + ,112.01 + ,110.59 + ,4 + ,4 + ,4 + ,4 + ,8 + ,1 + ,8 + ,4 + ,45 + ,39 + ,10 + ,117 + ,299.28 + ,105.04 + ,56.78 + ,4 + ,4 + ,0 + ,4 + ,9 + ,4 + ,12 + ,8 + ,26 + ,14 + ,12 + ,50 + ,170.84 + ,94.55 + ,86.38 + ,4 + ,-2 + ,-2 + ,4 + ,10 + ,8 + ,8 + ,0 + ,43 + ,29 + ,10 + ,175 + ,327.96 + ,116.07 + ,89.85 + ,4 + ,-2 + ,4 + ,4 + ,6 + ,4 + ,9 + ,7 + ,1 + ,0 + ,6 + ,148 + ,144.41 + ,73.88 + ,51.52 + ,4 + ,4 + ,4 + ,4 + ,7 + ,10 + ,10 + ,11 + ,58 + ,46 + ,28 + ,181 + ,392.62 + ,154.52 + ,80.91 + ,4 + ,4 + ,-2 + ,4 + ,10 + ,8 + ,14 + ,12 + ,21 + ,15 + ,8 + ,153 + ,154.13 + ,73.62 + ,55.9 + ,4 + ,-2 + ,4 + ,4 + ,10 + ,0 + ,8 + ,10 + ,26 + ,20 + ,28 + ,76 + ,261.71 + ,84.64 + ,81.66 + ,4 + ,0 + ,4 + ,0 + ,10 + ,6 + ,0 + ,4 + ,35 + ,34 + ,9 + ,75 + ,187.08 + ,86.4 + ,43.71 + ,4 + ,4 + ,-2 + ,4 + ,9 + ,5 + ,7 + ,6 + ,50 + ,40 + ,12 + ,125 + ,259.49 + ,115.13 + ,57.96 + ,4 + ,0 + ,-2 + ,-2 + ,10 + ,4 + ,11 + ,12 + ,58 + ,57 + ,26 + ,166 + ,296.36 + ,112.81 + ,83.42 + ,4 + ,-2 + ,4 + ,4 + ,4 + ,6 + ,11 + ,8 + ,25 + ,24 + ,18 + ,152 + ,506.74 + ,195.97 + ,111.89 + ,-2 + ,0 + ,4 + ,4 + ,6 + ,2 + ,12 + ,6 + ,76 + ,54 + ,26 + ,178 + ,553.98 + ,174.91 + ,107.6 + ,0 + ,0 + ,4 + ,-2 + ,10 + ,10 + ,14 + ,11 + ,47 + ,42 + ,21 + ,126 + ,257.89 + ,107.21 + ,72.74 + ,4 + ,-2 + ,4 + ,4 + ,0 + ,0 + ,0 + ,0 + ,42 + ,29 + ,15 + ,190 + ,474.79 + ,171.7 + ,122.74 + ,4 + ,4 + ,-2 + ,4 + ,6 + ,8 + ,8 + ,8 + ,3 + ,0 + ,17 + ,129 + ,309.51 + ,128.87 + ,76.59 + ,4 + ,0 + ,4 + ,4 + ,2 + ,2 + ,5 + ,2 + ,80 + ,73 + ,3 + ,107 + ,113.72 + ,66.61 + ,77.73 + ,4 + ,4 + ,-2 + ,4 + ,6 + ,6 + ,12 + ,5 + ,34 + ,33 + ,11 + ,97 + ,257.55 + ,88.36 + ,62.56 + ,4 + ,0 + ,4 + ,4 + ,4 + ,3 + ,9 + ,2 + ,60 + ,43 + ,6 + ,71 + ,303.42 + ,120.15 + ,79.75 + ,-2 + ,0 + ,0 + ,4 + ,5 + ,6 + ,11 + ,10 + ,43 + ,29 + ,2 + ,92 + ,88.7 + ,62.14 + ,69.51 + ,0 + ,0 + ,0 + ,4 + ,4 + ,3 + ,4 + ,4 + ,23 + ,22 + ,7 + ,145 + ,70.08 + ,54.52 + ,59.7 + ,0 + ,4 + ,0 + ,4 + ,5 + ,2 + ,4 + ,0 + ,52 + ,48 + ,5 + ,85 + ,148.39 + ,83.7 + ,45.19 + ,-2 + ,-2 + ,-2 + ,-2 + ,5 + ,2 + ,7 + ,0 + ,14 + ,10 + ,16 + ,159 + ,246.87 + ,114.4 + ,69.51 + ,4 + ,-2 + ,-2 + ,-2 + ,4 + ,0 + ,5 + ,3 + ,49 + ,36 + ,15 + ,197 + ,316.33 + ,114.45 + ,59.29 + ,4 + ,0 + ,4 + ,4 + ,10 + ,3 + ,8 + ,5 + ,53 + ,38 + ,7 + ,171 + ,169.74 + ,97.16 + ,76.27 + ,4 + ,-2 + ,4 + ,4 + ,6 + ,6 + ,9 + ,8 + ,52 + ,49 + ,10 + ,128 + ,283.99 + ,98.38 + ,67.38 + ,4 + ,4 + ,0 + ,4 + ,10 + ,10 + ,12 + ,5 + ,37 + ,28 + ,10 + ,140 + ,150.44 + ,73.91 + ,49.34 + ,4 + ,4 + ,-2 + ,4 + ,4 + ,5 + ,9 + ,5 + ,42 + ,35 + ,13 + ,159 + ,176.16 + ,75.88 + ,67.48 + ,4 + ,4 + ,-2 + ,4 + ,6 + ,2 + ,10 + ,0 + ,49 + ,43 + ,22 + ,180 + ,392.61 + ,158.84 + ,119.72 + ,4 + ,4 + ,4 + ,4 + ,10 + ,2 + ,10 + ,5 + ,40 + ,39 + ,9 + ,145 + ,235.59 + ,105.18 + ,95.47 + ,4 + ,4 + ,4 + ,4 + ,4 + ,2 + ,10 + ,11 + ,50 + ,33 + ,22 + ,162 + ,1632.54 + ,472.94 + ,291.98 + ,4 + ,4 + ,4 + ,4 + ,10 + ,10 + ,15 + ,14 + ,39 + ,37 + ,9 + ,144 + ,226.23 + ,96.33 + ,56.36 + ,4 + ,-2 + ,-2 + ,4 + ,10 + ,8 + ,8 + ,4 + ,39 + ,17 + ,8 + ,99 + ,245.68 + ,95.95 + ,61.84 + ,-2 + ,-2 + ,0 + ,4 + ,6 + ,2 + ,8 + ,0 + ,0 + ,0 + ,14 + ,85 + ,405.24 + ,136.41 + ,119.15 + ,-2 + ,-2 + ,-2 + ,4 + ,2 + ,0 + ,3 + ,5 + ,10 + ,10 + ,6 + ,100 + ,167.91 + ,74.16 + ,56.85 + ,4 + ,4 + ,4 + ,4 + ,6 + ,4 + ,13 + ,0 + ,41 + ,39 + ,16 + ,167 + ,361.29 + ,165.1 + ,120.98 + ,4 + ,-2 + ,-2 + ,4 + ,9 + ,2 + ,6 + ,6 + ,28 + ,22 + ,13 + ,160 + ,275.59 + ,98.39 + ,70.78 + ,4 + ,0 + ,0 + ,4 + ,6 + ,9 + ,6 + ,11 + ,43 + ,22 + ,10 + ,143 + ,287.48 + ,104.35 + ,58.13 + ,4 + ,4 + ,4 + ,4 + ,7 + ,5 + ,11 + ,8 + ,37 + ,34 + ,20 + ,169 + ,216.41 + ,96.09 + ,63.56 + ,4 + ,-2 + ,4 + ,-2 + ,10 + ,4 + ,8 + ,3 + ,15 + ,12 + ,23 + ,97 + ,376.42 + ,133.68 + ,51.49 + ,4 + ,4 + ,-2 + ,-2 + ,2 + ,2 + ,10 + ,3 + ,38 + ,29 + ,29 + ,152 + ,509.97 + ,160.12 + ,111.18 + ,4 + ,4 + ,4 + ,4 + ,9 + ,8 + ,14 + ,12 + ,22 + ,18 + ,13 + ,94 + ,289.48 + ,87.13 + ,75.79 + ,0 + ,0 + ,0 + ,0 + ,10 + ,4 + ,10 + ,7 + ,35 + ,33 + ,15 + ,137 + ,520.97 + ,163.58 + ,118.4 + ,-2 + ,-2 + ,4 + ,4 + ,6 + ,2 + ,9 + ,0 + ,7 + ,7 + ,0 + ,0 + ,0 + ,0 + ,0 + ,4 + ,-2 + ,-2 + ,4 + ,2 + ,5 + ,14 + ,0 + ,43 + ,31 + ,8 + ,183 + ,209.14 + ,91.17 + ,79.62 + ,4 + ,4 + ,4 + ,4 + ,7 + ,0 + ,10 + ,9 + ,29 + ,39 + ,9 + ,106 + ,194.69 + ,89.64 + ,64.99 + ,4 + ,0 + ,-2 + ,4 + ,8 + ,5 + ,12 + ,2 + ,50 + ,44 + ,19 + ,166 + ,323.56 + ,128.29 + ,84.88 + ,4 + ,4 + ,-2 + ,4 + ,10 + ,2 + ,14 + ,2 + ,44 + ,27 + ,28 + ,139 + ,507.98 + ,166.57 + ,70.66 + ,4 + ,0 + ,4 + ,4 + ,7 + ,8 + ,6 + ,9 + ,45 + ,26 + ,12 + ,155 + ,265.14 + ,109.83 + ,35.8 + ,-2 + ,4 + ,4 + ,4 + ,8 + ,0 + ,7 + ,6 + ,49 + ,35 + ,10 + ,123 + ,119.11 + ,72.36 + ,46.4 + ,4 + ,0 + ,4 + ,-2 + ,2 + ,4 + ,9 + ,3 + ,32 + ,13 + ,7 + ,28 + ,94.96 + ,67.4 + ,61.17 + ,4 + ,-2 + ,4 + ,4 + ,8 + ,4 + ,6 + ,3 + ,82 + ,42 + ,23 + ,141 + ,335.78 + ,126.3 + ,118.59 + ,4 + ,4 + ,4 + ,4 + ,7 + ,7 + ,9 + ,8 + ,59 + ,45 + ,19 + ,162 + ,430.35 + ,170.66 + ,113.88 + ,4 + ,0 + ,4 + ,4 + ,10 + ,0 + ,15 + ,6 + ,0 + ,0 + ,5 + ,66 + ,90.55 + ,58.46 + ,61.41 + ,4 + ,-2 + ,4 + ,4 + ,6 + ,2 + ,8 + ,4 + ,22 + ,14 + ,25 + ,95 + ,359.52 + ,119.05 + ,88.14 + ,4 + ,4 + ,4 + ,4 + ,7 + ,0 + ,2 + ,0 + ,6 + ,0 + ,8 + ,109 + ,137.46 + ,81.56 + ,67.17 + ,4 + ,4 + ,4 + ,4 + ,6 + ,7 + ,10 + ,8 + ,16 + ,16 + ,3 + ,26 + ,164.69 + ,85.39 + ,89.45 + ,4 + ,0 + ,-2 + ,4 + ,3 + ,6 + ,5 + ,0 + ,35 + ,23 + ,9 + ,141 + ,263.75 + ,101.27 + ,83.41 + ,4 + ,4 + ,4 + ,4 + ,6 + ,5 + ,13 + ,11 + ,41 + ,39 + ,16 + ,148 + ,353.12 + ,142.92 + ,70.79 + ,4 + ,4 + ,-2 + ,4 + ,7 + ,0 + ,12 + ,13 + ,44 + ,29 + ,8 + ,74 + ,332.74 + ,103.99 + ,70.21 + ,4 + ,4 + ,4 + ,4 + ,5 + ,7 + ,8 + ,4 + ,44 + ,36 + ,10 + ,167 + ,224.63 + ,97.87 + ,60.45 + ,4 + ,4 + ,-2 + ,4 + ,6 + ,5 + ,6 + ,2 + ,19 + ,19 + ,6 + ,68 + ,124.4 + ,71.46 + ,51.17 + ,4 + ,-2 + ,4 + ,4 + ,4 + ,0 + ,0 + ,0 + ,27 + ,18 + ,15 + ,177 + ,195.29 + ,85.98 + ,56.93 + ,4 + ,4 + ,4 + ,4 + ,6 + ,0 + ,8 + ,4 + ,0 + ,0 + ,16 + ,153 + ,305.43 + ,127.72 + ,108.1 + ,-2 + ,4 + ,-2 + ,-2 + ,2 + ,0 + ,10 + ,11 + ,36 + ,20 + ,22 + ,100 + ,341.47 + ,134.88 + ,92.82 + ,4 + ,4 + ,4 + ,-2 + ,4 + ,1 + ,0 + ,0 + ,71 + ,52 + ,8 + ,70 + ,154.54 + ,76.36 + ,78.28 + ,4 + ,4 + ,-2 + ,-2 + ,10 + ,9 + ,2 + ,5 + ,18 + ,14 + ,3 + ,44 + ,250.59 + ,75.61 + ,56.84 + ,4 + ,4 + ,4 + ,4 + ,4 + ,3 + ,6 + ,0 + ,48 + ,47 + ,12 + ,104 + ,449.02 + ,145.53 + ,74.74 + ,4 + ,4 + ,-2 + ,4 + ,7 + ,3 + ,9 + ,4 + ,62 + ,54 + ,22 + ,164 + ,649.48 + ,223.29 + ,153.23 + ,4 + ,4 + ,-2 + ,4 + ,9 + ,7 + ,12 + ,7 + ,42 + ,24 + ,9 + ,80 + ,260.15 + ,98.2 + ,51.23 + ,4 + ,4 + ,4 + ,4 + ,8 + ,4 + ,8 + ,0 + ,35 + ,22 + ,7 + ,165 + ,210.05 + ,86.05 + ,58.58 + ,4 + ,4 + ,4 + ,4 + ,6 + ,2 + ,8 + ,0 + ,69 + ,57 + ,12 + ,171 + ,372.24 + ,131.36 + ,110.18 + ,4 + ,4 + ,4 + ,4 + ,10 + ,10 + ,13 + ,13 + ,70 + ,38 + ,26 + ,157 + ,388.17 + ,127.1 + ,70.01 + ,0 + ,4 + ,4 + ,4 + ,10 + ,8 + ,13 + ,12 + ,69 + ,45 + ,11 + ,123 + ,181.49 + ,88.64 + ,58.72 + ,-2 + ,-2 + ,4 + ,0 + ,3 + ,2 + ,7 + ,2 + ,26 + ,19 + ,14 + ,169 + ,270.36 + ,104.02 + ,51.9 + ,4 + ,4 + ,-2 + ,-2 + ,10 + ,2 + ,8 + ,0 + ,56 + ,45 + ,12 + ,164 + ,309.63 + ,124.18 + ,111.09 + ,4 + ,-2 + ,4 + ,4 + ,10 + ,8 + ,13 + ,10 + ,40 + ,39 + ,23 + ,154 + ,325.01 + ,111.63 + ,62.88 + ,4 + ,0 + ,4 + ,4 + ,5 + ,6 + ,10 + ,0 + ,47 + ,39 + ,5 + ,154 + ,195.6 + ,91.82 + ,77.8 + ,4 + ,4 + ,4 + ,4 + ,10 + ,8 + ,13 + ,13 + ,65 + ,45 + ,15 + ,209 + ,272.61 + ,112.83 + ,122.3 + ,4 + ,4 + ,-2 + ,4 + ,6 + ,4 + ,8 + ,10 + ,48 + ,47 + ,18 + ,139 + ,187.55 + ,87.66 + ,66.52 + ,4 + ,4 + ,4 + ,4 + ,3 + ,0 + ,3 + ,4 + ,69 + ,51 + ,15 + ,103 + ,693.71 + ,171.29 + ,113.07 + ,4 + ,-2 + ,-2 + ,4 + ,8 + ,4 + ,12 + ,0) + ,dim=c(15 + ,125) + ,dimnames=list(c('x$numcomp' + ,'x$numreceived' + ,'x$mrt' + ,'x$nfm' + ,'x$afl' + ,'x$lpm' + ,'x$lpc' + ,'x$Q5' + ,'x$Q9' + ,'x$Q2' + ,'x$Q4' + ,'x$Q1_1' + ,'x$Q2_2' + ,'x$Q1_3' + ,'x$Q2_3') + ,1:125)) > y <- array(NA,dim=c(15,125),dimnames=list(c('x$numcomp','x$numreceived','x$mrt','x$nfm','x$afl','x$lpm','x$lpc','x$Q5','x$Q9','x$Q2','x$Q4','x$Q1_1','x$Q2_2','x$Q1_3','x$Q2_3'),1:125)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par18 = '' > par17 = '' > par16 = '0 1 0 0 1 0' > par15 = '0 1 0 1 0 0' > par14 = '0 1 1 0 0 0' > par13 = '0 1 0 0 0 0' > par12 = '1 0 0 0 0 0' > par11 = '0 0 0 0 0 0' > par10 = '' > par9 = '' > par8 = '14 15' > par7 = '12 13' > par6 = '10 11' > par5 = '8 9' > par4 = '3 4 5 6 7' > par3 = '1 2' > par2 = 'A A B B B B' > par1 = 'COMPUTE REVIEW EXAM0 EXAM1 EXAM2 EXAM3' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Dr. Ian E. Holliday > #To cite this work: Ian E. Holliday, 2009, YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: > #Technical description: > library(plspm) Loading required package: amap Loading required package: diagram Loading required package: shape > library(diagram) > y <- as.data.frame(t(y)) > is.data.frame(y) [1] TRUE > head(y) x$numcomp x$numreceived x$mrt x$nfm x$afl x$lpm x$lpc x$Q5 x$Q9 x$Q2 x$Q4 1 56 34 5 155 124.90 65.55 44.07 4 -2 -2 -2 2 70 40 16 137 550.95 191.95 113.60 4 4 4 -2 3 47 28 15 72 431.64 165.86 112.46 4 4 -2 4 4 12 10 0 0 0.00 0.00 0.00 -2 -2 -2 -2 5 31 20 3 75 111.97 62.78 72.72 0 -2 4 4 6 46 41 2 160 222.41 94.36 70.55 4 0 -2 4 x$Q1_1 x$Q2_2 x$Q1_3 x$Q2_3 1 5 8 11 3 2 8 2 11 3 3 4 4 7 6 4 8 0 12 8 5 6 10 7 6 6 10 8 12 12 > trim <- function(char) { + return(sub('s+$', '', sub('^s+', '', char))) + } > (latnames <- strsplit(par1,' ')[[1]]) [1] "COMPUTE" "REVIEW" "EXAM0" "EXAM1" "EXAM2" "EXAM3" > (n <- length(latnames)) [1] 6 > (L1 <- as.numeric(strsplit(par3,' ')[[1]])) [1] 1 2 > (L2 <- as.numeric(strsplit(par4,' ')[[1]])) [1] 3 4 5 6 7 > (L3 <- as.numeric(strsplit(par5,' ')[[1]])) [1] 8 9 > (L4 <- as.numeric(strsplit(par6,' ')[[1]])) [1] 10 11 > (L5 <- as.numeric(strsplit(par7,' ')[[1]])) [1] 12 13 > (L6 <- as.numeric(strsplit(par8,' ')[[1]])) [1] 14 15 > (L7 <- as.numeric(strsplit(par9,' ')[[1]])) numeric(0) > (L8 <- as.numeric(strsplit(par10,' ')[[1]])) numeric(0) > (S1 <- as.numeric(strsplit(par11,' ')[[1]])) [1] 0 0 0 0 0 0 > (S2 <- as.numeric(strsplit(par12,' ')[[1]])) [1] 1 0 0 0 0 0 > (S3 <- as.numeric(strsplit(par13,' ')[[1]])) [1] 0 1 0 0 0 0 > (S4 <- as.numeric(strsplit(par14,' ')[[1]])) [1] 0 1 1 0 0 0 > (S5 <- as.numeric(strsplit(par15,' ')[[1]])) [1] 0 1 0 1 0 0 > (S6 <- as.numeric(strsplit(par16,' ')[[1]])) [1] 0 1 0 0 1 0 > (S7 <- as.numeric(strsplit(par17,' ')[[1]])) numeric(0) > (S8 <- as.numeric(strsplit(par18,' ')[[1]])) numeric(0) > if (n==1) sat.mat <- rbind(S1) > if (n==2) sat.mat <- rbind(S1,S2) > if (n==3) sat.mat <- rbind(S1,S2,S3) > if (n==4) sat.mat <- rbind(S1,S2,S3,S4) > if (n==5) sat.mat <- rbind(S1,S2,S3,S4,S5) > if (n==6) sat.mat <- rbind(S1,S2,S3,S4,S5,S6) > if (n==7) sat.mat <- rbind(S1,S2,S3,S4,S5,S6,S7) > if (n==8) sat.mat <- rbind(S1,S2,S3,S4,S5,S6,S7,S8) > sat.mat [,1] [,2] [,3] [,4] [,5] [,6] S1 0 0 0 0 0 0 S2 1 0 0 0 0 0 S3 0 1 0 0 0 0 S4 0 1 1 0 0 0 S5 0 1 0 1 0 0 S6 0 1 0 0 1 0 > if (n==1) sat.sets <- list(L1) > if (n==2) sat.sets <- list(L1,L2) > if (n==3) sat.sets <- list(L1,L2,L3) > if (n==4) sat.sets <- list(L1,L2,L3,L4) > if (n==5) sat.sets <- list(L1,L2,L3,L4,L5) > if (n==6) sat.sets <- list(L1,L2,L3,L4,L5,L6) > if (n==7) sat.sets <- list(L1,L2,L3,L4,L5,L6,L7) > if (n==8) sat.sets <- list(L1,L2,L3,L4,L5,L6,L7,L8) > sat.sets [[1]] [1] 1 2 [[2]] [1] 3 4 5 6 7 [[3]] [1] 8 9 [[4]] [1] 10 11 [[5]] [1] 12 13 [[6]] [1] 14 15 > (sat.mod <- strsplit(par2,' ')[[1]]) [1] "A" "A" "B" "B" "B" "B" > res <- plspm(x=y, sat.mat, sat.sets, sat.mod, scheme='centroid', scaled=TRUE, boot.val=TRUE) > (r <- summary(res)) PARTIAL LEAST SQUARES PATH MODELING (PLS-PM) ---------------------------------------------------------- MODEL SPECIFICATION 1 Number of Cases 125 2 Latent Variables 6 3 Manifest Variables 15 4 Scale of Data Standardized Data 5 Weighting Scheme centroid 6 Tolerance Crit 1e-05 7 Max Num Iters 100 8 Convergence Iters 5 9 Paths by PLS-R FALSE 10 Bootstrapping TRUE 11 Bootstrap samples 100 ---------------------------------------------------------- BLOCKS DEFINITION Block Type NMVs Mode 1 S1 Exogenous 2 Reflective 2 S2 Endogenous 5 Reflective 3 S3 Endogenous 2 Formative 4 S4 Endogenous 2 Formative 5 S5 Endogenous 2 Formative 6 S6 Endogenous 2 Formative ---------------------------------------------------------- BLOCKS UNIDIMENSIONALITY Type.measure MVs C.alpha DG.rho eig.1st eig.2nd S1 Reflective 2 0.941 0.971 1.89 0.112 S2 Reflective 5 0.862 0.905 3.33 0.860 S3 Formative 2 0.000 0.000 1.18 0.815 S4 Formative 2 0.000 0.000 1.19 0.813 S5 Formative 2 0.000 0.000 1.30 0.703 S6 Formative 2 0.000 0.000 1.46 0.543 ---------------------------------------------------------- OUTER MODEL weights std.loads communal redundan S1 x$numcomp 0.554 0.976 0.952 0.0000 x$numreceived 0.475 0.967 0.935 0.0000 S2 x$mrt 0.133 0.644 0.415 0.0590 x$nfm 0.333 0.650 0.423 0.0601 x$afl 0.248 0.894 0.800 0.1137 x$lpm 0.262 0.937 0.878 0.1248 x$lpc 0.267 0.861 0.741 0.1053 S3 x$Q5 0.560 0.695 0.483 0.0185 x$Q9 0.731 0.835 0.697 0.0267 S4 x$Q2 0.699 0.810 0.657 0.0343 x$Q4 0.597 0.727 0.529 0.0276 S5 x$Q1_1 0.842 0.945 0.892 0.0759 x$Q2_2 0.344 0.594 0.353 0.0301 S6 x$Q1_3 0.588 0.855 0.730 0.2053 x$Q2_3 0.584 0.853 0.727 0.2043 ---------------------------------------------------------- CORRELATIONS BETWEEN MVs AND LVs S1 S2 S3 S4 S5 S6 S1 x$numcomp 0.9759 0.392 0.1435 0.0195 0.3231 0.3105 x$numreceived 0.9671 0.336 0.1227 -0.0257 0.2924 0.3123 S2 x$mrt 0.2134 0.644 0.1178 0.0810 0.0604 0.0976 x$nfm 0.3929 0.650 0.2333 0.1937 0.2164 0.3864 x$afl 0.2592 0.894 0.1118 0.1002 0.2623 0.3292 x$lpm 0.2831 0.937 0.1341 0.1151 0.2452 0.3434 x$lpc 0.2967 0.861 0.1484 0.1726 0.2324 0.2929 S3 x$Q5 0.1034 0.108 0.6953 0.1510 0.1581 0.0477 x$Q9 0.1091 0.185 0.8348 0.1264 0.0620 0.1411 S4 x$Q2 -0.0032 0.145 0.1074 0.8104 -0.0793 0.0082 x$Q4 0.0014 0.126 0.1710 0.7273 -0.0023 0.1209 S5 x$Q1_1 0.2610 0.264 0.0866 -0.0972 0.9446 0.3740 x$Q2_2 0.2851 0.143 0.1773 0.0730 0.5944 0.3936 S6 x$Q1_3 0.2918 0.301 0.0766 -0.0083 0.4191 0.8546 x$Q2_3 0.2550 0.368 0.1454 0.1417 0.3494 0.8525 ---------------------------------------------------------- INNER MODEL $S2 concept value 1 R2 0.1421 2 Intercept 0.0000 3 path_S1 0.3770 $S3 concept value 1 R2 0.0384 2 Intercept 0.0000 3 path_S2 0.1959 $S4 concept value 1 R2 0.0522 2 Intercept 0.0000 3 path_S2 0.1474 4 path_S3 0.1482 $S5 concept value 1 R2 0.0851 2 Intercept 0.0000 3 path_S2 0.2906 4 path_S4 -0.1081 $S6 concept value 1 R2 0.2811 2 Intercept 0.0000 3 path_S2 0.2907 4 path_S5 0.3714 ---------------------------------------------------------- CORRELATIONS BETWEEN LVs S1 S2 S3 S4 S5 S6 S1 1.0000 0.377 0.138 -0.0014 0.3179 0.3203 S2 0.3770 1.000 0.196 0.1764 0.2716 0.3916 S3 0.1377 0.196 1.000 0.1770 0.1339 0.1299 S4 -0.0014 0.176 0.177 1.0000 -0.0568 0.0779 S5 0.3179 0.272 0.134 -0.0568 1.0000 0.4504 S6 0.3203 0.392 0.130 0.0779 0.4504 1.0000 ---------------------------------------------------------- SUMMARY INNER MODEL LV.Type Measure MVs R.square Av.Commu Av.Redun AVE S1 Exogen Rflct 2 0.0000 0.944 0.0000 0.944 S2 Endogen Rflct 5 0.1421 0.651 0.0926 0.651 S3 Endogen Frmtv 2 0.0384 0.590 0.0226 0.000 S4 Endogen Frmtv 2 0.0522 0.593 0.0310 0.000 S5 Endogen Frmtv 2 0.0851 0.623 0.0530 0.000 S6 Endogen Frmtv 2 0.2811 0.729 0.2048 0.000 ---------------------------------------------------------- GOODNESS-OF-FIT GoF value 1 Absolute 0.2856 2 Relative 0.8180 3 Outer.mod 0.9944 4 Inner.mod 0.8226 ---------------------------------------------------------- TOTAL EFFECTS relationships dir.effects ind.effects tot.effects 1 S1->S2 0.377 0.00000 0.37698 2 S1->S3 0.000 0.07383 0.07383 3 S1->S4 0.000 0.06649 0.06649 4 S1->S5 0.000 0.10238 0.10238 5 S1->S6 0.000 0.14762 0.14762 6 S2->S3 0.196 0.00000 0.19585 7 S2->S4 0.147 0.02902 0.17638 8 S2->S5 0.291 -0.01906 0.27158 9 S2->S6 0.291 0.10087 0.39160 10 S3->S4 0.148 0.00000 0.14819 11 S3->S5 0.000 -0.01601 -0.01601 12 S3->S6 0.000 -0.00595 -0.00595 13 S4->S5 -0.108 0.00000 -0.10805 14 S4->S6 0.000 -0.04013 -0.04013 15 S5->S6 0.371 0.00000 0.37141 --------------------------------------------------------- BOOTSTRAP VALIDATION weights Original Mean.Boot Std.Error perc.05 perc.95 x$numcomp 0.554 0.553 0.0270 0.51438 0.598 x$numreceived 0.475 0.475 0.0266 0.43235 0.515 x$mrt 0.133 0.122 0.0677 -0.00617 0.214 x$nfm 0.333 0.330 0.0608 0.25184 0.440 x$afl 0.249 0.251 0.0274 0.20638 0.293 x$lpm 0.262 0.266 0.0179 0.23873 0.299 x$lpc 0.267 0.263 0.0444 0.19422 0.327 x$Q5 0.560 0.440 0.4099 -0.34797 0.946 x$Q9 0.731 0.673 0.3378 0.06778 1.007 x$Q2 0.699 0.526 0.3957 -0.27714 1.003 x$Q4 0.596 0.591 0.3818 -0.04943 1.019 x$Q1_1 0.842 0.729 0.2022 0.35023 0.999 x$Q2_2 0.344 0.449 0.2682 0.00188 0.833 x$Q1_3 0.588 0.561 0.1866 0.27121 0.916 x$Q2_3 0.584 0.588 0.2065 0.16170 0.853 loadings Original Mean.Boot Std.Error perc.05 perc.95 x$numcomp 0.976 0.976 0.00465 0.9682 0.983 x$numreceived 0.967 0.968 0.00676 0.9555 0.977 x$mrt 0.644 0.650 0.10631 0.4506 0.773 x$nfm 0.650 0.651 0.06605 0.5379 0.748 x$afl 0.894 0.893 0.02658 0.8375 0.924 x$lpm 0.937 0.937 0.01791 0.9053 0.959 x$lpc 0.861 0.846 0.05065 0.7473 0.915 x$Q5 0.695 0.560 0.37802 -0.1682 0.954 x$Q9 0.835 0.762 0.27821 0.3554 0.998 x$Q2 0.810 0.640 0.34526 -0.0256 0.999 x$Q4 0.727 0.685 0.34145 0.1044 0.994 x$Q1_1 0.945 0.859 0.12575 0.6493 0.997 x$Q2_2 0.594 0.660 0.21321 0.2781 0.943 x$Q1_3 0.855 0.825 0.10468 0.6319 0.978 x$Q2_3 0.852 0.840 0.12807 0.5796 0.972 paths Original Mean.Boot Std.Error perc.05 perc.95 S1->S2 0.377 0.3856 0.0669 0.28234 0.507 S2->S3 0.196 0.1727 0.1020 -0.02593 0.313 S2->S4 0.147 0.1445 0.0675 0.04471 0.241 S2->S5 0.291 0.2779 0.0818 0.13883 0.414 S2->S6 0.291 0.2953 0.0717 0.18501 0.404 S3->S4 0.148 0.1652 0.1140 -0.00154 0.321 S4->S5 -0.108 -0.0734 0.1585 -0.29170 0.163 S5->S6 0.371 0.3847 0.0880 0.24303 0.535 rsq Original Mean.Boot Std.Error perc.05 perc.95 S2 0.1421 0.1531 0.0522 0.07976 0.2574 S3 0.0384 0.0401 0.0303 0.00214 0.0982 S4 0.0522 0.0744 0.0413 0.02452 0.1470 S5 0.0851 0.1047 0.0437 0.03418 0.1873 S6 0.2811 0.3074 0.0797 0.18351 0.4425 total.efs Original Mean.Boot Std.Error perc.05 perc.95 S1->S2 0.37698 0.38556 0.0669 0.28234 0.5073 S1->S3 0.07383 0.06898 0.0446 -0.00986 0.1377 S1->S4 0.06649 0.06814 0.0330 0.01748 0.1203 S1->S5 0.10238 0.10265 0.0392 0.04530 0.1710 S1->S6 0.14762 0.15247 0.0370 0.09488 0.2054 S2->S3 0.19585 0.17274 0.1020 -0.02593 0.3134 S2->S4 0.17638 0.17431 0.0742 0.04976 0.2848 S2->S5 0.27158 0.26260 0.0821 0.13432 0.3922 S2->S6 0.39160 0.39736 0.0767 0.27498 0.5069 S3->S4 0.14819 0.16517 0.1140 -0.00154 0.3206 S3->S5 -0.01601 -0.00885 0.0307 -0.05638 0.0365 S3->S6 -0.00595 -0.00278 0.0120 -0.02098 0.0154 S4->S5 -0.10805 -0.07336 0.1585 -0.29170 0.1626 S4->S6 -0.04013 -0.02556 0.0623 -0.11189 0.0685 S5->S6 0.37141 0.38474 0.0880 0.24303 0.5349 > myr <- res$path.coefs > myind <- 1 > for (j in 1:(length(sat.mat[1,])-1)) { + for (i in 1:length(sat.mat[,1])) { + if (sat.mat[i,j] == 1) { + if (res$boot$path[myind,'perc.05'] < 0) { + myr[i,j] = 0 + } + myind = myind + 1 + } + } + } > postscript(file="/var/www/html/rcomp/tmp/1vwyt1283695593.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plotmat(round(myr,4), pos = NULL, curve = 0, name = rownames(res$path.coefs), + lwd = 1, box.lwd = 1, cex.txt = 1, box.type = 'circle', + box.prop = 0.5, box.cex = 1, arr.type = 'triangle', + arr.pos = 0.5, shadow.size = 0.01, prefix = '', arr.lcol = 'blue', + arr.col = 'blue', arr.width = 0.2, main = c('Inner Model', + 'Path Coefficients')) > dev.off() null device 1 > myr <- res$path.coefs > myind <- 1 > myi <- 1 > for (j in 1:(length(sat.mat[1,])-1)) { + for (i in 1:length(sat.mat[,1])) { + if (i > j) { + myr[i,j] = res$boot$total.efs[myi,'Original'] + myi = myi + 1 + if (res$boot$total.efs[myind,'perc.05'] < 0) { + myr[i,j] = 0 + } + myind = myind + 1 + } + } + } > postscript(file="/var/www/html/rcomp/tmp/2geeh1283695593.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plotmat(round(myr,4), pos = NULL, curve = 0, name = rownames(res$path.coefs), + lwd = 1, box.lwd = 1, cex.txt = 1, box.type = 'circle', + box.prop = 0.5, box.cex = 1, arr.type = 'triangle', + arr.pos = 0.5, shadow.size = 0.01, prefix = '', arr.lcol = 'blue', + arr.col = 'blue', arr.width = 0.2, main = c('Inner Model', + 'Total Effects')) > dev.off() null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'PARTIAL LEAST SQUARES PATH MODELING (PLS-PM)',2,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'MODEL SPECIFICATION',2,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Number of Cases',header=TRUE) > a<-table.element(a,r$xxx$obs) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Latent Variables',header=TRUE) > a<-table.element(a,n) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Manifest Variables',header=TRUE) > a<-table.element(a,length(y[1,])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Scaled?',header=TRUE) > a<-table.element(a,r$xxx$scaled) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Weighting Scheme',header=TRUE) > a<-table.element(a,r$xx$scheme) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Bootstrapping?',header=TRUE) > a<-table.element(a,r$xx$boot.val) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Bootstrap samples',header=TRUE) > a<-table.element(a,r$xx$br) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/35fbt1283695593.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'BLOCKS DEFINITION',4,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Block',header=TRUE) > a<-table.element(a,'Type',header=TRUE) > a<-table.element(a,'NMVs',header=TRUE) > a<-table.element(a,'Mode',header=TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,latnames[i],header=TRUE) + a<-table.element(a,r$input$Type[i]) + a<-table.element(a,r$unidim$MVs[i]) + a<-table.element(a,r$unidim$Type.measure[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/4qysy1283695593.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'BLOCKS UNIDIMENSIONALITY',7,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Block',header=TRUE) > a<-table.element(a,'Type.measure',header=TRUE) > a<-table.element(a,'MVs',header=TRUE) > a<-table.element(a,'eig.1st',header=TRUE) > a<-table.element(a,'eig.2nd',header=TRUE) > a<-table.element(a,'C.alpha',header=TRUE) > a<-table.element(a,'DG.rho',header=TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,latnames[i],header=TRUE) + a<-table.element(a,r$unidim$Type.measure[i]) + a<-table.element(a,r$unidim$MVs[i]) + a<-table.element(a,r$unidim$eig.1st[i]) + a<-table.element(a,r$unidim$eig.2nd[i]) + a<-table.element(a,r$unidim$C.alpha[i]) + a<-table.element(a,r$unidim$DG.rho[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/5ugq41283695593.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'OUTER MODEL',5,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Block',header=TRUE) > a<-table.element(a,'weights',header=TRUE) > a<-table.element(a,'std.loads',header=TRUE) > a<-table.element(a,'communal',header=TRUE) > a<-table.element(a,'redundan',header=TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,latnames[i],5,header=TRUE) + a<-table.row.end(a) + for (j in 1:length(r$outer.mod[[i]][,1])) { + a<-table.row.start(a) + a<-table.element(a,rownames(r$outer.mod[[i]])[j],header=T) + a<-table.element(a,r$outer.mod[[i]][j,1]) + a<-table.element(a,r$outer.mod[[i]][j,2]) + a<-table.element(a,r$outer.mod[[i]][j,3]) + a<-table.element(a,r$outer.mod[[i]][j,4]) + a<-table.row.end(a) + } + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/6xzps1283695593.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'CORRELATIONS BETWEEN MVs AND LVs',n+1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Block',header=TRUE) > for (iii in 1:n) { + a<-table.element(a,latnames[iii],header=TRUE) + } > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,latnames[i],n+1,header=TRUE) + a<-table.row.end(a) + for (j in 1:length(r$outer.cor[[i]][,1])) { + a<-table.row.start(a) + a<-table.element(a,rownames(r$outer.cor[[i]])[j],header=T) + for (iii in 1:n) { + a<-table.element(a,r$outer.cor[[i]][j,iii]) + } + a<-table.row.end(a) + } + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/70z5y1283695593.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'INNER MODEL',3,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Block',header=TRUE) > a<-table.element(a,'Concept',header=TRUE) > a<-table.element(a,'Value',header=TRUE) > a<-table.row.end(a) > for (i in 1:(length(labels(r$inner.mod)))) { + a<-table.row.start(a) + print (paste('i=',i,sep='')) + a<-table.element(a,labels(r$inner.mod)[i],3,header=TRUE) + a<-table.row.end(a) + for (j in 1:length(r$inner.mod[[i]][,1])) { + print (paste('j=',j,sep='')) + a<-table.row.start(a) + a<-table.element(a,rownames(r$inner.mod[[i]])[j],header=T) + a<-table.element(a,r$inner.mod[[i]][j,1],header=T) + a<-table.element(a,r$inner.mod[[i]][j,2]) + a<-table.row.end(a) + } + } [1] "i=1" [1] "j=1" [1] "j=2" [1] "j=3" [1] "i=2" [1] "j=1" [1] "j=2" [1] "j=3" [1] "i=3" [1] "j=1" [1] "j=2" [1] "j=3" [1] "j=4" [1] "i=4" [1] "j=1" [1] "j=2" [1] "j=3" [1] "j=4" [1] "i=5" [1] "j=1" [1] "j=2" [1] "j=3" [1] "j=4" > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/8m0441283695593.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'CORRELATIONS BETWEEN LVs',n+1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'',header=TRUE) > for (iii in 1:n) { + a<-table.element(a,latnames[iii],header=TRUE) + } > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,latnames[i],header=T) + for (j in 1:n) { + a<-table.element(a,r$latent.cor[i,j]) + } + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/97jls1283695593.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'SUMMARY INNER MODEL',8,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'',header=TRUE) > a<-table.element(a,'LV.Type',header=TRUE) > a<-table.element(a,'Measure',header=TRUE) > a<-table.element(a,'MVs',header=TRUE) > a<-table.element(a,'R.square',header=TRUE) > a<-table.element(a,'Av.Commu',header=TRUE) > a<-table.element(a,'Av.Redun',header=TRUE) > a<-table.element(a,'AVE',header=TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,latnames[i],header=T) + a<-table.element(a,r$inner.sum[i,1]) + a<-table.element(a,r$inner.sum[i,2]) + a<-table.element(a,r$inner.sum[i,3]) + a<-table.element(a,r$inner.sum[i,4]) + a<-table.element(a,r$inner.sum[i,5]) + a<-table.element(a,r$inner.sum[i,6]) + a<-table.element(a,r$inner.sum[i,7]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/10t1jy1283695593.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'GOODNESS-OF-FIT',2,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'GoF',header=TRUE) > a<-table.element(a,'Value',header=TRUE) > a<-table.row.end(a) > for (i in 1:4) { + a<-table.row.start(a) + a<-table.element(a,r$gof[i,1],header=T) + a<-table.element(a,r$gof[i,2]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/11maij1283695593.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'TOTAL EFFECTS',4,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'relationships',header=TRUE) > a<-table.element(a,'dir.effect',header=TRUE) > a<-table.element(a,'ind.effect',header=TRUE) > a<-table.element(a,'tot.effect',header=TRUE) > a<-table.row.end(a) > for (i in 1:length(r$effects[,1])) { + a<-table.row.start(a) + a<-table.element(a,r$effects[i,1],header=T) + a<-table.element(a,r$effects[i,2]) + a<-table.element(a,r$effects[i,3]) + a<-table.element(a,r$effects[i,4]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/12ptzo1283695593.tab") > dum <- r$boot$weights > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'BOOTSTRAP VALIDATION - WEIGHTS',length(colnames(dum))+1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'',header=TRUE) > for (i in 1:length(colnames(dum))) { + a<-table.element(a,colnames(dum)[i],header=TRUE) + } > a<-table.row.end(a) > for (i in 1:length(rownames(dum))) { + a<-table.row.start(a) + a<-table.element(a,rownames(dum)[i],header=T) + for (j in 1:length(colnames(dum))) { + a<-table.element(a,dum[i,j]) + } + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/13abxu1283695593.tab") > dum <- r$boot$loadings > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'BOOTSTRAP VALIDATION - LOADINGS',length(colnames(dum))+1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'',header=TRUE) > for (i in 1:length(colnames(dum))) { + a<-table.element(a,colnames(dum)[i],header=TRUE) + } > a<-table.row.end(a) > for (i in 1:length(rownames(dum))) { + a<-table.row.start(a) + a<-table.element(a,rownames(dum)[i],header=T) + for (j in 1:length(colnames(dum))) { + a<-table.element(a,dum[i,j]) + } + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/14l3ff1283695593.tab") > dum <- r$boot$paths > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'BOOTSTRAP VALIDATION - PATHS',length(colnames(dum))+1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'',header=TRUE) > for (i in 1:length(colnames(dum))) { + a<-table.element(a,colnames(dum)[i],header=TRUE) + } > a<-table.row.end(a) > for (i in 1:length(rownames(dum))) { + a<-table.row.start(a) + a<-table.element(a,rownames(dum)[i],header=T) + for (j in 1:length(colnames(dum))) { + a<-table.element(a,dum[i,j]) + } + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/15oldl1283695593.tab") > dum <- r$boot$rsq > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'BOOTSTRAP VALIDATION - RSQ',length(colnames(dum))+1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'',header=TRUE) > for (i in 1:length(colnames(dum))) { + a<-table.element(a,colnames(dum)[i],header=TRUE) + } > a<-table.row.end(a) > for (i in 1:length(rownames(dum))) { + a<-table.row.start(a) + a<-table.element(a,rownames(dum)[i],header=T) + for (j in 1:length(colnames(dum))) { + a<-table.element(a,dum[i,j]) + } + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/16zuu61283695593.tab") > dum <- r$boot$total.efs > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'BOOTSTRAP VALIDATION - TOTAL EFFECTS',length(colnames(dum))+1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'',header=TRUE) > for (i in 1:length(colnames(dum))) { + a<-table.element(a,colnames(dum)[i],header=TRUE) + } > a<-table.row.end(a) > for (i in 1:length(rownames(dum))) { + a<-table.row.start(a) + a<-table.element(a,rownames(dum)[i],header=T) + for (j in 1:length(colnames(dum))) { + a<-table.element(a,dum[i,j]) + } + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/17kvtc1283695593.tab") > try(system("convert tmp/1vwyt1283695593.ps tmp/1vwyt1283695593.png",intern=TRUE)) character(0) > try(system("convert tmp/2geeh1283695593.ps tmp/2geeh1283695593.png",intern=TRUE)) character(0) > > #-SERVER-193.190.124.10:1001 > > > > proc.time() user system elapsed 15.369 0.278 15.688