一个统计问题：怎么评判哪种测量方法好 - 未名空间MITBBS历史存档

国际科技财经博客移民网络热点娱乐民生时事公众号

Redian新闻

>未名空间

>Biology - 生物学

一个统计问题：怎么评判哪种测量方法好

一个统计问题：怎么评判哪种测量方法好# Biology - 生物学

c*02015-08-04 07:08

1 楼

从小特羡慕会弹钢琴的父母，咱两口子都是那种五音不全的主。儿子快四岁，特调皮，
没一刻坐的住。
小子从娘胎就喜欢音乐，并且越是摇滚越喜欢。不管在家还是在daycare，从小就喜欢
随着音乐摇摆。
去朋友家发现大部分家庭都摆着一架钢琴，主人到底弹的怎么样，一次也没欣赏过。唯
一对弹钢琴的印象是10多年前在国内时，带当初才十三四岁的一大姐的儿子去逛商场。
在钢琴部摆着一架上百万的钢琴，所有客人都在那望洋兴叹，每一个敢试的，有几个和
我一样五音不全的也只是把那几万十几万的钢琴噼里啪啦按两下。
大姐的儿子长的很帅，从小上贵族学校，那小子往那琴凳上一坐，马上一曲施特劳斯还
是谁的名曲就从他指下流畅出来。
打那时起，就幻想着将来也要自己的孩子当众这样炫一下。
但我家这小子横看竖看顶多也就一个搞运动的料，长的黑黑壮壮的。怎么要他去练钢琴
呀。加上我俩，别说钢琴，就是口琴，当初花了好长时间也没吹出一个1234567来。
现在特向狠点心，买架好点的钢琴摆家里，强迫这小子坐旁边。btw，以前也是因为这
个幻想，给这小子买过一架电子琴，他那时太小，也就偶尔碰一下。大多是我自己玩。
我不会谈，只会用它放音乐。
我有个音乐家朋友，在维也纳金色大厅演出过，然后在纽约的林肯中心和卡耐基表演过
几年。搞民乐的，在这边不吃香。现在老公有点事业，自己也不用再去忙活了。平时没
事就喝喝咖啡，好像钢琴谈的不错，但她本人特讨厌音乐，在家就喜欢静静的躺着。准
备等儿子过了四岁后要他帮忙随便调教一下，改良一下我家祖传的基因。人家家里基因
好，音乐世家。
不知大伙认为我值不值得狠点心，买架钢琴当家具摆家里熏陶一下这小子。
万一有客人看见我家有钢琴，要我们夫妇谁来一曲，那是不是囧的很呀！

c*e2015-08-04 07:08

2 楼

我2个月前因为没有找到工作，file了F1 to F2 的application, 现在还在processing
过程中。但是就在昨天，收到一个
offer，一个公司愿意sponsor我H1B，我能不能现在申请h1b? 还是需要等f2批下来之后
才能申请h1b? 这个情况需要请律
师么？

r*e2015-08-04 07:08

3 楼

老虎的屁股得这么摸 .gif

A*H2015-08-04 07:08

4 楼

有一百个样品。每一个样品的实际值是1，2，3，或者4。观察值范围变动大概是0.5到
12左右，也就是说，对于实际值是2的样品，有可能观察到1.8或者5.1。我想要的是一
个方法能有效的把样品分到四组，对应于实际值1，2，3，4
有几种测量方法，每一种都会得到这么一个样子的数据
实际值观测值
2 3.4133
2 3.3333
2 3.0588
3 6.9932
4 7.0595
3 6.5312
2 6.4757
1 1.2419
4 10.337
2 6.0229
2 7.2705
4 8.1558
1 0.7025
怎么评判哪种测量方法好
大概可以看到观测值随着实际值增大而增大，correlation也支持这一点
补充：
我可以使用几种不同的测量方法（或者说计算方法），比如说甲乙丙三种方法，我想知
道哪种方法能让我更好的判断某个样品的实际值属于四组里面的哪一个小组（也就是实
际值是1，2，3，或者4）
给出的数据是一部分，真正的变化有点大，比如说
方法甲会产生
0.5～3.1观测值for实际值1
1.1～6.2观测值for实际值2
1.0～8.1观测值for实际值3
2.5～12.1观测值for实际值4
再：
不知道multinomial goodness of fit分析是不是合适
看了下，似乎不对。。。
再2：
举一个极端的例子
如果方法丙能产生
0.8～1.2观测值for实际值1
1.8～2.2观测值for实际值2
2.8～3.2观测值for实际值3
3.8～4.2观测值for实际值4
那这个方法肯定比上面说的方法甲要好
贴两个方法的data summary:
Method X:
> summary(r1$V2) #--- group for expected value = 1
Min. 1st Qu. Median Mean 3rd Qu. Max.
1.019 1.694 2.095 2.110 2.451 3.435
> summary(r2$V2) #--- group for expected value = 2
Min. 1st Qu. Median Mean 3rd Qu. Max.
1.892 3.235 3.865 4.215 4.917 9.557
> summary(r3$V2) #--- group for expected value = 3
Min. 1st Qu. Median Mean 3rd Qu. Max.
5.122 6.438 7.031 7.491 8.155 10.760
> summary(r4$V2) #--- group for expected value = 4
Min. 1st Qu. Median Mean 3rd Qu. Max.
6.732 7.269 8.644 8.975 10.110 12.370
Method Y:
> summary(r1$V2)
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.7025 1.2420 1.2720 1.6030 1.8330 3.1200
> summary(r2$V2)
Min. 1st Qu. Median Mean 3rd Qu. Max.
1.210 2.344 3.181 3.633 4.895 10.270
> summary(r3$V2)
Min. 1st Qu. Median Mean 3rd Qu. Max.
1.574 3.784 4.801 5.182 6.380 11.590
> summary(r4$V2)
Min. 1st Qu. Median Mean 3rd Qu. Max.
5.232 6.876 7.547 8.473 8.701 15.580

