由买买提看人间百态

topics

全部话题 - 话题: obs
1 2 3 4 5 6 7 8 9 10 下页 末页 (共10页)
w**s
发帖数: 26
1
来自主题: Statistics版 - SAS Regression Excluding OBS i
请教如何run regression by groups excluding obs i in the group to get the
parameter estimates and then get the predicted value for i based on these
estimates. 比如说一组里有20个obs, 那么for the predicted value of obs 1, use
obs 2-20; for obs 2, use obs 1, 3-20. etc.
基本的model是这样的:
proc reg data=test noprint;
model y=x1 x2 x3 /r;
by group;
output out=reg r=residual;
run;
问题是如何run the above regression excluding obs i in the group.
多谢!!!
S******y
发帖数: 1123
2
/* Macro to compare 2 SAS data sets using set operator: except */
%macro cmpar_d2(lib1=,mem1=,lib2=,mem2=);
/* Get the obs in lib1.mem1 but not in lib2.mem2 */
/* and store OBS and VARS in ONEXTWO Dataset */
proc sql;
create table onextwo as
select *
from &lib1..&mem1
except corr all
select *
from &lib2..&mem2
;
quit;
/* Get the obs in lib2.mem2 but not in lib1.mem1 */
/* and store OBS and VARS in TWOXONE Dataset */
proc sql;
create table t
l******1
发帖数: 86
3
谢谢,大家帮忙。是我表达的不够清楚。
其实因为是两个人进行data entry。我用proc compare可以查到有哪些variable在哪个
obs不一样,还可以告诉# of obs in each data set。但是,不能告诉我缺少哪个obs.
s******r
发帖数: 1524
4
data score;
do score=1 to 88;
output;
end;
run;
%macro subgrp(set=,var=, grp=);
proc sort data=&set;by &var;run;
data _null_;set &set end=eof;
if eof then call symput ('total',_n_);
run;
%let step=%eval(&total/&grp);
%if %eval(&step*&grp) < &total %then %let step=%eval(&step+1);
%do i =0 %to %eval(&grp-1);
%let first_obs=%eval(&i*&step+1);
%let obs=%eval(&first_obs+&step -1);
data &set&i;set &set(firstobs=&first_obs obs=&obs);run;
%end;
%mend;
%subgrp(set=score,var=score, grp=5);
bz pls.
C***i
发帖数: 486
5
来自主题: DataSciences版 - 如何用不同来源的Obs来fit 一个 Dist.
一个project 中的一个小环节,有个问题向各位高人请教下。问题抽象出来是,客户找
到一些历史数据 (数据源=20)。因为是历史数据,而且报告地点比较离散,数据质量
很差。具体来讲,数据的样本不同,观测条件也不一致。。。 大概看起来如下:
Study 1: # of obs 100, [1,2,3,...., 100], location A, sample condition I,
etc...
Study 2: # of obs 10, [5, 7, 9,....], location B, sample condition I, etc...
Study 3: # of obs 50, [20, 25, 30, 35, ...], , location C, sample condition
II, etc...
客户想把这些数据都pool 在一起 fit 一个distribution. 我认为因为不同study, 样
本数量不同,不能简单的混在一起。因为一定要给出一个dist., 我的想法是:
1. 根据现有的conditions, merge 类似的数据,比如合并 在同一个... 阅读全帖
p********a
发帖数: 5352
6
来自主题: Statistics版 - [合集] 多obs算出一个值怎么弄呢?
☆─────────────────────────────────────☆
daisytang (daisy) 于 (Mon Mar 10 10:42:07 2008) 提到:
数据大约是这样的
obs month value
1 1 112
1 2 323
1 3 233
1 4 554
1 5 355
2 1 343
2 2 344
2 3 444
2 4 666
2 5 554
3 1 665
3 2 545
3 3 545
3 4 554
3 5 454
4 1 545
4 2 767
4 3 562
4 4 977
4 5 542
算月的value的平均值,最后数据应该是这样的
obs meanofvalue
1 344
2 546
3 767
4 445
怎么用sas做
s*******2
发帖数: 791
7
来自主题: Statistics版 - 怎样得到OBS的值?
用proc contents可以在output里看到OBS;印象中有一个方程,可以把这个output里面
的OBS值直接输出。
谁知道怎么得到这个值?什么方程吗?
A*******s
发帖数: 3942
8
suppose we use a procedure to get a table from output mode, for example:
output out=dataset ...;
or use ODS output table like:
ods output ods_table_name=dataset;
However in neither way can i use obs= dataset options to restrict the number of observations in output dataset, the following code doesnt work
ods output ods_table_name=dataset(obs=1);
But other dataset options seems OK, say where=, keep=, drop=, ...
couldnt understand why
z******e
发帖数: 302
9
单张面值:
70K-130K
可接受的价格(必须明码标价!):
$11 for 1K, OBS
物品新旧要求:
邮寄方式要求:
EMAIL
买卖双方谁承担邮寄损失(Required if not code only):
付款方式说明:
paypal
其他补充说明:
广告的有效期:
UNTIL GONE
物品来源:
各种
我的联系方式:
pm.
email: s*******[email protected]
z******e
发帖数: 302
10
单张面值:
70K-130K
可接受的价格(必须明码标价!):
$11 for 1K, OBS
物品新旧要求:
邮寄方式要求:
EMAIL
买卖双方谁承担邮寄损失(Required if not code only):
付款方式说明:
paypal
其他补充说明:
广告的有效期:
UNTIL GONE
物品来源:
各种
我的联系方式:
pm.
email: s*******[email protected]
p********a
发帖数: 5352
11
来自主题: Database版 - How to put OBS number into the data?
I need to create data with the obs # as a new variable. Is there any way to
do that?
Thanks
D******n
发帖数: 2836
12
now i am confused....sigh...

