M*********g
发帖数: 178 | 1 版权所有：castor_v_pollux 原作 提交时间：16:19:21 08月12日
第四章 白云深处
一
应该说，玻尔关于原子结构的新理论出台后，是并不怎么受到物理学家们的欢迎的。这个
理论，在某些人的眼中，居然怀有推翻麦克斯韦体系的狂妄意图，本身就是大逆不道的。
瑞利爵士（我们前面提到过的瑞利-金斯线的发现者之一）对此表现得完全不感兴趣，J.J
.汤姆逊，玻尔在剑桥的导师，拒绝对此发表评论。另一些不那么德高望重的人就直白多
了，比如一位物理学家在课堂上宣布：“如果这些要用量子力学才能解释的话，那么我情
愿不予解释。”另一些人则声称，要是量子模型居然是真实的话，他们从此退出物理学界
。即使是思想开放的人，比如爱因斯坦和波恩，最初也觉得完全接受这一理论太勉强了一
些。
但是量子的力量超乎任何人的想象。胜利来得如此之快之迅猛，令玻尔本人都几乎茫然而
不知所措。首先，玻尔的推导完全符合巴耳末公式所描述的氢原子谱线，而从W2-W1 = h
ν这个公式，我们可以倒过来推算ν的表述，从而和巴耳末的原始公式ν＝R（1/2^2 -
1/n^2）对比，计算出里德伯常数R的理论值来。而事实上，玻尔理 |
|
f*******n
发帖数: 70 | 2 Right. More precisely, discretize the SDE as:
S(t)-S(t-1)=r*S(t-1)*dt+v*S(t-1)*sqrt(dt)*Normal(0,1),
simulate the N(0,1) r.v. and proceed step by step.
Actually what I usually do is to simulate log(St) first by:
d(lnS)=(r-v^2/2)dt+vdW, and then take the exp, a bit easier to implement.
If W1, W2 are not independent, you need to simulate two correlated N(0,1) r.v.
Oops, it seems to me the "v" in the first equation is the same as the "vt" in
the 2nd(volatility?), then you need to simulate vt first, |
|
l*****i
发帖数: 43 | 3 how to change the output of proc freq
proc freq data=a;
table a *b;
weight w1;
run;
the total freq is like 1.21E7 (1.21*10 7th). The actual number is 12151244.
How SAS put the actual number?
Thanks, |
|
|
d*******1
发帖数: 293 | 5 这个正是我要找的。
如果用SAS的话， 这个sharpiro test 的W 是不是先要存在dataset里面，然后读出来
。我想2000组测试完， 产生这样的输出：
# W_value
1 W1
2 W2
3 W3
......
....
这个不知道好不好弄 |
|
T*******I
发帖数: 5138 | 6 看来，在那些抽象的概念上纠缠永远不能有结果。让我说得具体一点。
给定一个两分法的样本（假定X是自变量而Y是因变量，两段都是简单线性模型，且临界点是在X上）。现行算法及分段模型组的基本表述如下：
hat\y_1 = a1+b1X if X<=t
hat\y_2 = a2+b2X if X>t
ID X Y M CR
1 x1 y1 m1 cr1
2 x2 y2 m2 cr2
3 x3 y3 m3 cr3
4 x4 y4 m4 cr4
5 x5 y5 m5 cr5
6 x6 y6 m6 cr6
7 x7 y7 m7* cr7 min(.)
8 x8 y8 m8 cr8
9 x9 y9 m9 cr9
0 x0 y0 m0 cr0
其中，M是由分段模型组的系数构成的矩阵，CR是分段模型的合并残差。*表示根据最小
CR选定的分段模型，如果我们有 ... 阅读全帖 |
|
T*******I
发帖数: 5138 | 7 这几天版上总有人找我茬。相信他们无一不是数学背景出来搞统计的。他们以为自己掌
握了一点数学技能就在统计学里自命不凡。如果他们不是孬种，就请接受我的以下挑战
，并回答我在最后提出的简单问题。回答不了的，或不敢回答的，就请他/她滚回数学里
去讨饭吃，别仗着自己那份高深莫测的数学理论继续在统计学里胡说八道。为了不再继
续为版上添乱，我想请seattleren, ningyan, kaleege等人接受我的挑战。当然，我也欢
迎任何人参与严肃的讨论。不能说出个一二三四的，就请自动回避，免得自讨没趣（我
想对pp65说的是，我对你感到抱歉，因为本段最后的话对你来说说得太晚了）。
给定一个两分法的样本（假定X是自变量而Y是因变量，两段都是简单线性模型，且临界
点是在X上）。现行算法及分段模型组的基本表述如下：
hat\y_1 = a1+b1X if X<=t
hat\y_2 = a2+b2X if X>t
ID X Y M CR
1 x1 y1 m1 cr1
2 x2 y2 m2 cr2
3 x3 y3 ... 阅读全帖 |
|
T*******I
发帖数: 5138 | 8 你是真看不懂还是故意装14？要不我把问题重复在这里：
给定一个两分法的样本（假定X是自变量而Y是因变量，两段都是简单线性模型，且临界
点是在X上）。现行算法及分段模型组的基本表述如下：
hat\y_1 = a1+b1X if X<=t
hat\y_2 = a2+b2X if X>t
ID X Y M CR
1 x1 y1 m1 cr1
2 x2 y2 m2 cr2
3 x3 y3 m3 cr3
4 x4 y4 m4 cr4
5 x5 y5 m5 cr5
6 x6 y6 m6 cr6
7 x7 y7 m7* cr7 min(.)
8 x8 y8 m8 cr8
9 x9 y9 m9 cr9
0 x0 y0 m0 cr0
我对上述方法进行了如下改造：
hat\y_1 = a1+b1X if X<=t_bar (t: Threshold)
... 阅读全帖 |
|
c*********r
发帖数: 1802 | 9 处理数据中，按其中一个变量的数值大小，分成10组，标记成w1-w10。
共168个obs，没有缺失值，用proc sort将该变量排了序，
然后用if, else if, else分组，（检查了N遍，确信中间没有断档）
结果只得到8组数据，最后2个组的数据丢了，run了很多遍，不知道为何？
抓狂中。。。
紧急求助，包子答谢！ |
|
c*********r
发帖数: 1802 | 10 没有人遇到过类似的问题吗？
data category;
set cate;
if TFW< 5 then Tlevel="W1";
else if 5<=TFW<8 then Tlevel="W2";
else if 8<=TFW<11 then Tlevel="W3";
else if 11<=TFW<14 then Tlevel="W4";
else if 14<=TFW<17 then Tlevel="W5";
else if 17<=TFW<20 then Tlevel="W6";
else if 20<=TFW<23 then Tlevel="W7";
else if 23<=TFW<26 then Tlevel="W8";
else if 26<=TFW<29 then Tlevel="W9";
else if 29<=TFW<32 then Tlevel="W10";
else if 32<=TFW<50 then Tlevel="W11";
else Tlevel="W12";
run;
proc print data=category;
by Tlevel;
ID ... 阅读全帖 |
|
G*******s
发帖数: 10605 | 11 呵呵，这么搞W10,W11的都去W1了，定义字符串变量第一次赋值会确定长度，这种情况
你要先定义长度，length Tlevel $3 就可以了 |
|
b****w
发帖数: 71 | 12 多谢大家回帖，我针对我之前没有说清楚的地方进行了一点补充。
面试官很Nice，一步步引导我啊，可惜我还是不知道如何解。。又悲剧了一个phone
interview啊。下面是两道新鲜出炉的面试题，请教大家
1.针对同一个data set,现有两个model A,B,可能是用不同的方法做出来的，DV是一样的，IV没说。 已知outcome(A)=100,outcome(B)=80, var(A)=25,var(B)=49,并且
两个Model的error都是iid的，先要把两个model结合起来,赋不同的weight, minimize
variance.怎么弄？我当时说是不是就是求Min var（w1*A+w2*B），他说不是。我真不
知道是如何结合。而且我觉得应该是要用到error iid这个条件
2.这个简单点，但当时我已经秀逗了，完全转不了脑子了，嗨。。只能怪经验太少以求
安慰了
具体我就不说了，反正我得到最后的一步就是binominal, (20,n)*(0.15)^n*(0.85)^(
20-n)=0.12,求n.那个（20，n）就是20 choose n，我不知道怎么答。这个... 阅读全帖 |
|
p******s
发帖数: 229 | 13 please see the log as follows:
40 proc import datafile= 'C:Projectsrawdata.xls' out= w1
41 dbms= excel ;
The SAS/ACCESS Interface to PC Files is not installed. Please install this
module in order to IMPORT/EXPORT to these file types.
NOTE: The SAS System stopped processing this step because of errors.
Then I checked if SAS/ACCESS is installed. please see the following:
42 proc setinit ; run ;
NOTE: PROCEDURE SETINIT used (Total process time):
System birthday: 27FEB2013.
Operating Sys... 阅读全帖 |
|
c**********e
发帖数: 2007 | 14 Suppose one of the inputs is a macro variable. The macro variable is a list
of variable name separated by space.
For example, in the macro, the list of predictors is the input:
%macro regression(ind_var);
proc reg data=mydata;
model y = &ind_var.;
run;
%mend;
A call of the macro takes the following form:
%regression(x1 x2 x3 w1 w2);
My question is how to get the number of predictors. In the above macro call,
I want to get the number 5, which is the number of independent variables.
How to do... 阅读全帖 |
|
a******1
发帖数: 201 | 15 Then how good is your statistics teacher?
In the simple example I gave, to calculate the r as OP did, you need to
calculate variance first. What is the variance in my example? The first
sample is (h1, h2), and the second is (w1, w2), they don't even measure the
same thing, you can't get any meaningful variance from these two "samples". |
|
a******1
发帖数: 201 | 16 Actually, what is the mean, (h1+w1)/2 and h2+w2)/2?
the |
|
|
k**********g
发帖数: 33 | 18 1、美国德州太平洋集团出资50亿元人民币设立的德太中国西部成长基金，并将 中
国西部总部设在重庆两 江新区；
2、两江新区开发投资集团与上海张江（集团）有限公司将在水土高新技术产业 园
内携手打造2平方公里的 “重庆张江信息科技产业园”，投资额为200亿元人民币；
3、西南地区灾备指挥、后勤物资保障中心落户两江；
4、产品广泛用于美国潜艇与航母的挪威RPR科技公司在渝设立亚太地区总部；( ^: @%
e7 f/ A: U9 b
5、英国诺实科技公司投资10亿英镑，建设“两江低碳产业园”；
6、北大方正医药总部基地落户两江，投资50亿，成为西南地区规模最大的生物 医药
总部基地；
7、工商银行、农业银行、中国银行、国家开发银行、华夏银行、重庆农村商业 银
行、西南证券、重庆银 行等15家重点金融机构，已经在江北嘴CBD落户设立区域总部
；* {5 p2 m9 b" W, Z* @
8、总投资20亿元、年物资周转量超过500亿元的中国移动西南大区物流基地在两 江
新区空港功能区开建；, @) W4 j$ A* p# Y9 ?
9、国家海上风电工程技术研究... 阅读全帖 |
|
z*******n
发帖数: 1034 | 19 现如今， 三大作业系统内核的调度单元皆基于线程，多线程场景更是近乎100%，
静态后端的开发需要log一些信息到file上，以助于多线程调试，有log信息查看总是幸
福。
这么多日子， c/c++始终没有提供类似NSLog的thread safe function call在标准库中，
无奈之下，需要轮子，
静态后端评估轮子的第一要求就是这个公司或组织有人在标准委员会中，
Bjarne Stroustrup抱怨boost模块化差，不能同意更多，排除之外，找到google的glog，
接续下来，评估glog |
|
W*********n
发帖数: 1694 | 20 【 以下文字转载自 WaterWorld 讨论区 】
发信人: laminin (冲击宇宙), 信区: WaterWorld
标 题: 色字头上一把刀，德国ONS钓鱼门
关键字: 一夜情,开元论坛,钓鱼门,ONS
发信站: BBS 未名空间站 (Mon Feb 15 22:43:16 2010, 美东)
德国的开元论坛，上面有个板块叫“一夜情缘”，是给欧洲的XDJM找乐子的。现在有这
么一个上海MM，自称“小眼睛”，发了一个帖子找GG来happy，众多欧洲WSN们蜂拥而至
，纷纷发了自己的靓照和情话给这位MM。这位MM也够狠，充分发挥了钓鱼精神，把众多
GG的来信和pp全部发在了自己的MSN space上面供人瞻仰。色字头上一把刀啊，北美的
众位WSN们要引以为戒啊。
以下为这位MM的原帖
〔一夜情缘] 大家好，我是慕尼黑的小眼睛
( Z本帖最后由 aiyawawa 于 11.12.2009 02:16 编辑
刚刚来慕尼黑，听朋友说有个开元论坛，上来一看吓了我一条，国外真的好开放呀，可
以光明正大的2 W! B# Z7 w1 V Q9 P4 V1 }
其实还没有做好心理准备，先看看 |
|
|
|
|
|
w*******y
发帖数: 60932 | 25 Wired Magazine just rolled out a new subscription model where you "subscribe
once and access anywhere". Previously, you were limited to paying for
individual issues (e.g. iPad app), but now your physical subscription
carries over. This is now LIVE (already done it myself).
Before the trolls come, each issue costs $4 a pop if on the iPad w/no
subscription. This amounts to about $24 savings a year at least over buying
individual issues (and in my case about $44 since I got in on one of those
cheap... 阅读全帖 |
|
|
|
|
|
|
|
n******e
发帖数: 623 | 32 直接把www.aerith.net的Weekly Information抄过来了：
It is brightening much faster than expected. It is already so bright as 6.2
mag (May 14, Marco Goiato). The apparent size is so large as 20 arcmin. It
will pass 0.85 A.U. from the sun in late June, and it may reach up to 4.5
mag. In the Northern Hemipshere, it is observable until around May 25, when
the comet will be 6 mag. It will be unobservable for one and a half month
around the perihelion passage. But it will appear in the morning sky again
at 5.5 m |
|