i*e2015-08-04 07:08

5 楼

不缺银两的话，我支持你买。反正你的目的是“买架钢琴当家具摆家里熏陶一下这小子
”，这是最保值的家具了。
客人来了也没关系的，你说自己其实不会弹，就是为了给懂音乐的客人准备的。你自己
肯定是不囧的，客人是不是囧就不知道了。 :)

O*d2015-08-04 07:08

6 楼

Youtube is here
http://www.youtube.com/watch?v=1AZn5nWIj_g

D*a2015-08-04 07:08

7 楼

啥叫“哪种测量方法好“？你是要测量方法好还是要分四组？
观察值范围变动大概是0.5到
12左右，也就是说，对于实际值是2的样品，有可能观察到1.8或者5.1
---这个变动范围不是3么，怎么是12？到底是对于实际值是2的样品有可能观察到14还
是？

k*n2015-08-04 07:08

8 楼

看得我好囧:))
"囧"这个字怎么念阿?拼音打不出来

【在 i******e 的大作中提到】

: 不缺银两的话，我支持你买。反正你的目的是“买架钢琴当家具摆家里熏陶一下这小子
: ”，这是最保值的家具了。
: 客人来了也没关系的，你说自己其实不会弹，就是为了给懂音乐的客人准备的。你自己
: 肯定是不囧的，客人是不是囧就不知道了。 :)

f*22015-08-04 07:08

9 楼

这个youtube的top comment太逗了

【在 O*******d 的大作中提到】

: Youtube is here
: http://www.youtube.com/watch?v=1AZn5nWIj_g

A*H2015-08-04 07:08

10 楼

我可以使用几种不同的测量方法（或者说计算方法），比如说甲乙丙三种方法，我想知
道哪种方法能让我更好的判断某个样品的实际值属于四组里面的哪一个小组（也就是实
际值是1，2，3，或者4）
给出的数据是一部分，真正的变化有点大，比如说
方法甲会产生
0.5～3.1观测值for实际值1
1.1～6.2观测值for实际值2
1.0～8.1观测值for实际值3
2.5～12.1观测值for实际值4
举一个极端的例子
如果方法丙能产生
0.8～1.2观测值for实际值1
1.8～2.2观测值for实际值2
2.8～3.2观测值for实际值3
3.8～4.2观测值for实际值4
那这个方法肯定比上面说的方法甲要好

【在 D*a 的大作中提到】

: 啥叫“哪种测量方法好“？你是要测量方法好还是要分四组？
: 观察值范围变动大概是0.5到
: 12左右，也就是说，对于实际值是2的样品，有可能观察到1.8或者5.1
: ---这个变动范围不是3么，怎么是12？到底是对于实际值是2的样品有可能观察到14还
: 是？

i*e2015-08-04 07:08

11 楼

我也不会念，但意思明白。我copy+paste的。
我是很认真的！钢琴真的可以当家具，好的真的保值。

【在 k**n 的大作中提到】

: 看得我好囧:))
: "囧"这个字怎么念阿?拼音打不出来

f*22015-08-04 07:08

12 楼

这个youtube的top comment太逗了

【在 O*******d 的大作中提到】

: Youtube is here
: http://www.youtube.com/watch?v=1AZn5nWIj_g

D*a2015-08-04 07:08

