p********0
发帖数: 186 | 1 OK,
let me ask the same question the other way around
Why the average interest rate over 2 years is the simple average of
short term rate over(time0, time1, time2)???
If I understand correctly, interest rate at time0 are from now till 1 year.
interest rate at time1 is from 1 year till 2 year.
So if I want to get the average interest rate over 2 years horizon, why
isnot it the
average of (r(time0)+t(time1))? |
|
t******l
发帖数: 32 | 2 Why can't you treat X＿TIME1and Y＿TIME2 as time varying covariate with X＿
TIME1=0 for t>TIME1 and Y＿TIME2=0 for t |
|
p*******r
发帖数: 1951 | 3 首先要明确你要做什么比较？看你说的意思你的null hypothesis是不是
H0: treatment effect = time2 - time1 = 0 ?
如果是，那你的dependent variable 完全可以用 effect = time2 - time1。而不是单
独的time1, time2。
其次，你的数据显示不同的subject有完全不同的时间点， 看你的SAS程序里面也根本
没有用到时间变量week。所以你想做的不是longitudinal data analysis？你只是想比
较不同treatment, pre-, post- response的差别。这样的话用repeated measure
ANOVA就成了。
很久没有用过mixed model了。说的不对的地方大家指教。 |
|
m*****r
发帖数: 51 | 4 我有observations with different starting points. 比如， 第一个人1月开始赚钱
， 第二个人2月开始赚钱， 第三个人4月开始赚钱。 现在我需要把所有的人在第一个
月赚的钱放到一个variable （time1)下， 第二个月（2月，3月，5月）的放到同一个
variable (time2)下, 依次类推， 对每个人都要看开始赚钱后接下来的24个月。 有没
有什么办法create time1－time24. 如果24个时间点一个个的来， 太慢了。
例子来了：
原data
ID Jan Feb March April May Jun ....Dec
1 20 30 40 50 0 50
2 10 15 20 25 30
3 30 30 40
我需要的是
ID time1 time2 time3 time4 ...tiimeN
1 20 30 40 50
2 10 15 ... 阅读全帖 |
|
发帖数: 1 | 5 47 subjects
每个subject测量了2次（time1 and time2)，想知道这2次的测量值有没有difference
我知道最简单的就是跑个paired t test
问题的关键是这47个subjects的measure是在三个不同的地方采集的（采集的protocol还
不太一样），而且年龄，性别分布都不一样，
所以我能想到的是
Y=Y0+age+sex+site_index来remove掉 age/sex/site_index的影响
但是之后我就不知道怎么来比较time1 vs. time2
我之前想的是拿Y=Y0+age+sex+site_index这个fitting之后的residual值来
做paired t test,但是time1的fitting和time2的fitting不一样啊，能这么比较吗？
思路有点乱，哪位朋友帮忙理一下，谢谢
祝新年快乐！ |
|
g****t
发帖数: 31659 | 6 我十几年前MS Word的VBA很熟。花了不少时间把上海市政府一个部门所有文件做了分析，
提取出来名词术语。尽量让“下岗”这样的词的用法无矛盾。干的多了，我发现写程序
其实应该用word写。但是word虽好，却不支持这个功能。
好在emacs／vim支持。
对编辑器进行程序化操作，可以人机深度融合，改造working flow，提升编程的维度。
编辑器最后其实就是meta－meta－meta－。。。－programming
我举个最简单的例子。例如你有个python pandas s，其中一列为time。
现在想转成一个数组名字为time。
time=np.array(s['time'])
现在想几十列都这么办，怎么最快的做到？在python的空间里这是不容易的。
和之前我提问过的那个找到变量名字的问题有点关系。不是熟练工恐怕搞不定。
而且还容易出bug难维护。python 2 升级这种事再来一次，程序说不定就废了
但是在org－mode里，上面一段程序自动生成一长条程序，几十个
time=np.array(s['time'])
time1=np.array(s['time1'])
t... 阅读全帖 |
|
m******r
发帖数: 1033 | 7 最简单的方法生成如下格式
time=np.array(s['time'])
time1=np.array(s['time1'])
time2=np.array(s['time2'])
只需在excel里第一列先输入1,然后输入A1 +1, 然后下拉， 第二列用concat函数写你
要的东西，然后下拉，最后copy ,paste到编辑器即可。 |
|
n******g
发帖数: 2201 | 8 哈哈 可以 但是如果你在server 上干活 没有excel可用
这种情况常见
[在 magliner (magliner) 的大作中提到：]
:最简单的方法生成如下格式
:time=np.array(s['time'])
:time1=np.array(s['time1'])
:time2=np.array(s['time2'])
:只需在excel里第一列先输入1,然后输入A1 +1, 然后下拉， 第二列用concat函数写你
:要的东西，然后下拉，最后copy ,paste到编辑器即可。 |
|
D******6
发帖数: 6211 | 9 我现在的原始数据是一个数值变量，一个相对SAS的60年1月1日的妙数，巨大的数字。
我现在想要把它转换成真正的时间，有无现成函数可以把这个变量转换成年月日小时分
钟妙？
比如我的原始数据是31536001，就是说距离1960年1月1日0点0分0妙是31536001妙，也
就是1961年1月1日0点0分1妙。
我希望 time1= 31536001, time2= 01-01-1961 00:00:01
我现在只有time1，没有time2。有无现成的函数可以生成time2？
谢谢解答，包子答谢。 |
|
m******3
发帖数: 62 | 10 我用SAS run
proc mixed data=mix method=reml noclprint;
class ID treatment;
model score=treatment treatment*time1 treatment*time2
/solution noint;
random intercept time1 time2/subject=id;
repeated/subject=id;
weight wgt;
run;
结果DF 为0， 无法得到p-vlaue. 那位大牛给s说说，问题在哪呢？
谢谢 |
|
l****u
发帖数: 529 | 11 If time1 and time2 are random effects, treatment*time1 and treatment*time2
should also be random effects. Then both of them should
be put in the random statement. |
|
S******J
发帖数: 30 | 12 1. Most of the time you should not use Random and Repeated cards at the same
time.
2. You only need time1 and time2 in the Random card, not treatment*time1 and
treatment*time2.
3. You may want to specify the variance structure as well. |
|
j******o
发帖数: 127 | 13 Not fancy, but should works.
------------------------------------
data have;
input time $;
datalines;
10175996
10175925
;
run;
data obtain;
set have;
if input(substr(time, 7,2), best12.)<=30 then time1=substr(time, 1, 6);
else do ;
x=input(substr(time, 1,6), best12.)+1;
time1=put(x,z6.);
end;
run;
------------------------------------ |
|
m***c
发帖数: 118 | 14 data a;
input ID $ jan feb mar apr may jun jul aug sep oct;
cards;
1 20 30 50 0 50 40 20 70 55 25
2 . . 30 50 0 90 40 0 40 20
3 . . . 60 20 0 50 10 0 40
;
data b(keep=ID time1-time10); set a;
length tot $ 100;
array am(10) jan--oct;
do i=1 to dim(am);
if am(i)<=0 then am(i)=.;
tot=compress(catx('',tot,am(i)),,'p');
end;
array at(10) time1-time10;
do j=1 to 10;
at(j)=scan(tot,j);
end;
run; |
|
m******a
发帖数: 84 | 15 there is a data flow, find the median of the data in a given time interval.
有什么好的算法么？比如数据不停的流进来，然后给N查询区间[time1, tim2]，求区间
中位数 |
|
L*******y
发帖数: 23 | 16 我是186.00 买的Tesla Motors Inc (TSLA) 50股
Real-Time: 197.22 at 2/10/2014 1:55:30 PM -- Source BATS Exchange
Last Price197.22
Today's Change10.69 (5.73%)
Quote Time1:55 PM
Volume505,718
Open188.67
Bid197.13
Ask197.44
IndustryCONSUMER GOODSDay's Range188.55 - 199.29
52 wk Range34.25 - 194.50
Market Cap$23B
P/En/a
EPS-3.69
Div & Yieldn/a - (n/a)
Div Last Daten/a
Beta1.24 |
|
|
|
|
b*****e
发帖数: 853 | 20 Gooog已经准备好获胜演讲稿。
那这里拷贝粘贴一下，截至目前的投票情况。是不是马甲，是不是鱼版ID，是不是鱼版
常驻ID，大家自己识别吧。
竞选谁上谁下，怎么计票，那是老刑说了算。这里就是拷贝粘贴一下。
192/6742 GoooG
[申请]一天到晚上钩的鱼(Fishing) 版主(375b)
楼层失踪 楼层 鱼版常驻ID 鱼版ID ID IP Vote1 TIME1
2 zillionaire 24 支持 Wed Mar 12 15:31:19 2014
3 vanfossil 209 支持 Wed Mar 12 15:32:18 2014
4 kz80 184 支持 Wed Mar 12 15:36:53 2014
5 mako 168 灌水帖 Wed Mar 12 15:36:56 2014
6 LPisMM 75 支持 Wed Mar 12 15:52:36 2014
7 panangler 198 支持 Wed Mar 12 15:56:37 2014
8 hotlood 132 支持 Wed Mar 12 15:56:54 2014
9 ... 阅读全帖 |
|
b*****e
发帖数: 853 | 21 倒票是竞选规则允许的吧。不论是非，纯粹选举，很难。一想起阿扁肚皮上的枪伤，我
就一脑门浆糊。比起那些来，这里的江湖平津和气正气多了。Anyway，看看另一边的情
况。
Layer ID IP TIME1 给Willybass的投票 给gooog的投票
2 Dreamerx 107 Sun Mar 9 12:42:50 2014 支持 #N/A
3 CatchGodLine 108 Sun Mar 9 12:43:47 2014 支持 No
4 sonatina99 198 Sun Mar 9 13:10:28 2014 支持 #N/A
5 julia1223 71 Sun Mar 9 13:11:42 2014 支持 #N/A
6 XBaldwin 68 Sun Mar 9 13:12:55 2014 支持 No
7 VVV555 24 ... 阅读全帖 |
|
n********6
发帖数: 1511 | 22 第一次学在96年。大概part time1年，每周2小时，跟室友学，看书。
然后中断。
最近两周又在网上下。 |
|
A********2
发帖数: 546 | 23 Strong MM, you have a beautiful mind and graceful heart. Bless your Dad, he
IS thankful than you may think, for HE IS THERE watching you ALL THE TIME1.
YOU can live your life peacefully and thankfully. Treat your husband as
tender as possiple, for that is YOUR LIFE. |
|
s**********r
发帖数: 28 | 24 1．
那天在麦当劳吃饭，碰到张鸣剑，说了两句话就彼此很不快。这个狗嘴里吐不出象牙的人
，张嘴就说：“凌安你是我所见到过的人中运气最好的，跟了Steve那个超级混混混，日
子真是轻松，只要会吹牛老板就开心。”我心里很不淑女的骂了一句娘。心想你要说我是
一个小混混就直说好了，不用这么拐弯抹角的，你这个该死的娘娘腔。我迅速的把这个结
论扩大到广大的留美男生的群体，并且毫不留情的把跟我八百辈子都没有关系的前男友顺
遍再骂了一遍。
我是这样一个人，我可以说自己长得不好看，不聪明，喜欢口不择言，我想这也是实情，
但是并不等于别人可以说或者可以附和我骂自己的话。我当时的反映就是看着他笑，一边
笑一边在心里骂他白痴，他被我盯的发毛。我就一个字一个字的说：“嗯，这至少说明我
的口语比你的好，找工作没有问题。”他的脸色就变得又青又白的，十分的多姿多彩。这
是他的痛处，张鸣剑一年半硕士就毕业，找了一个工作一直都很顺利，可是onsite
interview的时候，他的一口结巴英语把hr给吓坏了，两个人对着说了半个小时，居然发
现彼此都听不懂对方讲的是什么。在最初的几分钟hr以为是自己得了语言混乱症，等到她
清 |
|
c**t
发帖数: 2744 | 25 awstats doesn't track clickstream.
For example:
Visitor A:
[time1] page 1 (entry page)
[time2] page 2
...
[timeN] page N (exit page)
awstats only count how many entry pages and exit pages, no break down.
Do you know any of software to track the breakdown? |
|
s********k
发帖数: 6180 | 26 能不能直接对JSON中的field做query，比如我的表里面有field是一个JSON的stream
[{time1:value1},{time2:value2}...],希望对特定的time或者特定比如大于某个值做
query |
|
s********k
发帖数: 6180 | 27 比如我把一个JSON的stream存在MongoDB
[{time1：value1},{time2:valu2}...]
CREATE TABLE json_test (
id serial primary time,
data jsonb
);
然后做
SELECT * FROM json_test WHERE Value>threhold?
这样的可以 |
|
c***c
发帖数: 21374 | 28 【 以下文字转载自 Internet 讨论区 】
【 原文由 cynic 所发表 】
用aw－stat，日志格式总是设定不对
日志如下
61.155.235.238 - - [17/Nov/2003:01:57:50 -0500] "POST /cgi-bin/regcheck.pl HTT
P/1.1" 200 228 etsscore.com "http://etsscore.com/regcheck.htm" "Mozilla/4.0 (c
ompatible; MSIE 6.0; Windows 98)" "-"
我这么写错了
%host %other %time1 %methodurl %code %bytesd %virtualname %refererquot %uaquot
%other
请问应该怎么写呀？》
谢谢 |
|
p********0
发帖数: 186 | 29 Hi,
I am using simulation to estimate the swaption value.
I have 500 unique paths of short rate over the next 10 years.
I am trying to use the Black model the valuate the 1 year into 1 year
european option.
In order to use the Black model, I need to figure out the Forward swap rate
volatility.
For each path, The forward swap rate can be backed from FSR(1, 2) = D(0, 1)
- D(0, 2)/Annuity(0, 1, 2).
Since I have short rate for time1, time2, so I can easily get the D(0,1)=EXP
(-AVG(r0, r1)*1). So I c |
|
h******a
发帖数: 198 | 30 > print(summary(lme4_pre_wm1))
Linear mixed model fit by REML
Formula: log(fat) ~ time + time1 + gender + per_income + weight + t2 +
(time + gender + weight | subject)
Data: nurt_miss_n1989
AIC BIC logLik deviance REMLdev
117809 117973 -58886 117673 117773
Random effects:
Groups Name Variance Std.Dev. Corr
subject (Intercept) 3.7478e-01 0.6121924
time 4.8376e-04 0.0219946 -0.216
gender1 |
|
g******c
发帖数: 116 | 31 thanks,
问题就是 X＿TIME1=０（第２点）会有解释上的歧义。X＿TIME１在 TIME２就不存在（
不是缺失而是没有发生），如果给它一个０值，那么就是说它在TIME２有值并且为０，
这个不合情况。
我不知道有什么模型能考虑到变量只发生某时间点的情况 |
|
P******e
发帖数: 75 | 32 我们从Literature里收集了很多关于牛奶中18种氨基酸的含量的数据，每个数据点是文
献当中的平均值。不同文章的牛奶在不同时间收集，不同地方来源.
数据整理后,我们想看一下不同因子有没有影响。
从三个因素考虑，分别是time,region,term.数据是不平衡的.
下表标出了其中两个因素。每个格子里面是多少个观测点。可见是很不平衡的。
Region1 Region2 Region3 Region4
time1 . 6 2 1
time2 . 7 2 3
time3 1 8 7 5
time4 1 11 4 5
time5 3 10 . 2
我们有几个问题请教：
1. 由于我们的每个数据点是每篇文献当中的平均值,能用ANOVA么?数据好像不是很正态
分布。要检测Assumptions?
2. 对我们这个unbalanced的数据处理,可以做Three-way ANOVA么?
3. 用proc GLM, Model AA1 AA2=tim |
|
d*******o
发帖数: 493 | 33 /*SET UP A RESULT DATASET*/
data result;
do id = 'a', 'b', 'c' ;
output;
end;
run;
/*SPECIFIC INTERVALS ARE ASSIGNED*/
%macro summary2(time1, time2, time3, time4);
%do i=1 %to 4;
proc sql;
create table result as
select a.*, b.first&&time&i
from result as a, (
select id, sum(sale) as first&&time&i
from test
where day le &&time... 阅读全帖 |
|
c**********5
发帖数: 653 | 34 First all,if if time1 and time2 in the interaction term,shoud it in the
model first? |
|
l****u
发帖数: 529 | 35 大侠，能给解释一下吗？
他把time1 and time2 设为random effect，那么treatment*time也应当为random
effect
呀，或者至少不能放在model中啊。
多谢了
same
and |
|
d**********r
发帖数: 24123 | 36 你的 treatment*time1 或者 treatment*time2 只有一个对应的sample了吧。
每个effects没有>=2 的 repition, 你的DF 就会为0了。 |
|
s*****e
发帖数: 157 | 37 在以下dataset中，每个id在不同的星期可以重复参与不同的treatment（A或B),最后参
与的模式包括A,B, AB,AA,BB,ABB,AAB等等，同时每个id在同一个tr下分别在time1和
time2测量了一次.所以一个id可能有多个repeated measure.
我想建个mixed model比较tr对response的影响
proc mixed;
class id tr time;
model response=tr|time;
repeated /subject=id;
run;
程序无法运行，可能因为同一个id在同一个time有多个observation,这种情况怎么处理
呢？请大家支招，多谢！
dataset:
id week tr time response
01 1 A 1 xx
01 1 A 2 xx
01 2 A 1 xx
01 2 A 2 xx
01 3 B 1 xx
01 3 B 2 xx
02 1 ... 阅读全帖 |
|
s*****e
发帖数: 157 | 38 time1是pre, time2是post,所以要比较。week是随机对应两种tr的，我不知道还要不要
考虑？
structrure随便给了一个。error具体不记得了，但似乎是由于多个repeated measure
导致的 |
|
s*****e
发帖数: 157 | 39 time1 和 time2 各一次，这个不是导致错误的重复阿 |
|
j*******y
发帖数: 58 | 40 you have the hell 2 weeks at time1 and 2 weeks at time2. look at the hell
subject 4. |
|
s*****e
发帖数: 157 | 41 但是我要比较time1(pre)和time2(post)response的不同，如果time和week混起来，怎
么比较呢？ |
|
j******4
发帖数: 6090 | 42 variables: Age(in month) Female(dummy), Smoke(dummy)
outcome(counts): time1 time2 time3 time4
Reshaped into long form then the data looks like this:
ID time Y Age Female Smoke
1 1 0 30 1 1
1 2 1 30 1 1
1 3 5 30 1 1
1 4 2 30 1 1
...
可以用GEE或者mixed model，但是据我观察，在 time=1 的时候，会有很多 Y = 0，我
在想 Zero-Inflated-Poisson(ZIP) 是否可以用来分析这种repeated measures的data，要
如何用呢？
多谢各位~ |
|
w******8
发帖数: 59 | 43 Dear all, I am running this procedure for time dependent covariates. I also
have two time INDEPENDENT covariates: Age at baseline which is continuous;
HIVPatient ( a status indicating whether having HIV) which is binary.
The following codes works fine without adding HIVPatient, but I don't
understand why there is nothing estimated for HIVPatient once I included it.
I think this is weird. Could someone help me?
Thanks!
proc phreg data = OnePer;
class HIVPatient(PARAM=Reference Ref='0');
mod... 阅读全帖 |
|
s**y
发帖数: 358 | 44 这是目前想到的方法，但觉得不够精炼，有没有更好的办法？
求大牛出来指点一下！
data one;
input @1 subject @5 date mdyampm25.2;
format date mdyampm25.2;
datalines;
1 01/08/2014 08:21:00 AM
1 01/08/2014 11:44:00 AM
2 01/11/2014 04:57:00 AM
2 01/11/2014 09:24:00 AM
3 01/07/2014 03:41:00 PM
3 01/07/2014 04:52:00 PM
...
;
format mdyampm25.2;
run;
data two;
input Subject Measure @13 Time mdyampm25.2;
format time mdyampm25.2;
datalines;
Subject Measure Time
1 10.50 01/07/2014 02:27:26 PM
1 7.10 01/08/2014 ... 阅读全帖 |
|
s*********h
发帖数: 6288 | 45 SAS就是处理这种东西的时候特别麻烦。
你的原始数据的列数应该有最大值吧？比如说是N
那建一个array ORG {N} 把原来的那些列都扔进去。
然后建一个 array NEW {24} time1-time24
设一个indicator A来标记是不是开始算钱了
A=0;
k=1;
do i=1 to 24;
if A = 0 and org[i] ^=. then do;
new[1]=org[i];
A = 1;
k + 1;
end;
else if A = 1 and k <= 24 then do;
new[k] = org[i];
k + 1;
end;
end;
然后keep 这个array new好了。
大致上思路差不多。前提是如果开始赚钱了后面有一个月没赚钱也得算进去
|
|
g*********d
发帖数: 24 | 46 一家小破公司，开始没当回事，后来发过来几道题，都sql，让做。
然后给跪了，很受打击！
最后一道题到现在没想明白是不是不用cursor根本就做不出来。
id1 time1 1
id2 time2 1
id3 time3 1
id4 time4 2
id5 time4 1
id6 time5 2
id7 time6 2
....
求所有同样数字连续超过三次的纪录。大概就这个意思吧，比如上面的例子中，1就是
符合条件的，因为连续时间段内出现了三次，但是2就不行，因为三个2中间插了一个1. |
|
a****k
发帖数: 117 | 47 If using mysql, we can use self-join.
id, time, num
id1 time1 1
id2 time2 1
id3 time3 1
SELECT DISTINCT (a.num) FROM table a JOIN table b ON a.num=b.num JOIN table
c ON a.num=c.num WHERE a.id-b.id = 1 and b.id-c.id=1 |
|
