如何从宽观测子集中获得具有新观测数的正态分布

Question

我有一个变量的以下观察值

vars
  [1] 1053 1018 1048 1040 1046 1038 1014 1020 1038 1006 1016 1025 1011 1018 1014 1021 1022 1000 1003 1018 1037 1053 1033 1047 1018
 [26] 1022 1014 1048 1048 1026 1028 1056 1015 1030 1008 1035 1007 1026 1021 1008 1053 1052 1009 1015 1020 1013 1023 1015 1018 1038
 [51] 1007 1018 1035 1030 1016 1034 1028 1018 1018 1024 1024 1030 1028 1027 1016 1012 1023 1022 1045 1012 1019 1010 1086 1016 1034
 [76] 1043 1027 1013 1022 1030 1022 1059 1021 1008 1023 1010 1016 1020 1062 1013 1064 1052 1018 1033 1000 1025 1011 1016 1017 1048
[101] 1030 1022 1018 1020 1011 1095 1038 1018 1042 1007 1048 1038 1024 1016 1008 1051 1014 1020 1024 1008 1046 1045 1064 1075 1035
[126] 1010 1010 1010 1012 1012 1013 1016 1039 1030 1003 1013 1025 1005 1000 1029 1028 1008 1014 1026 1020 1062 1023 1008 1032 1024
[151] 1026 1011 1018 1028 1031 1042 1057 1054 1010 1012 1014 1011 1008 1036 1060 1018 1102 1015 1022 1010 1018 1018 1016 1028 1035
[176] 1047 1026 1050 1013 1018 1014 1018 1027 1033 1028 1024 1034 1009 1009 1054 1015 1014 1018 1016 1027 1012 1014 1027 1010 1046
[201] 1000 1018 1028 1068 1022 1012 1026 1024 1000 1022 1012 1030 1024 1002 1016 1034 1050 1062 1005 1010 1011 1010 1006 1023 1012
[226] 1018 1018 1032 1042 1015 1054 1048 1010 1016 1018 1042 1049 1004 1042 1047 1014 1028 1008 1043 1008 1010 1026 1014 1028 1037
[251] 1024 1010 1045 1015 1010 1035 1036 1020 1017 1022 1064 1021 1025 1008 1012 1022 1024 1008 1027 1070 1034 1032 1012 1006 1036
[276] 1028 1006 1014 1018 1022 1020 1024 1020 1027 1019 1017 1010 1050 1036 1008 1022 1058 1030 1020 1000 1022 1000 1008 1000 1059
[301] 1013 1004 1019 1014 1000 1018 1020

如何通过仅选择这些观察值中的 60 个来获得正态分布？ 使用以下代码是否正确？

seq = seq(1, 60, 1)
y = rnorm(seq, mean(vars), sd(vars))
y
hist(y)

谢谢

Answer 1

为了扩展我的评论，我加载了该数据并构建了密度图：

vars <- read.table(text="
  [1] 1053 1018 1048 1040 1046 1038 1014 1020 1038 1006 1016 1025 1011 1018 1014 1021 1022 1000 1003 1018 1037 1053 1033 1047 1018
 [26] 1022 1014 1048 1048 1026 1028 1056 1015 1030 1008 1035 1007 1026 1021 1008 1053 1052 1009 1015 1020 1013 1023 1015 1018 1038
 [51] 1007 1018 1035 1030 1016 1034 1028 1018 1018 1024 1024 1030 1028 1027 1016 1012 1023 1022 1045 1012 1019 1010 1086 1016 1034
 [76] 1043 1027 1013 1022 1030 1022 1059 1021 1008 1023 1010 1016 1020 1062 1013 1064 1052 1018 1033 1000 1025 1011 1016 1017 1048
[101] 1030 1022 1018 1020 1011 1095 1038 1018 1042 1007 1048 1038 1024 1016 1008 1051 1014 1020 1024 1008 1046 1045 1064 1075 1035
[126] 1010 1010 1010 1012 1012 1013 1016 1039 1030 1003 1013 1025 1005 1000 1029 1028 1008 1014 1026 1020 1062 1023 1008 1032 1024
[151] 1026 1011 1018 1028 1031 1042 1057 1054 1010 1012 1014 1011 1008 1036 1060 1018 1102 1015 1022 1010 1018 1018 1016 1028 1035
[176] 1047 1026 1050 1013 1018 1014 1018 1027 1033 1028 1024 1034 1009 1009 1054 1015 1014 1018 1016 1027 1012 1014 1027 1010 1046
[201] 1000 1018 1028 1068 1022 1012 1026 1024 1000 1022 1012 1030 1024 1002 1016 1034 1050 1062 1005 1010 1011 1010 1006 1023 1012
[226] 1018 1018 1032 1042 1015 1054 1048 1010 1016 1018 1042 1049 1004 1042 1047 1014 1028 1008 1043 1008 1010 1026 1014 1028 1037
[251] 1024 1010 1045 1015 1010 1035 1036 1020 1017 1022 1064 1021 1025 1008 1012 1022 1024 1008 1027 1070 1034 1032 1012 1006 1036
[276] 1028 1006 1014 1018 1022 1020 1024 1020 1027 1019 1017 1010 1050 1036 1008 1022 1058 1030 1020 1000 1022 1000 1008 1000 1059
[301] 1013 1004 1019 1014 1000 1018 1020", fill=TRUE)
vars <- unlist(vars[-1]) ; remove  extraneous item numbering

现在绘制非 NA 值

png();  plot(density(na.omit(vars)))
 dev.off()

所以基础分布是右偏的，随机抽奖不会是“正态”