13 楼

那这种变动符合正态分布么？比如把所有数据都放一起能否看到几个明显的峰？

【在 A**H 的大作中提到】

: 我可以使用几种不同的测量方法（或者说计算方法），比如说甲乙丙三种方法，我想知
: 道哪种方法能让我更好的判断某个样品的实际值属于四组里面的哪一个小组（也就是实
: 际值是1，2，3，或者4）
: 给出的数据是一部分，真正的变化有点大，比如说
: 方法甲会产生
: 0.5～3.1观测值for实际值1
: 1.1～6.2观测值for实际值2
: 1.0～8.1观测值for实际值3
: 2.5～12.1观测值for实际值4
: 举一个极端的例子

k*n2015-08-04 07:08

14 楼

我也是这么觉得的,其实楼主是为自己买琴,打着孩子的旗号.
或者,如果不买那个几百万的,买个几千块的摆在家里真的没什么,就是自己消费罢了.买
个花瓶和买个琴,不就是花钱多少的问题吗.以后不要逼孩子练就是了,当个玩具就成.

【在 i******e 的大作中提到】

: 我也不会念，但意思明白。我copy+paste的。
: 我是很认真的！钢琴真的可以当家具，好的真的保值。

z*i2015-08-04 07:08

15 楼

zan

【在 r*********e 的大作中提到】

: 老虎的屁股得这么摸 .gif

A*H2015-08-04 07:08

16 楼

应该不太能看到，因为样品里面的有实际值1，2，3，4的数目是不一样的，比如说30个
1，20个2，5个3，15个4

【在 D*a 的大作中提到】

: 那这种变动符合正态分布么？比如把所有数据都放一起能否看到几个明显的峰？

i*e2015-08-04 07:08

17 楼

是真买琴么？我怎么觉得在挖坑。

【在 k**n 的大作中提到】

: 我也是这么觉得的,其实楼主是为自己买琴,打着孩子的旗号.
: 或者,如果不买那个几百万的,买个几千块的摆在家里真的没什么,就是自己消费罢了.买
: 个花瓶和买个琴,不就是花钱多少的问题吗.以后不要逼孩子练就是了,当个玩具就成.

A*H2015-08-04 07:08

18 楼

Ultimately I need develop a method that can assign samples based on observed
values into groups (with expected values of 1, 2, 3, or 4). Of course the
more correctly assigned samples, the better.

c*02015-08-04 07:08

19 楼

你太敏感了。
是真的想买件家具，现在在看yamaha u1, u3。
贵的买不起，买了对我们这种土包子也是暴殄天物。
太差的话，怕人家钢琴老师弹不下去。

【在 i******e 的大作中提到】

: 是真买琴么？我怎么觉得在挖坑。

u*h2015-08-04 07:08

20 楼

讨论一下：
假设：
对于 1 观察值分布是否是正态分布？
对于 2, 3, 4 是不是也这样？
如果是这样，是不是可以用multinomial logistic?

【在 A**H 的大作中提到】

: 有一百个样品。每一个样品的实际值是1，2，3，或者4。观察值范围变动大概是0.5到
: 12左右，也就是说，对于实际值是2的样品，有可能观察到1.8或者5.1。我想要的是一
: 个方法能有效的把样品分到四组，对应于实际值1，2，3，4
: 有几种测量方法，每一种都会得到这么一个样子的数据
: 实际值观测值
: 2 3.4133
: 2 3.3333
: 2 3.0588
: 3 6.9932
: 4 7.0595

m*a2015-08-04 07:08

21 楼

显然不用给儿子买啊
当家具就给爸妈买了好了
估计儿子不会学滴

A*H2015-08-04 07:08

22 楼

应该不是很正态分布，不过我想也可以大约到正态分布吧
我查查看multinomial logistic去

【在 u*h 的大作中提到】

: 讨论一下：
: 假设：
: 对于 1 观察值分布是否是正态分布？
: 对于 2, 3, 4 是不是也这样？
: 如果是这样，是不是可以用multinomial logistic?

t*r2015-08-04 07:08

23 楼

在钢琴坑前果断 pass 。。。

A*H2015-08-04 07:08

24 楼

