R version 2.6.1 (2007-11-26) Copyright (C) 2007 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 <- c(33.02 + ,22.69 + ,24.54 + ,23.56 + ,23.56 + ,22.83 + ,23.35 + ,23.12 + ,23.44 + ,22.37 + ,25.28 + ,27.80 + ,24.24 + ,23.08 + ,27.79 + ,28.00 + ,32.73 + ,38.45 + ,40.39 + ,29.76 + ,23.75 + ,23.25 + ,23.09 + ,23.82 + ,26.58 + ,23.05 + ,21.07 + ,22.63 + ,23.58 + ,28.57 + ,34.14 + ,24.92 + ,27.25 + ,23.40 + ,25.26 + ,26.81 + ,27.11 + ,34.33 + ,26.47 + ,27.25 + ,30.02 + ,29.35 + ,27.75 + ,23.69 + ,26.94 + ,49.43 + ,26.31 + ,37.37 + ,31.05 + ,23.25 + ,25.36 + ,28.40 + ,23.33 + ,27.90 + ,35.25 + ,56.71 + ,23.51 + ,22.43 + ,29.86 + ,23.37 + ,51.48 + ,54.06 + ,59.83 + ,57.61 + ,68.69 + ,64.71 + ,67.66 + ,74.05 + ,79.17 + ,80.95 + ,79.87 + ,75.87 + ,76.70 + ,75.19 + ,79.22 + ,157.44 + ,75.95 + ,78.05 + ,73.46 + ,71.45 + ,74.47 + ,69.40 + ,48.26 + ,32.60 + ,55.52 + ,62.77 + ,58.41 + ,67.76 + ,71.95 + ,60.07 + ,48.67 + ,55.05 + ,54.77 + ,44.14 + ,49.53 + ,47.04 + ,46.29 + ,46.76 + ,39.34 + ,34.29 + ,33.79 + ,31.54 + ,29.49 + ,31.06 + ,34.10 + ,30.28 + ,29.72 + ,40.37 + ,29.89 + ,29.77 + ,26.16 + ,28.51 + ,25.44 + ,24.05 + ,25.33 + ,24.23 + ,24.20 + ,24.03 + ,31.21 + ,25.54 + ,28.84 + ,27.27 + ,27.77 + ,26.33 + ,23.78 + ,25.75 + ,26.31 + ,26.29 + ,30.44 + ,41.44 + ,62.59 + ,60.95 + ,40.48 + ,27.29 + ,31.17 + ,24.72 + ,26.24 + ,28.55 + ,22.83 + ,24.00 + ,24.18 + ,22.97 + ,23.68 + ,23.90 + ,25.67 + ,21.30 + ,24.94 + ,25.35 + ,21.25 + ,31.30 + ,23.74 + ,22.69 + ,22.67 + ,24.14 + ,23.72 + ,24.15 + ,22.41 + ,32.10 + ,25.77 + ,24.18 + ,25.52 + ,24.11 + ,24.59 + ,23.68 + ,27.12 + ,24.12 + ,24.16 + ,35.98 + ,29.45 + ,24.48 + ,22.71 + ,25.93 + ,23.31 + ,23.03 + ,24.12 + ,25.01 + ,23.52 + ,24.53 + ,30.98 + ,23.25 + ,23.45 + ,24.45 + ,24.67 + ,24.37 + ,23.82 + ,24.71 + ,24.85 + ,25.19 + ,23.24 + ,23.20 + ,25.08 + ,25.11 + ,24.41 + ,26.35 + ,22.36 + ,28.30 + ,25.40 + ,24.60 + ,24.06 + ,25.41 + ,23.97 + ,23.66 + ,25.64 + ,25.27 + ,27.33 + ,24.05 + ,24.38 + ,26.85 + ,25.47 + ,39.26 + ,25.11 + ,23.78 + ,23.81 + ,23.68 + ,23.73 + ,23.54 + ,23.00 + ,33.34 + ,23.06 + ,24.45 + ,24.80 + ,26.26 + ,23.19 + ,24.05 + ,24.41 + ,26.12 + ,23.10 + ,36.31 + ,28.20 + ,23.96 + ,24.57 + ,24.03 + ,25.29 + ,24.38 + ,23.70 + ,23.82 + ,24.62 + ,24.83 + ,30.52 + ,25.27 + ,26.30 + ,24.03 + ,54.01 + ,27.94 + ,23.45 + ,23.90 + ,25.69 + ,24.58 + ,28.93 + ,24.43 + ,22.08 + ,28.27 + ,23.63 + ,23.79 + ,23.35 + ,109.07 + ,22.73 + ,28.87 + ,23.87 + ,24.73 + ,25.50 + ,23.91 + ,90.53 + ,37.44 + ,26.55 + ,22.55 + ,22.96 + ,23.00 + ,24.67 + ,30.38 + ,33.02 + ,22.59 + ,24.47 + ,23.44 + ,24.08 + ,23.60 + ,25.47 + ,36.44 + ,26.04 + ,25.50 + ,24.33 + ,24.88 + ,25.38 + ,24.30 + ,24.74 + ,43.78 + ,24.93 + ,37.83 + ,27.86 + ,23.82 + ,26.54 + ,25.01 + ,24.23 + ,25.24 + ,28.45 + ,26.91 + ,24.83 + ,27.67 + ,35.54 + ,26.44 + ,90.67 + ,24.22 + ,26.43 + ,26.12 + ,25.28 + ,26.43 + ,26.69 + ,26.02 + ,24.95 + ,24.67 + ,25.47 + ,25.10 + ,24.45 + ,25.55 + ,25.83 + ,23.30 + ,25.41 + ,26.85 + ,23.97 + ,25.73 + ,25.37 + ,26.56 + ,27.92 + ,24.63 + ,26.58 + ,25.12 + ,27.26 + ,27.37 + ,24.81 + ,32.11 + ,27.15 + ,24.81 + ,26.13 + ,24.97 + ,25.32 + ,27.18 + ,24.92 + ,33.41 + ,25.04 + ,23.40 + ,24.95 + ,26.03 + ,24.65 + ,24.88 + ,24.90 + ,39.91 + ,26.50 + ,42.79 + ,28.31 + ,23.10 + ,25.27 + ,26.72 + ,24.93 + ,25.42 + ,24.07 + ,26.92 + ,26.71 + ,26.75 + ,31.35 + ,23.36 + ,28.27 + ,27.33 + ,26.14 + ,32.22 + ,27.90 + ,24.07 + ,28.13 + ,31.32 + ,31.25 + ,26.70 + ,23.68 + ,28.06 + ,25.80 + ,26.18 + ,25.73 + ,33.13 + ,24.58 + ,26.75 + ,29.28 + ,28.97 + ,26.71 + ,25.75 + ,25.23 + ,26.86 + ,28.74 + ,24.48 + ,26.20 + ,25.97 + ,23.07 + ,32.97 + ,37.59 + ,28.01 + ,25.52 + ,26.86 + ,28.98 + ,28.41 + ,25.43 + ,32.20 + ,25.55 + ,27.05 + ,27.23 + ,27.09 + ,27.78 + ,27.01 + ,29.35 + ,36.90 + ,35.84 + ,42.25 + ,32.05 + ,26.03 + ,27.82 + ,25.66 + ,26.16 + ,25.11 + ,29.28 + ,25.03 + ,28.80 + ,37.79 + ,39.84 + ,31.73 + ,38.47 + ,25.49 + ,30.06 + ,28.85 + ,55.95 + ,60.64 + ,66.31 + ,76.91 + ,92.23 + ,63.35 + ,67.86 + ,64.38 + ,47.54 + ,21.94 + ,25.69 + ,95.92 + ,24.01 + ,24.01 + ,23.22 + ,22.83 + ,24.49 + ,25.01 + ,26.80 + ,23.47 + ,27.91 + ,25.94 + ,24.42 + ,24.11 + ,25.08 + ,31.09 + ,43.91 + ,26.33 + ,24.22 + ,23.38 + ,26.18 + ,24.44 + ,24.91 + ,35.21 + ,25.22 + ,24.06 + ,23.70 + ,23.07 + ,25.04 + ,26.13 + ,26.06 + ,38.14 + ,23.89 + ,40.01 + ,30.95 + ,26.24 + ,26.34 + ,25.28 + ,29.11 + ,23.59 + ,23.31 + ,22.34 + ,24.42 + ,24.02 + ,31.98 + ,22.65 + ,40.52 + ,24.23 + ,23.14 + ,23.74 + ,24.75 + ,23.71 + ,23.16 + ,23.92 + ,26.70 + ,25.31 + ,26.39 + ,25.17 + ,24.09 + ,23.24 + ,28.57 + ,27.02 + ,25.52 + ,23.12 + ,24.05 + ,23.32 + ,23.93 + ,23.86 + ,23.33 + ,24.75 + ,26.36 + ,23.55 + ,22.97 + ,24.77 + ,24.34 + ,32.40 + ,38.15 + ,26.37 + ,24.14 + ,23.75 + ,25.38 + ,23.67 + ,26.29 + ,33.24 + ,24.80 + ,27.46 + ,23.48 + ,25.38 + ,24.20 + ,25.94 + ,25.88 + ,45.06 + ,27.68 + ,41.42 + ,29.04 + ,24.75 + ,22.78 + ,24.23 + ,25.06 + ,24.57 + ,24.69 + ,23.29 + ,23.09 + ,24.12 + ,29.50 + ,23.06 + ,29.28 + ,24.11 + ,25.65 + ,25.31 + ,23.20 + ,26.11 + ,25.10 + ,26.38 + ,22.73 + ,22.52 + ,23.81 + ,22.09 + ,24.02 + ,23.92 + ,23.99 + ,35.54 + ,23.62 + ,25.79 + ,23.98 + ,24.52 + ,24.95 + ,24.76 + ,23.32 + ,24.27 + ,29.36 + ,23.77 + ,25.93 + ,26.71 + ,24.17 + ,30.42 + ,37.25 + ,25.03 + ,22.77 + ,24.95 + ,24.87 + ,24.22 + ,23.64 + ,31.30 + ,23.90 + ,22.64 + ,26.09 + ,24.57 + ,24.45 + ,24.54 + ,24.62 + ,34.16 + ,25.13 + ,43.14 + ,29.11 + ,25.56 + ,23.50 + ,22.99 + ,23.60 + ,23.94 + ,24.32 + ,24.09 + ,23.66 + ,22.60 + ,32.42 + ,24.78 + ,34.30 + ,49.68 + ,21.73 + ,20.67 + ,22.99 + ,23.47 + ,22.28 + ,20.89 + ,20.87 + ,22.71 + ,21.02 + ,21.37 + ,23.24 + ,21.31 + ,20.72 + ,118.83 + ,22.51 + ,22.29 + ,22.89 + ,24.11 + ,22.97 + ,22.56 + ,23.75 + ,23.10 + ,23.40 + ,22.09 + ,23.58 + ,22.06 + ,21.63 + ,30.67 + ,26.64 + ,25.96 + ,30.74 + ,27.81 + ,28.83 + ,29.33 + ,28.30 + ,38.93 + ,30.18 + ,29.54 + ,30.33 + ,30.97 + ,27.83 + ,30.45 + ,30.09 + ,44.98 + ,28.21 + ,43.11 + ,33.51 + ,30.59 + ,29.61 + ,30.21 + ,31.14 + ,31.15 + ,30.37 + ,31.11 + ,27.85 + ,29.61 + ,35.47 + ,30.94 + ,47.28 + ,29.61 + ,35.16 + ,28.16 + ,28.90 + ,32.52 + ,32.79 + ,31.14 + ,31.65 + ,30.24 + ,28.98 + ,29.04 + ,30.12 + ,29.59 + ,29.54 + ,50.00 + ,29.51 + ,28.45 + ,29.91 + ,29.39 + ,28.51 + ,27.95 + ,29.27 + ,29.85 + ,28.65 + ,27.68 + ,27.94 + ,27.29 + ,28.48 + ,35.56 + ,46.72 + ,28.60 + ,29.07 + ,30.00 + ,31.48 + ,27.34 + ,29.56 + ,36.96 + ,27.17 + ,28.69 + ,27.52 + ,25.08 + ,25.79 + ,23.41 + ,23.51 + ,36.14 + ,25.11 + ,42.16 + ,25.08 + ,24.84 + ,23.04 + ,23.20 + ,24.85 + ,24.92 + ,25.85 + ,21.47 + ,22.69 + ,24.74 + ,29.32 + ,24.73 + ,48.24 + ,25.58 + ,24.66 + ,24.55 + ,25.12 + ,23.44 + ,26.64 + ,24.96 + ,23.85 + ,23.16 + ,24.00 + ,25.11 + ,22.40 + ,22.00 + ,25.46 + ,36.69 + ,24.06 + ,24.33 + ,23.29 + ,27.14 + ,24.18 + ,23.06 + ,23.02 + ,23.99 + ,25.91 + ,21.24 + ,22.89 + ,23.17 + ,23.49 + ,32.60 + ,38.70 + ,25.01 + ,24.20 + ,24.94 + ,25.41 + ,25.21 + ,22.22 + ,33.11 + ,24.77 + ,24.60 + ,23.99 + ,26.30 + ,25.13 + ,23.96 + ,26.95 + ,46.07 + ,25.32 + ,38.00 + ,27.74 + ,26.24 + ,22.50 + ,25.60 + ,24.27 + ,24.41 + ,24.11 + ,24.27 + ,27.99 + ,29.50 + ,37.08 + ,29.19 + ,33.76 + ,28.66 + ,29.73 + ,31.02 + ,30.64 + ,32.35 + ,34.43 + ,24.98 + ,25.07 + ,23.81 + ,24.70 + ,24.97 + ,25.03 + ,25.92 + ,25.37 + ,127.19 + ,25.67 + ,26.23 + ,24.32 + ,25.55 + ,25.89 + ,24.43 + ,24.85 + ,25.11 + ,25.96 + ,24.72 + ,23.76 + ,23.34 + ,26.15 + ,32.80 + ,41.38 + ,24.62 + ,25.26 + ,25.84 + ,24.96 + ,25.48 + ,24.65 + ,34.63 + ,25.24 + ,25.84 + ,26.69 + ,23.52 + ,25.77 + ,24.43 + ,24.80 + ,39.50 + ,25.09 + ,39.38 + ,29.78 + ,22.47 + ,26.07 + ,24.67 + ,24.10 + ,25.51 + ,25.08 + ,22.38 + ,25.99 + ,25.01 + ,32.52 + ,24.82 + ,44.82 + ,24.75 + ,23.00 + ,24.11 + ,28.28 + ,26.48 + ,24.33 + ,25.23 + ,23.17 + ,22.44 + ,23.85 + ,24.30 + ,25.31 + ,23.28 + ,25.81 + ,36.40 + ,24.40 + ,24.49 + ,25.55 + ,24.65 + ,22.81 + ,23.95 + ,25.26 + ,26.68 + ,24.42 + ,25.05 + ,25.34 + ,23.43 + ,25.30 + ,30.55 + ,27.34 + ,24.75 + ,26.40 + ,24.81 + ,24.03 + ,23.24 + ,24.66 + ,31.27 + ,22.27 + ,24.86 + ,22.78 + ,23.90 + ,24.59 + ,24.15 + ,24.46 + ,35.00 + ,25.12 + ,41.65 + ,24.57 + ,23.47 + ,24.03 + ,23.26 + ,26.07 + ,23.66 + ,25.36 + ,23.49 + ,26.18 + ,25.03 + ,30.25 + ,24.12 + ,38.68 + ,27.19 + ,24.04 + ,23.53 + ,25.48 + ,24.35 + ,25.20 + ,24.57 + ,24.51 + ,24.98 + ,24.80 + ,25.31 + ,24.27 + ,24.81 + ,25.34 + ,27.43 + ,24.78 + ,23.59 + ,24.96 + ,24.02 + ,22.60 + ,23.18 + ,24.13 + ,24.60 + ,26.28 + ,24.18 + ,25.47 + ,24.57 + ,24.28 + ,30.99 + ,26.39 + ,25.65 + ,25.70 + ,25.40 + ,24.31 + ,23.40 + ,26.72 + ,32.78 + ,24.93 + ,24.43 + ,23.41 + ,22.42 + ,24.73 + ,25.93 + ,23.13 + ,35.25 + ,23.38 + ,40.88 + ,28.32 + ,23.57 + ,23.41 + ,24.22 + ,24.48 + ,23.39 + ,23.90 + ,24.06 + ,23.98 + ,23.35 + ,30.08 + ,24.60 + ,40.25 + ,25.28 + ,25.39 + ,27.25 + ,24.18 + ,24.10 + ,24.28 + ,24.64 + ,25.40 + ,23.46 + ,25.95 + ,24.10 + ,23.07 + ,25.35 + ,26.99 + ,119.44 + ,26.21 + ,24.50 + ,24.00 + ,24.87 + ,24.01 + ,25.20 + ,23.71 + ,23.79 + ,24.83 + ,22.76 + ,22.06 + ,23.13 + ,21.35 + ,32.54 + ,39.56 + ,22.65 + ,24.01 + ,24.02 + ,24.31 + ,24.43 + ,25.00 + ,33.02 + ,24.28 + ,24.29 + ,24.41 + ,22.83 + ,24.81 + ,24.86 + ,24.98 + ,32.02 + ,24.22 + ,41.42 + ,25.85 + ,24.48 + ,26.19 + ,26.29 + ,24.46 + ,24.05 + ,24.31 + ,23.51 + ,22.07 + ,24.42 + ,33.90 + ,22.97 + ,36.10 + ,23.40 + ,24.55 + ,23.11 + ,24.12 + ,22.41 + ,22.07 + ,25.95 + ,23.17 + ,23.29 + ,23.21 + ,24.85 + ,22.47 + ,22.69 + ,26.25 + ,44.49 + ,24.01 + ,22.41 + ,24.52 + ,25.98 + ,23.39 + ,22.42 + ,25.22 + ,23.27 + ,24.71 + ,24.37 + ,24.32 + ,24.77 + ,24.21 + ,31.63 + ,33.89 + ,22.99 + ,25.41 + ,28.41 + ,23.66 + ,25.53 + ,26.41 + ,35.66 + ,23.53 + ,26.21 + ,25.59 + ,23.05 + ,24.95 + ,25.74 + ,24.34 + ,41.66 + ,30.53 + ,40.20 + ,29.46 + ,24.17 + ,26.15 + ,25.81 + ,26.15 + ,24.66 + ,22.57 + ,24.41 + ,27.28 + ,26.09 + ,32.05 + ,22.78 + ,48.54 + ,24.95 + ,24.68 + ,24.52 + ,24.70 + ,24.75 + ,24.48 + ,23.85 + ,24.47 + ,23.92 + ,23.34 + ,23.12 + ,23.33 + ,24.15 + ,23.61 + ,24.92 + ,22.85 + ,27.01 + ,24.48 + ,26.61 + ,24.23 + ,23.58 + ,23.91 + ,24.09 + ,27.00 + ,22.73 + ,23.35 + ,22.42 + ,22.39 + ,31.67 + ,35.98 + ,24.01 + ,25.49 + ,23.38 + ,24.52 + ,25.35 + ,23.80 + ,31.94 + ,23.00 + ,23.19 + ,28.10 + ,24.64 + ,22.16 + ,25.43 + ,25.36 + ,29.13 + ,22.24 + ,39.26 + ,27.16 + ,22.89 + ,23.81 + ,25.21 + ,23.92 + ,24.05 + ,24.16 + ,25.57 + ,22.45 + ,23.89 + ,32.89 + ,24.05 + ,24.32 + ,24.15 + ,25.49 + ,23.64 + ,24.60 + ,25.23 + ,25.98 + ,23.49 + ,25.11 + ,24.62 + ,25.46 + ,25.66 + ,26.53 + ,24.62 + ,25.70 + ,112.49 + ,24.27 + ,23.64 + ,22.63 + ,25.09 + ,24.54 + ,24.44 + ,23.37 + ,24.69 + ,28.17 + ,24.81 + ,25.15 + ,26.81 + ,23.69 + ,34.99 + ,42.49 + ,23.80 + ,24.36 + ,24.91 + ,25.47 + ,24.31 + ,24.46 + ,32.81 + ,23.47 + ,25.34 + ,25.02 + ,23.89 + ,23.34 + ,20.79 + ,23.04 + ,43.64 + ,21.54 + ,37.68 + ,28.26 + ,23.74 + ,23.82 + ,24.97 + ,23.59 + ,22.71 + ,24.18 + ,23.51 + ,23.78 + ,22.68 + ,30.98 + ,23.03 + ,37.22 + ,22.62 + ,25.52 + ,25.53 + ,24.80 + ,23.14 + ,23.63 + ,22.73 + ,24.23 + ,23.79 + ,24.16 + ,22.39 + ,24.12 + ,24.41 + ,25.53 + ,28.23 + ,23.88 + ,25.05 + ,25.90 + ,26.73 + ,25.36 + ,24.12 + ,24.60 + ,25.08 + ,26.81 + ,21.93 + ,23.48 + ,25.64 + ,22.44 + ,29.75 + ,37.70 + ,25.14 + ,25.25 + ,24.27 + ,24.96 + ,23.65 + ,22.66 + ,33.33 + ,23.81 + ,23.26 + ,23.82 + ,23.11 + ,22.33 + ,23.37 + ,24.36 + ,36.80 + ,22.77 + ,38.10 + ,27.83 + ,22.39 + ,23.27 + ,25.00 + ,23.32 + ,22.82 + ,23.35 + ,23.78 + ,23.38 + ,24.04 + ,32.13 + ,24.01 + ,24.85 + ,24.42 + ,23.79 + ,23.49 + ,21.75 + ,25.33 + ,22.28 + ,22.64 + ,23.45 + ,23.17 + ,23.81 + ,23.29 + ,25.43 + ,24.29 + ,23.96 + ,26.56 + ,23.70 + ,22.14 + ,23.47 + ,24.77 + ,22.80 + ,23.12 + ,23.26 + ,23.22 + ,24.69 + ,24.90 + ,24.10 + ,22.04 + ,24.24 + ,31.59 + ,32.61 + ,24.75 + ,24.84 + ,24.64 + ,22.64 + ,23.34 + ,23.72 + ,34.58 + ,23.17 + ,22.86 + ,23.96 + ,21.18 + ,25.47 + ,23.22 + ,23.34 + ,35.82 + ,23.42 + ,37.05 + ,24.81 + ,25.93 + ,23.30 + ,22.74 + ,21.89 + ,24.90 + ,23.96 + ,23.43 + ,24.70 + ,25.96 + ,30.91 + ,23.17 + ,35.99 + ,24.62 + ,22.05 + ,23.29 + ,23.40 + ,24.31 + ,22.80 + ,23.09 + ,22.99 + ,21.74 + ,23.21 + ,23.25 + ,23.57 + ,22.42 + ,24.25 + ,119.19 + ,24.32 + ,23.87 + ,24.42 + ,23.47 + ,23.81 + ,24.92 + ,24.01 + ,23.82 + ,23.82 + ,21.51 + ,22.90 + ,21.82 + ,21.00 + ,28.35 + ,30.94 + ,23.03 + ,21.31 + ,21.62 + ,21.71 + ,21.88 + ,21.07 + ,32.17 + ,23.10 + ,23.45 + ,23.54 + ,23.15 + ,23.86 + ,21.97 + ,22.28 + ,35.86 + ,23.09 + ,34.68 + ,28.99 + ,23.19 + ,21.81 + ,25.44 + ,23.31 + ,23.44 + ,22.41 + ,23.39 + ,23.27 + ,21.33 + ,30.12 + ,24.10 + ,24.62 + ,23.40 + ,24.54 + ,24.15 + ,22.76 + ,23.76 + ,24.69 + ,23.65 + ,22.92 + ,21.98 + ,23.79 + ,22.22 + ,24.02 + ,24.04 + ,22.75 + ,24.16 + ,24.50 + ,21.95 + ,23.73 + ,24.44 + ,22.75 + ,21.48 + ,22.84 + ,23.25 + ,23.90 + ,23.23 + ,22.67 + ,22.66 + ,21.37 + ,31.02 + ,24.06 + ,24.25 + ,21.07 + ,23.80 + ,22.09 + ,22.29 + ,23.90 + ,30.94 + ,22.17 + ,22.97 + ,22.52 + ,22.80 + ,22.08 + ,22.63 + ,23.09 + ,23.69 + ,21.95 + ,36.23 + ,25.55 + ,22.54 + ,22.97 + ,22.94 + ,22.90 + ,24.44 + ,23.35 + ,22.77 + ,22.37 + ,23.39 + ,31.38) > x <-sort(x[!is.na(x)]) > q1 <- function(data,n,p,i,f) { + np <- n*p; + i <<- floor(np) + f <<- np - i + qvalue <- (1-f)*data[i] + f*data[i+1] + } > q2 <- function(data,n,p,i,f) { + np <- (n+1)*p + i <<- floor(np) + f <<- np - i + qvalue <- (1-f)*data[i] + f*data[i+1] + } > q3 <- function(data,n,p,i,f) { + np <- n*p + i <<- floor(np) + f <<- np - i + if (f==0) { + qvalue <- data[i] + } else { + qvalue <- data[i+1] + } + } > q4 <- function(data,n,p,i,f) { + np <- n*p + i <<- floor(np) + f <<- np - i + if (f==0) { + qvalue <- (data[i]+data[i+1])/2 + } else { + qvalue <- data[i+1] + } + } > q5 <- function(data,n,p,i,f) { + np <- (n-1)*p + i <<- floor(np) + f <<- np - i + if (f==0) { + qvalue <- data[i+1] + } else { + qvalue <- data[i+1] + f*(data[i+2]-data[i+1]) + } + } > q6 <- function(data,n,p,i,f) { + np <- n*p+0.5 + i <<- floor(np) + f <<- np - i + qvalue <- data[i] + } > q7 <- function(data,n,p,i,f) { + np <- (n+1)*p + i <<- floor(np) + f <<- np - i + if (f==0) { + qvalue <- data[i] + } else { + qvalue <- f*data[i] + (1-f)*data[i+1] + } + } > q8 <- function(data,n,p,i,f) { + np <- (n+1)*p + i <<- floor(np) + f <<- np - i + if (f==0) { + qvalue <- data[i] + } else { + if (f == 0.5) { + qvalue <- (data[i]+data[i+1])/2 + } else { + if (f < 0.5) { + qvalue <- data[i] + } else { + qvalue <- data[i+1] + } + } + } + } > lx <- length(x) > qval <- array(NA,dim=c(99,8)) > mystep <- 25 > mystart <- 25 > if (lx>10){ + mystep=10 + mystart=10 + } > if (lx>20){ + mystep=5 + mystart=5 + } > if (lx>50){ + mystep=2 + mystart=2 + } > if (lx>=100){ + mystep=1 + mystart=1 + } > for (perc in seq(mystart,99,mystep)) { + qval[perc,1] <- q1(x,lx,perc/100,i,f) + qval[perc,2] <- q2(x,lx,perc/100,i,f) + qval[perc,3] <- q3(x,lx,perc/100,i,f) + qval[perc,4] <- q4(x,lx,perc/100,i,f) + qval[perc,5] <- q5(x,lx,perc/100,i,f) + qval[perc,6] <- q6(x,lx,perc/100,i,f) + qval[perc,7] <- q7(x,lx,perc/100,i,f) + qval[perc,8] <- q8(x,lx,perc/100,i,f) + } > postscript(file="/var/www/html/rcomp/tmp/15iav1203954872.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > myqqnorm <- qqnorm(x,col=2) > dev.off() null device 1 > load(file='/var/www/html/rcomp/createtable') > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Percentiles - Ungrouped Data',9,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'p',1,TRUE) > a<-table.element(a,hyperlink('http://www.xycoon.com/method_1.htm', 'Weighted Average at Xnp',''),1,TRUE) > a<-table.element(a,hyperlink('http://www.xycoon.com/method_2.htm','Weighted Average at X(n+1)p',''),1,TRUE) > a<-table.element(a,hyperlink('http://www.xycoon.com/method_3.htm','Empirical Distribution Function',''),1,TRUE) > a<-table.element(a,hyperlink('http://www.xycoon.com/method_4.htm','Empirical Distribution Function - Averaging',''),1,TRUE) > a<-table.element(a,hyperlink('http://www.xycoon.com/method_5.htm','Empirical Distribution Function - Interpolation',''),1,TRUE) > a<-table.element(a,hyperlink('http://www.xycoon.com/method_6.htm','Closest Observation',''),1,TRUE) > a<-table.element(a,hyperlink('http://www.xycoon.com/method_7.htm','True Basic - Statistics Graphics Toolkit',''),1,TRUE) > a<-table.element(a,hyperlink('http://www.xycoon.com/method_8.htm','MS Excel (old versions)',''),1,TRUE) > a<-table.row.end(a) > for (perc in seq(mystart,99,mystep)) { + a<-table.row.start(a) + a<-table.element(a,round(perc/100,2),1,TRUE) + for (j in 1:8) { + a<-table.element(a,round(qval[perc,j],6)) + } + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/2ddrn1203954872.tab") > > system("convert tmp/15iav1203954872.ps tmp/15iav1203954872.png") > > > proc.time() user system elapsed 1.670 0.303 1.696