obs.
l******1
发帖数: 86
13
呵呵,又被我搞糊涂了? 两个人按着CRF同时输入数据。然后我要比较两个data sets
有什么不同,为了确保entry正确。
用了proc compare; by id; 我可以看到哪个obs的哪个variable不一样,然后重新看
CRF确认。但是我找不到哪个id在一个data 中有,另一个没有。
大家有什么好提议的?
A*******s
发帖数: 3942
14
那就是大傻根说的A-B和B-A
发信人: DaShagen (Unbearable lightness), 信区: Statistics
标 题: Re: 弱问SAS:如何找出两个data sets中missing的obs
发信站: BBS 未名空间站 (Thu May 27 09:27:19 2010, 美东)
lz means A-B and B-A i guess, (set subtraction),
can use sql except

sets
o****o
发帖数: 8077
15
check the manual:
OBS= data step option:
Specifies when to stop processing observations
1. Valid in: DATA step and PROC steps
2. Category: Observation Control
3. Restriction: Use with input data sets only <----********
4. Restriction: Cannot use with PROC SQL views
5. Default MAX

number of observations in output dataset, the following code doesnt work
D******n
发帖数: 2836
16
其实也没什么必要,就是有时候一个score value占了很多个obs, 一个score value占
两个bin,
用proc rank就合并了,其实也合理,可是有人就喜欢看着一个score value占两个bin。
l**********8
发帖数: 305
17
来自主题: Statistics版 - SAS macro Debug

************************************************************
上个完整的LOG, 谢谢大家了
****************************************************************
NOTE: PROCEDURE PRINTTO used (Total process time):
real time 0.00 seconds
cpu time 0.01 seconds