贴两个方法的data summary:
Method X:
> summary(r1$V2) #--- group for expected value = 1
Min. 1st Qu. Median Mean 3rd Qu. Max.
1.019 1.694 2.095 2.110 2.451 3.435
> summary(r2$V2) #--- group for expected value = 2
Min. 1st Qu. Median Mean 3rd Qu. Max.
1.892 3.235 3.865 4.215 4.917 9.557
> summary(r3$V2) #--- group for expected value = 3
Min. 1st Qu. Median Mean 3rd Qu. Max.
5.122 6.438 7.031 7.491 8.155 10.760
> summary(r4$V2) #--- group for expected value = 4
Min. 1st Qu. Median Mean 3rd Qu. Max.
6.732 7.269 8.644 8.975 10.110 12.370
Method Y:
> summary(r1$V2)
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.7025 1.2420 1.2720 1.6030 1.8330 3.1200
> summary(r2$V2)
Min. 1st Qu. Median Mean 3rd Qu. Max.
1.210 2.344 3.181 3.633 4.895 10.270
> summary(r3$V2)
Min. 1st Qu. Median Mean 3rd Qu. Max.
1.574 3.784 4.801 5.182 6.380 11.590
> summary(r4$V2)
Min. 1st Qu. Median Mean 3rd Qu. Max.
5.232 6.876 7.547 8.473 8.701 15.580

i*e2015-08-04 07:08

25 楼

既然真不是挖坑。我就认真答一下你原帖最后一句的concern.
不管怎么说，学琴的孩子多过大人。我们家那“家具”摆了至少十年了，（都记不清那
年买的了）。从来没有客人要求我们夫妇弹一曲，大家都默认是我家孩子在学琴。同样
，我到别人家，除非知道主人是音乐人，或没小孩，我也默认是那家孩子在弹琴。
就算是真有人叫你弹一曲，就直说你不会弹，是儿子要学就好了，或是儿子也不学了，
现在是100%家具了。我觉得这一点也不难堪。能买得其6千刀的家具摆着，说明你有一
定的经济势力，这挺光荣的事儿。

【在 c*********0 的大作中提到】

: 你太敏感了。
: 是真的想买件家具，现在在看yamaha u1, u3。
: 贵的买不起，买了对我们这种土包子也是暴殄天物。
: 太差的话，怕人家钢琴老师弹不下去。

y*s2015-08-04 07:08

26 楼

ANCOVA, 测量方法 being the independent variable, 观测值 being the outcome,
实际值being the covariate.

c*02015-08-04 07:08

27 楼

大家说的对，我打心里是为了满足自己的心理需求。
早些年人穷，什么都没学上。什么乐器都不懂，但我还是特喜欢听音乐。只是命比较苦
，这辈子只有当观众或者听众的命了。
我不想儿子将来也和我一样，凡事都是观众加听众。碰到什么他喜欢的玩具，都会帮他
买下来。自己儿时未能圆的梦，没能尝的心愿。也希望儿子能代替自己什么都走一趟。

【在 i******e 的大作中提到】

: 既然真不是挖坑。我就认真答一下你原帖最后一句的concern.
: 不管怎么说，学琴的孩子多过大人。我们家那“家具”摆了至少十年了，（都记不清那
: 年买的了）。从来没有客人要求我们夫妇弹一曲，大家都默认是我家孩子在学琴。同样
: ，我到别人家，除非知道主人是音乐人，或没小孩，我也默认是那家孩子在弹琴。
: 就算是真有人叫你弹一曲，就直说你不会弹，是儿子要学就好了，或是儿子也不学了，
: 现在是100%家具了。我觉得这一点也不难堪。能买得其6千刀的家具摆着，说明你有一
: 定的经济势力，这挺光荣的事儿。

s*e2015-08-04 07:08

28 楼

感觉可以用每个方法得出的数据做cluster analysis，结果和真正值用chi square去比
，看哪个跟真正值最吻合。

t*r2015-08-04 07:08

29 楼

楼上的，考虑买 Yamaha Disklavier，让钢琴自己弹。
其实我也有这个想法，不过我家已经有了 Neato Robotic Vacuum，暂时
对买 Disklavier 很犹豫。主要是担心家里没人的时候，俩机器人自己打
起来。。。楼主有没有啥好点的建议。。。

s*j2015-08-04 07:08

30 楼

就LZ给的数据, 你这样算出来的”最吻合”结果.
你敢信吗?

【在 s**e 的大作中提到】

: 感觉可以用每个方法得出的数据做cluster analysis，结果和真正值用chi square去比
: ，看哪个跟真正值最吻合。

c*02015-08-04 07:08

31 楼

别吓唬我了，我是经常喜欢到商店看上百万的钢琴，但说搬到家，咱可是想都不敢想，
这么大个，看着都吓着我了。这种东西，到酒店，商场摸摸就行了。

【在 t*******r 的大作中提到】

