R_Language 板


LINE

我的测试档案 ----------------tmp.txt------------------ name,time,price A,2000/01/31,-1 A,2000/02/29,2 A,2000/03/31,-3.3 A,2010/12/31,5 E,2000/01/31,-4 E,2000/02/29,5 E,2000/03/31,-1.1 E,2010/12/31,5 Z,2000/01/31,-10 Z,2000/02/29,7 Z,2000/03/31,-1.8 Z,2010/12/31,-2 ---------------end of file---------------- R code: x=read.table("tmp.txt",sep=",", header=TRUE) x[,3] = abs(x[,3]) x.wide=reshape(x,direction="wide", v.names=c("price"), timevar="name", idvar="time") output: time price.A price.E price.Z 1 2000/01/31 -1.0 -4.0 -10.0 2 2000/02/29 2.0 5.0 7.0 3 2000/03/31 -3.3 -1.1 -1.8 4 2010/12/31 5.0 5.0 -2.0 matlab code ---------------new.txt------------- name time price A 2000/01/31 -1 A 2000/02/29 2 A 2000/03/31 -3.3 A 2010/12/31 5 E 2000/01/31 -4 E 2000/02/29 5 E 2000/03/31 -1.1 E 2010/12/31 5 Z 2000/01/31 -10 Z 2000/02/29 7 Z 2000/03/31 -1.8 Z 2010/12/31 -2 ------------------------------------ new.txt之间是tab,是空白的话,下面'\t'请改成' ' Dat = readtable('new.txt', 'Format', '%s%s%f', 'delimiter', '\t'); Dat.price = abs(Dat.price); W = unstack(Dat, 'price', 'name'); Matlab的速度很快,可以试试看。 --



※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 36.238.81.92
1F:→ goodparent1:有办法用r跑吗,MATLAB 档案很大 很难下载。 11/01 17:44
R的code就在上面,只是用sep=",",如果你的间隔是空白就改成sep=" " ※ 编辑: celestialgod 来自: 218.164.79.5 (11/01 19:09)
2F:→ goodparent1:我是指负数转正数的code 如果我要用r跑 该怎麽写 11/02 11:45
3F:→ goodparent1:不好意思,我看到了。 11/02 11:47
4F:→ goodparent1:可是我用ABS跑,只能跑到1800多列,剩下的档案会遗失 11/03 22:41
5F:→ celestialgod:我没那麽大的档案 我没法测试= = 11/03 23:17
6F:→ goodparent1:https://www.dropbox.com/s/6chu02wylo3owmq/test.txt 11/04 12:11
我改成这样,结果如下: x=read.table("test.txt", header=TRUE) x[,3] = abs(as.numeric(x[,3])) x.wide=reshape(x,direction="wide", v.names=c("Adjprc"), timevar="name", idvar="time") > dim(x.wide) [1] 504 1986 ※ 编辑: celestialgod 来自: 36.239.250.66 (11/04 19:03)
7F:→ goodparent1:可是这样转出来,我的时间列就不见了.... 11/05 20:22
8F:推 Wush978:push 11/07 22:17
9F:推 Wush978:my code: 11/07 22:18
10F:→ Wush978:x=read.table("test.txt", header=TRUE) 11/07 22:18
11F:→ Wush978:x$time=as.Date(as.character(x$time), format="%Y%m%d") 11/07 22:18
12F:→ Wush978:x$Adjprc=abs(x$Adjprc) 11/07 22:19
13F:→ Wush978:reshape部分相同, 最後得到如: `head(x.wide[,1:5])` 11/07 22:19
14F:→ Wush978:第一行就是时间(日期) 11/07 22:19







like.gif 您可能会有兴趣的文章
icon.png[问题/行为] 猫晚上进房间会不会有憋尿问题
icon.pngRe: [闲聊] 选了错误的女孩成为魔法少女 XDDDDDDDDDD
icon.png[正妹] 瑞典 一张
icon.png[心得] EMS高领长版毛衣.墨小楼MC1002
icon.png[分享] 丹龙隔热纸GE55+33+22
icon.png[问题] 清洗洗衣机
icon.png[寻物] 窗台下的空间
icon.png[闲聊] 双极の女神1 木魔爵
icon.png[售车] 新竹 1997 march 1297cc 白色 四门
icon.png[讨论] 能从照片感受到摄影者心情吗
icon.png[狂贺] 贺贺贺贺 贺!岛村卯月!总选举NO.1
icon.png[难过] 羡慕白皮肤的女生
icon.png阅读文章
icon.png[黑特]
icon.png[问题] SBK S1安装於安全帽位置
icon.png[分享] 旧woo100绝版开箱!!
icon.pngRe: [无言] 关於小包卫生纸
icon.png[开箱] E5-2683V3 RX480Strix 快睿C1 简单测试
icon.png[心得] 苍の海贼龙 地狱 执行者16PT
icon.png[售车] 1999年Virage iO 1.8EXi
icon.png[心得] 挑战33 LV10 狮子座pt solo
icon.png[闲聊] 手把手教你不被桶之新手主购教学
icon.png[分享] Civic Type R 量产版官方照无预警流出
icon.png[售车] Golf 4 2.0 银色 自排
icon.png[出售] Graco提篮汽座(有底座)2000元诚可议
icon.png[问题] 请问补牙材质掉了还能再补吗?(台中半年内
icon.png[问题] 44th 单曲 生写竟然都给重复的啊啊!
icon.png[心得] 华南红卡/icash 核卡
icon.png[问题] 拔牙矫正这样正常吗
icon.png[赠送] 老莫高业 初业 102年版
icon.png[情报] 三大行动支付 本季掀战火
icon.png[宝宝] 博客来Amos水蜡笔5/1特价五折
icon.pngRe: [心得] 新鲜人一些面试分享
icon.png[心得] 苍の海贼龙 地狱 麒麟25PT
icon.pngRe: [闲聊] (君の名は。雷慎入) 君名二创漫画翻译
icon.pngRe: [闲聊] OGN中场影片:失踪人口局 (英文字幕)
icon.png[问题] 台湾大哥大4G讯号差
icon.png[出售] [全国]全新千寻侘草LED灯, 水草

请输入看板名称,例如:BuyTogether站内搜寻

TOP