25
26 *-----------------------------------------------------;
27 %macro process(year,yr,file,first,numobs);
28
29 data cartemp;
30 set loc.carl&year.
31 (firsto... 阅读全帖
k*****u
发帖数: 1688
18
dont forget to check ERROR WARNING in the log file.
send out output to email.
merge/inner join/ left join/ right join, don't forget to check if they are
correct. check if the key variable should be unique or not.
Join a small table with a large table, index them .
join a large table with a large table(I once did a 3M table left join
1billion table), proc sql left join takes more than 3 days and not finished.
Then I tried split the 2M table to a small table with 5k obs each, then use
macro and lo... 阅读全帖
A******s
发帖数: 77
19
来自主题: MedicalCareer版 - @@@@ 2015 IM IV/Rej list@@@@
from USMLE-FORUMS as of last night. That thread is way better organized than
USMleforum.com As you can see, alot of folks with strong package arent
getting anything from anywhere.
http://www.usmle-forums.com/img-residency-match-forum/89681-201
Green card / Citizens:
usmleuk- 266/276/pass/2010/2m obsie/4 pubs/GC
Viper- 261/269/pass/2008/no USCE/3 pubs/MPH (US)/GC
sobored- 257/252/pass/2012/1.5 years clinical research, 31 pubs/Citizen
prtysha- 256/257/pass/2013/1m USCE/citizen - Reading, Univ Buff... 阅读全帖
a***r
发帖数: 420
20
来自主题: Statistics版 - 请教一个数据分类的问题
好好整理了一下思路,可以转化为这样的问题
我真正有的dataset是 list.dat和 LD.dat
list.dat格式:
V1
A
B
C
D
LD.dat格式
V1V2
A A
A C
B B
B D
C A
C C
D B
D D
1。从list.dat的第一个obs开始,将第一个obs存入一个新的database:independent.
dat
2。从list.dat中删除在LD.dat中与这个V1对应的所有V2的值(没有则只删除这个V1本
身)
3。删除后的database覆盖原list.dat
重复step1,2,3
直到list.dat的obs数为0
over
在这个举例的dataset里,A和C是同质的,B和D是同质的
从list的第一个obs开始
把A存入independent.dat,并把A,C从list中删去
然后剩下B D
从B开始,把B 存入independent.dat,并把B,D从list中删去
此时list obs为0
independent.dat 里存有A B
因为C和A是同质的,D和B是同质的
这两个obs即可代表这四个
杯具是,这个我也不会
a******n
发帖数: 11246
21
来自主题: Statistics版 - 请教一个bootstrapping的问题。
请教各位朋友一个关于bootstrapping方法本身的问题。
背景是这样的,N个observations,要fit一个model,假设model形式是Y=a+bX。
我理解的bootstrapping做法是,从这N个obs里re-sample N个obs, with replacement
,得到新的一组N obs,反复re sample,得到上千组N个的obs。每组fit一次model,最
后把上千个a和b平均一下。请问我这么做对么?
另外一种做法是,re-sample少量几次,比如10次吧,然后把这些re-sampled obs和原
来的N个合并在一起,这样一共有11*N个obs,直接fit model。这样做对吗?我总觉得
这么做和我理解的是完全不一样的,根本不能算bootstrapping,谁帮我confirm一下呀
s******s
发帖数: 2837
22
来自主题: Statistics版 - 请教两个SAS ADV问题
我觉得这两题是这样理解的。
第一题:D里面第二个dataset sum以后colume name变成了sum(Cost),out union corr
以后就出现了三列:id,cost,sum(Cost).A里面cost大小写确实不一样,不过sas识别
var是case insensitive,所以至少不是本题考查的点,不要在意这些细节。。。
第二题:multiple set option是依序读每个dataset的obs,读到最少obs的那个dataset
就不读了。所以只用set one;set two的话,读完第一个obs就不会继续读了,因为第二
个dataset只有一个obs。
用if _n_=1,可以读取并retain第二个dataset里面那个obs的值.PDV里面存有SumY并且
一直是36.当_n_=2的时候就没第二个dataset的事情了,只读第一个dataset的obs。所
以答案是D。
A*********u
发帖数: 8976
23
来自主题: Statistics版 - 问一道sas base题
first program
first obs -> PDV
then
looping starts
then
capital=5000 capital=10000 capital=15000
then
looping ends
then
PDV output to first obs of dataset
name="first", rate=0.0718, capital=15000
(the default output happens when SAS meets "run;")
then
2nd obs -> PDV. in PDV, name and rate are overwriten and
capital is retained b/c the sum function
then
looping starts
then
capital=20000 capital=25000 capital=30000
then
.
.
.
second program
looping starts
first obs -> PDV
capital=5000
second obs -
A*********u
发帖数: 8976
24
来自主题: Statistics版 - 问一道sas base题

first program
first obs -> PDV
then
looping starts
then
capital=5000 capital=10000 capital=15000
then
looping ends
then
PDV output to first obs of dataset
name="first", rate=0.0718, capital=15000
(the default output happens when SAS meets "run;")
then
2nd obs -> PDV. in PDV, name and rate are overwriten and
capital is retained b/c the sum function
then
looping starts
then
capital=20000 capital=25000 capital=30000
then
.
.
.
second program
looping starts
first obs -> PDV
capital=5000
second obs
p********1
发帖数: 1011
25
一共6个independent variables。从这个output上看,哪些variables是好的
predictors?
开始的那些note“!= 0 predicts failure perfectly”
是什么意思? 怎么解决?
logistic y i.var1 i.var2 i.var3 i.expect i.pwar i.open
note: 1.var1 != 0 predicts failure perfectly
1.var1 dropped and 14539 obs not used
note: 2.var1 != 0 predicts failure perfectly
2.var1 dropped and 4903 obs not used
note: 4.var2 != 0 predicts failure perfectly
4.var2 dropped and 96 obs not used
note: 4.var3 != 0 predicts failure perfectly
4.var3 dro... 阅读全帖
r*********8
发帖数: 4
26
来自主题: Statistics版 - 问一个SAS 的问题
我现在有2列char 的var, 我的目的是先把var中的每个“word" 分成一个个独立的
field, 然后看这两个var中有多少match的"word".
example:
原始数据:
var1 var2
looks like this looks like
disney disney cruise
----> 分解成:
var1_1 var1_2 var2_1 var2_2 var2_3
looks like this looks like
disney disney cruise
在var中每个“word"是用空格隔开的。我现在遇到的问题是,每个var含的"word”是不
一样的,可是我的做法好像就以最后一个obs的var的word个数来定loop的长度了。比如
上面的例子中,第一个obs,var1有2个word,第二个obs,var1只有一个word,结果do
loop的长度就成了1,然后就没有var1_2产生了。... 阅读全帖
f*r
发帖数: 3968
27
来自主题: Astronomy版 - 中国天文台一览表
中国天文台一览表
摘自洪韵芳主编的《天文爱好者手册》 扫描整理 覃 育
(如本文作者认为本站侵犯了版权,请马上通知本站,本站会在接到您的通知后的两天内
将本文撤下。)
台名
所在地
坐 标
海拔(米)
主要仪器设备
经度
h m
纬度
° '
上海天文台
Shanghai Obs.
上海市徐家汇
东经8 05.8
北纬31 12
95
1.56米天体测量望远镜,口径25米射电望远镜VLBI站,60厘米入卫激光测距仪,40厘米
双筒望远镜, Ⅱ型光电等高仪,氢原子钟等
上海市松江县余山 东经8 04.3 北纬31 06
云南天文台
Yunnan Obs.
云南省昆明市
东经6 51.1
北纬25 02
2002
1米RCC望远镜、太阳精细结构望远镜,SBG天文大地自动照相机,1.2米地平式激光测距
仪,60厘米望远镜。色球双筒望远镜,口径10米太阳射电望远镜等
陕西天文台
Shanxi Obs.
陕西省临潼
东经7 18.2
北纬34 57
497
长波授时台,短波授时台,水平式子午环,I型全自动光电等高仪等
北京天文台
Beijing Obs.
北京沙河工作站(共有沙河
p********a
发帖数: 5352
28
来自主题: Statistics版 - 请教sas code问题
the retain statement retains the values of y1/y2 once you define them.
Before _N_=25, y1=. y2=.
After _N_=25, y1=x(obs=25), y2=.
When _N_>=50, y1=x(obs=25), y2=x(obs=50)
so use _N_=50 to keep only one obs left.
Actually, you can choose any number >=50. eg: "if _N_=60;" or "if _N_=100;"
p********a
发帖数: 5352
29
☆─────────────────────────────────────☆
ifif (babe\'s babe Q~) 于 (Thu Mar 6 21:38:39 2008) 提到:
产生i从0到10,每次+0.1。将i=1.1的obs set到一个新data set b 里,可是log说没有
读到obs.将i=0.6的obs set到一个新data set c 里,这回读到一个obs了. 我用的SAS
9.0
为什么b里没有observation, c里有?
为什么b里没有observation, c里有?
data a;
do i=0 to 10 by 0.1;
output;
end;
run;
data b;
set a;
if i = 1.1;
run;
data c;
set a;
if i=0.6;
run;
☆─────────────────────────────────────☆
ifif (babe\'s babe Q~) 于 (Thu Mar 6 21:46:17 2008) 提到:
顶一个,为什么b里没有data,c里有?
a********a
发帖数: 346
30
来自主题: Statistics版 - help in R
I have a program as following,
beta1=2
beta2=8
#for (i in 1:2){
obs=2
x1=matrix(NA,obs,1)
x2=matrix(NA,obs,1)
for (g in 1:obs){
x1[g,]=runif(1,1,2)
x2[g,]=rnorm(1,1)
}
data=cbind(x1,x2)
data
#d=rbind(data[i])
#}
I run the program 2 times, and get the data like following each time,
> data
[,1] [,2]
[1,] 1.452696 0.2718456
[2,] 1.514587 1.2971293
> data
[,1] [,2]
[1,] 1.172726 -0.9674824
[2,] 1.159079 0.5483838
how can I combine these data by row, i.e. I want to get d
s*****n
发帖数: 2174
31
来自主题: Statistics版 - help in R
beta1=2
beta2=8
result <- NULL
for (i in 1:2){
obs=2
x1=matrix(NA,obs,1)
x2=matrix(NA,obs,1)
for (g in 1:obs){
x1[g,]=runif(1,1,2)
x2[g,]=rnorm(1,1)
}
data=cbind(x1,x2)
result <- rbind(result, data)
}
result
x***u
发帖数: 962
32
来自主题: Statistics版 - 问一个control chart的问题
比较郁闷 underlying data 不是normal 分布的, 而且一个时间点只有一个
observation,不知道该用哪种control chart了。其实这些observation也是算出来的, 所以连续的obs之间不是独立的。比如说,第一个obs是用1-5算出来的,第二个obs是用2-6算出来的,第三个obs是用3-6算出来的。实在想不出该用哪种chart了,而且要写成sas code,最好不要太复杂。哪位有经验 帮帮忙 谢谢!
R*******s
发帖数: 136
33
来自主题: Statistics版 - 问一个control chart的问题
移动平均?
minitab里面有

的, 所以连续的obs之间不是独立的。比如说,第一个obs是用1-5算出来的,第二个
obs是用2-6算出来的,第三个obs是用3-6算出来的。实在想不出该用哪种chart了,而
且要写成sas code,最好不要太复杂。哪位有经验 帮帮忙 谢谢!
t********y
发帖数: 469
34
来自主题: Statistics版 - how do you deal with sparse data?
sparse data
你应该给那些outcome不同的weight
再来run 那个proc logistic
比如说 yes的obs是1个
no的obs是100个
可以把no的outcome先group一下,只取10的obs
再给yes的weight=1
no的obs的weight =10
S******y
发帖数: 1123
35
来自主题: Statistics版 - nodupkey 删除的observations怎么查看?
%macro cmpar_d2(lib1=,mem1=,lib2=,mem2=);
/* Get the obs in lib1.mem1 but not in lib2.mem2 */
/* and store OBS and VARS in ONEXTWO Dataset */
proc sql;
create table onextwo as
select *
from &lib1..&mem1
except corr all
select *
from &lib2..&mem2
;
quit;
/* Get the obs in lib2.mem2 but not in lib1.mem1 */
/* and store OBS and VARS in TWOXONE Dataset */
proc sql;
create table twoxone as
select *
from &lib2..&mem2
except corr all
h********o
发帖数: 103
36
DATA ONE;
INPUT A @@;
CARDS;
5 4 7 8 2 1 6 8
;
DATA ONE;
SET ONE;
OBS + 1;
RUN;
PROC SORT DATA = ONE;
BY DESCENDING OBS;
RUN;
DATA ONE;
SET ONE;
B = LAG(A);
RUN;
PROC SORT DATA = ONE OUT = FINAL (DROP = OBS);
BY OBS;
RUN;
k*****u
发帖数: 1688
37
来自主题: Statistics版 - 请教一个SAS _n_的问题
为什么“if _n_=1 then do until (last)”这一段中要用if _n_=1来判断呢?如果不
用的话,sum_x sum_y都是可以正确求出来,但是obs就是missing
请教一下这儿‘if _n_=1’起的什么作用?
data temp;
input x y@;
datalines;
1 2
3 4
5 6
. 9
6 7
7 .
1 8
6 3
;
run;
data q3;
if _n_=1 then do until (last);
set temp nobs=obs end=last;
sum_x+x;
sum_y+y;
end;
set temp;
if x=. then x=sum_x/obs;
if y=. then y=sum_y/obs;
run;
s******o
发帖数: 283
38
来自主题: Statistics版 - 请帮忙解答两个SAS base考题
Q1, For this question, why the answer is C not A, the month is descending order
is not July,June and May, Why? Thanks very much for your kind help!
Given the SAS data set WORK.TEMPS:
Day Month Temp
--- ----- ----
1 May 75
15 May 70
15 June 80
3 June 76
2 July 85
14 July 89
The following program is submitted:
proc sort data=WORK.TEMPS;
by descending Month Day;
run;
proc print data=WORK.TEMPS;
run;
Which output is correct?
... 阅读全帖
t********m
发帖数: 939
39
来自主题: Statistics版 - 请教SAS中如何如果flag的问题
谢谢你这么快回复。我一开始也是这么写的,可是结果不对。按照这样写出来的结果是:
Obs obs v1 flag
1 1 1 0
2 2 1 1
3 3 1 0
4 4 2 1
5 5 3 0
6 6 3 0
7 7 4 1
8 8 5 1
而我想要的结果是:
Obs obs v1 flag
1 1 1 0
2 2 1 0
3 3 1 0
4 4 2 1
5 5 3 0
6 6 3 0
7 7 4 1
8 8 5 0
D******n
发帖数: 2836
40
Let me share with you some tricks or techniques i use to debug or facilitate
debugging SAS codes in Linux environment.
> proc print with obs= option. (ya very simple very naive)
> endsas statement.
> /* */ to comment out codes.
> sas -obs ### xxxx.sas (SAS options need not to be hard written in the
script, you can call it when batch running it)
> use put in data step to monitor SAS variable.
> use mprint, mlogic, symbolgen, source, source2 to debug macros.
> use %put to monitor macro variable.
>... 阅读全帖
b*********l
发帖数: 5
41
来自主题: Statistics版 - 请教版上高人一个SAS编程问题
我现在有一个sas data set像如下这样, 其中有三个char variable (id, A, B),数据
中空白的都是missing value.
obs id A B
1 a x1
2 x2
3 b y3
4 y4
5 y5
6 c x6
7 y7
现在我想把id这个column的missing value补完整。比如x2这个obs对应的id是a, y4 y5
这两个obs对应的id是b,y7这个obs对应的id是c。
请问有什么方法呢?多谢大家给点提示了!
c*******n
发帖数: 892
42
y是binary,positive/negative
x是一个categorical variable,有20个level,每个level有30-50个obs不等
现在想比较x小于40个obs的那组negative rate是不是比大于40个obs的那组大,或者是
说x每个level的negative rate会不会随着obs的增多而减小?
请教一下,能不能把x每个level的rate和frequency拿出来作为一个新的data,根据
frequency大于40小于40分组做ttest 或者直接rate和frequency做regression?
如果不能,那应该怎么做呢?
n*****5
发帖数: 61
43
来自主题: Statistics版 - sas base 题一问
27.Given the SAS data set WORK.TEMPS:
Day Month Temp
— —– —-
1 May 75
15 May 70
15 June 80
3 June 76
2 July 85
14 July 89
The following program is submitted:
proc sort data=WORK.TEMPS;
by descending Month Day;
run;
proc print data=WORK.TEMPS;
run;
Which output is correct?
A.
Obs Day Month Temp
— — —– —-
1 2 July 85
2 14 July 89
3 3 June 76
4 15 June 80
5 1 May 75
6 15 May... 阅读全帖
s*****r
发帖数: 223
44
来自主题: Statistics版 - 最大值和最小值输出
刚学sas不久, 需要写macro 输出指定var 的最大值或最小值的obs(可能不只一个)。
exampledata
obs var1 var2
1 3 2
2 9 9
3 4 6
4 9 2
我知道可以
* sort to find the max or min;
proc sort data = exampledata out = sortedexampledata ;
by ascending/descending var1/var2;
run;
*determine how many obs to keep in the data set;
data sortedexampledata;
set sortedexampledata ( obs = N);
run;
现在的问题我不知道怎么coding 找到 这个 N, 我知道 N 就是 var1 or var 2 的
最大值或者最小值的 frequency, 怎么 输出这个frequency 然后传递给 N。
谢谢
s*****r
发帖数: 223
45
来自主题: Statistics版 - 最大值和最小值输出
刚学sas不久, 需要写macro 输出指定var 的最大值或最小值的obs(可能不只一个)。
exampledata
obs var1 var2
1 3 2
2 9 9
3 4 6
4 9 2
我知道可以
* sort to find the max or min;
proc sort data = exampledata out = sortedexampledata ;
by ascending/descending var1/var2;
run;
*determine how many obs to keep in the data set;
data sortedexampledata;
set sortedexampledata ( obs = N);
run;
现在的问题我不知道怎么coding 找到 这个 N, 我知道 N 就是 var1 or var 2 的
最大值或者最小值的 frequency, 怎么 输出这个frequency 然后传递给 N。
谢谢
t**********y
发帖数: 374
46
I did some test runs using built-in dataset and got different outputs:
1. ===
> test=lmer(Reaction~Days+(Days|Subject),sleepstudy)
> summary(test)
Linear mixed model fit by REML ['lmerMod']
Formula: Reaction ~ Days + (Days | Subject)
Data: sleepstudy
REML criterion at convergence: 1743.6
Scaled residuals:
Min 1Q Median 3Q Max
-3.9536 -0.4634 0.0231 0.4634 5.1793
Random effects:
Groups Name Variance Std.Dev. Corr
Subject (Intercept) 612.09 24.740
... 阅读全帖
f*r
发帖数: 3968
47
来自主题: _Astronomy版 - 中国天文台一览表
中国天文台一览表
摘自洪韵芳主编的《天文爱好者手册》 扫描整理 覃 育
(如本文作者认为本站侵犯了版权,请马上通知本站,本站会在接到您的通知后的两天内
将本文撤下。)
台名
所在地
坐 标
海拔(米)
主要仪器设备
经度
h m
纬度
° '
上海天文台
Shanghai Obs.
上海市徐家汇
东经8 05.8
北纬31 12
95
1.56米天体测量望远镜,口径25米射电望远镜VLBI站,60厘米入卫激光测距仪,40厘米
双筒望远镜, Ⅱ型光电等高仪,氢原子钟等
上海市松江县余山 东经8 04.3 北纬31 06
云南天文台
Yunnan Obs.
云南省昆明市
东经6 51.1
北纬25 02
2002
1米RCC望远镜、太阳精细结构望远镜,SBG天文大地自动照相机,1.2米地平式激光测距
仪,60厘米望远镜。色球双筒望远镜,口径10米太阳射电望远镜等
陕西天文台
Shanxi Obs.
陕西省临潼
东经7 18.2
北纬34 57
497
长波授时台,短波授时台,水平式子午环,I型全自动光电等高仪等
北京天文台
Beijing Obs.
北京沙河工作站(共有沙河
I******a
发帖数: 3812
48
http://www.wenxuecity.com/news/2011/09/03/1461780.html
2010年3月26日晚8时22分许,韩国“天安”号警戒舰在正常巡逻至白翎岛以南水域(这
里有朝韩经常发生冲突的“北方限界线”)时遭袭击沉没。根据两个月后韩国联合调查
团发布的最终结果,认定攻击源来自秘密渗透到韩国水域的朝鲜鲑鱼级潜艇,韩国因此
中断与朝鲜的一切交流。
韩国海军反潜警戒舰“天安”号沉没事件过去一年后,韩国《统一新闻》最新爆料称,
这起悲剧的“主犯”其实是美国,而造成“天安”号爆炸沉没的正是被美军秘密邀请来
参加美韩演习的以色列潜艇。
图为“天安”号舰首残骸被打捞出水。
图为“天安”号舰首残骸被打捞出水。
图为“天安”号舰首残骸被打捞出水。
韩方称被吊出水面的“天安舰”舰首呈现遭受强力冲击的痕迹(韩联社)
神秘鱼雷 来自何方的神圣
韩国联合调查团发布的这份调查结论存在诸多疑点,韩国国防部又一直遮遮掩掩,对民
众的质疑不予解释。韩国民主党议员朴英善指出,根据韩国联合参谋本部的作战序列,
属于韩国海军第二舰队的“天安”号平时是不准进入白翎岛以南水域的,因为那里平均
水深... 阅读全帖
s****i
发帖数: 243
49
来自主题: NextGeneration版 - 剖腹产后11个月怀孕,真的很危险么?
I am pregnant with 2nd baby 6 months after C-section. I asked three
different obs about the risk, they all told me there is no problem as long
as the second delivery is a planned C. One of the obs even suggested that if
I want, I can try VBAC. The other one told me that within more than 30
years of her practice, she has never encountered any mother who had
complications just because of a second pregnancy within one year of the last
one.
Please kindly help your sister-in-laws to seek more obs opi... 阅读全帖
b*s
发帖数: 82482
50
来自主题: LeisureTime版 - 本年度最让我期待的好莱坞电影
amour:
( / m (r)/ , Fr. /a mur/ )
Forms: Also ME amoure, 15–16 amor.
Etymology: < Old French amur, amour < Latin amōr-em love, < amā-re... (
Show More)
Thesaurus
1. gen. Love, affection, friendship. Obs.
c1300 K. Alis. 4573 Alisaunder..wolde him, with gret honour, Have y-
fonge in his amour.
c1330 Florice & Bl. 521 Tho spak Clarice to Blauncheflour Wordes ful of
fin amour.
c1386 Chaucer Knight's Tale 297 For par amour I loved hire first or thou.
1660 R. Burney Κ ρδι... 阅读全帖
1 2 3 4 5 6 7 8 9 10 下页 末页 (共10页)