: 楼上的，考虑买 Yamaha Disklavier，让钢琴自己弹。
: 其实我也有这个想法，不过我家已经有了 Neato Robotic Vacuum，暂时
: 对买 Disklavier 很犹豫。主要是担心家里没人的时候，俩机器人自己打
: 起来。。。楼主有没有啥好点的建议。。。

A*H2015-08-04 07:08

32 楼

数据和和ANCOVA的例子倒是长得有点像
不过ANCOVA是分析variance或者test hypothesis，不知道怎么用来评判我的方法中哪
个更好
我再学学去。。。

【在 y***s 的大作中提到】

: ANCOVA, 测量方法 being the independent variable, 观测值 being the outcome,
: 实际值being the covariate.

n*h2015-08-04 07:08

33 楼

想买钢琴当家具，　有便宜的钢琴．　熏陶孩子，　买钢琴是做不到的．　很多家里买
了琴，　孩子也不学，　就当个家具．　真要熏陶，　你给孩子放ＣＤ就行了．

【在 c*********0 的大作中提到】

u*h2015-08-04 07:08

34 楼

一个建议：
choose values for a1 < a2 < a3 as separator values,
for example, a1 = 1.3, a2 = 2.8, a3=8.6;
For all the sample, based on the observation value:
if obv < a1, Y_hat = 1
if a1 < obv < a2, Y_hat = 2
if a2 < obv < a3, Y_hat = 3
if a3 < obv, Y_hat = 4
count the sum for Y_real <> Y_hat, call it error_sum
for each method, adjust a1, a2, a3 so the error_sum is minimum,
To compare between methods, the smaller error_sum, the better method.
note: each method should have it own a1, a2, a3 values.

【在 A**H 的大作中提到】

t*r2015-08-04 07:08

35 楼

其实 Disklavier 根本不贵，就比对应的三角琴贵上两三千左右好像。
如果家里地方小点，没关系，就这点小事儿能难道 Yamaha 的码工么？
Yamaha 也有 Disklavier 的立式钢琴版，样子跟立式钢琴一模一样。
价钱也就比对应的立式钢琴稍稍贵一点点：Yamaha DU1E3
http://usa.yamaha.com/products/musical-instruments/keyboards/di

【在 c*********0 的大作中提到】

: 别吓唬我了，我是经常喜欢到商店看上百万的钢琴，但说搬到家，咱可是想都不敢想，
: 这么大个，看着都吓着我了。这种东西，到酒店，商场摸摸就行了。

s*s2015-08-04 07:08

36 楼

exactly. lz这数据再怎么model都没用，就三个separator，最小化错误率就行了。

【在 u*h 的大作中提到】

: 一个建议：
: choose values for a1 < a2 < a3 as separator values,
: for example, a1 = 1.3, a2 = 2.8, a3=8.6;
: For all the sample, based on the observation value:
: if obv < a1, Y_hat = 1
: if a1 < obv < a2, Y_hat = 2
: if a2 < obv < a3, Y_hat = 3
: if a3 < obv, Y_hat = 4
: count the sum for Y_real <> Y_hat, call it error_sum
: for each method, adjust a1, a2, a3 so the error_sum is minimum,

c*02015-08-04 07:08

37 楼

这玩意不便宜，买回家不仅仅是我们压力大，儿子的鸭梨那肯定是山大了。
我们被客人笑话不要紧，老子不能的肯定要转嫁到儿子身上。
要不客人还不偷偷说某某一家都那基因，家里还摆个4万多的大钢琴。

【在 t*******r 的大作中提到】

: 其实 Disklavier 根本不贵，就比对应的三角琴贵上两三千左右好像。
: 如果家里地方小点，没关系，就这点小事儿能难道 Yamaha 的码工么？
: Yamaha 也有 Disklavier 的立式钢琴版，样子跟立式钢琴一模一样。
: 价钱也就比对应的立式钢琴稍稍贵一点点：Yamaha DU1E3
: http://usa.yamaha.com/products/musical-instruments/keyboards/di

c*02015-08-04 07:08

38 楼

DGB1KE3
这个还不错，三角琴，估计两万以内都能拿下，也是你说的那个机器人。
要是儿子真能试出点天分，这个投资肯定值。

t*r2015-08-04 07:08

39 楼

要是有钱想买家具，这个挺好，吃饭时也可以顺便放段啥屎特拉闻司机
的啥曲子。。。
其实您老要是有兴趣自己写/改编曲子，自己搞一段，那就真显牛鼻了。
客人来了 show off 自己的最新作品，搞不好是真高尚。。。

【在 c*********0 的大作中提到】

: DGB1KE3
: 这个还不错，三角琴，估计两万以内都能拿下，也是你说的那个机器人。
: 要是儿子真能试出点天分，这个投资肯定值。