Redian新闻
>
新手请教一个 R 作图的问题
avatar
d*8
2
住在SAN GABRIEL, CA. 想给大宝小宝找一个好的钢琴老师。 现在的钢琴老师太nice
, 孩子没啥
压力, 就是离家近。 想找个有经验, 严格点的。
听说很多老师每年都会安排学生参加一些表演, 这样家长也能了解, 老师和孩子都有
目标。 我们家俩
宝快5岁了, 跟现在的老师学了半年了, 能看懂简单的五线谱, 但手型和节奏很不规
矩。
听国内一教钢琴的人说, 4岁孩子3个月应该能弹曲子了, 我们的好像有点慢呀!
大家都有什么经验呀?
如有好的老师, 请推荐!
先谢谢了!
avatar
c*7
3
【 以下文字转载自 Joke 讨论区 】
发信人: linvk (aglee), 信区: Joke
标 题: 莫言最伤人的83句话,非常经典
发信站: BBS 未名空间站 (Fri Oct 12 01:38:46 2012, 美东)
看了不准哭
avatar
h*t
4
我刚刚开始自学R,作图的时候遇到了一些问题。可能很基本,但是搜索了半天也没有找
到答案。这里高手多,请大家帮我看看。
我需要做一个 scatterplot 的matrix,找到了 asbio package。 我是这么做的:
data library(asbio)
attach(data)
panel.cor.res(x, y, digits = 2, meth = "spearman", cex.cor=1)
panel.lm(x, y, col =1, bg = NA, pch = 16, cex = 1,
col.line = 2, lty = 1)
pairs(data, cex.labels=1, cex=.95, gap=.1, lower.panel=panel.cor.res, upper.
panel=panel.lm)
panel.cor.res 和 panel.lm 是asbio 里的function,可是每次总是出错:
> panel.cor.res(x, y, digits = 2, meth = "spearman", cex.cor=1)
Error in is.data.frame(y) : object 'y' not found
> panel.lm(x, y, col =1, bg = NA, pch = 16, cex = 1,
+ col.line = 2, lty = 1)
Error in points(x, y, pch = pch, col = col, bg = bg, cex = cex) :
object 'x' not found
实在不知道怎么解决了。大家能帮帮我么。
avatar
c*d
5
国内太push了。。。
不用care这个 国内能坚持下来的没多少 虽然现在吵得热热闹闹的 贵在坚持贵在坚持
不过 弹曲子 没问问是弹得什么?除了音阶 练习曲也算曲子 难道你家娃还在弹单个音
符?

住在SAN GABRIEL, CA. 想给大宝小宝找一个好的钢琴老师。 现在的钢琴老师太nice
, 孩子没啥
压力, 就是离家近。 想找个有经验, 严格点的。
听说很多老师每年都会安排学生参加一些表演, 这样家长也能了解, 老师和孩子都有
目标。 我们家俩
宝快5岁了, 跟现在的老师学了半年了, 能看懂简单的五线谱, 但手型和节奏很不规
矩。
听国内一教钢琴的人说, 4岁孩子3个月应该能弹曲子了, 我们的好像有点慢呀!
大家都有什么经验呀?
如有好的老师, 请推荐!
先谢谢了!

【在 d**********8 的大作中提到】
: 住在SAN GABRIEL, CA. 想给大宝小宝找一个好的钢琴老师。 现在的钢琴老师太nice
: , 孩子没啥
: 压力, 就是离家近。 想找个有经验, 严格点的。
: 听说很多老师每年都会安排学生参加一些表演, 这样家长也能了解, 老师和孩子都有
: 目标。 我们家俩
: 宝快5岁了, 跟现在的老师学了半年了, 能看懂简单的五线谱, 但手型和节奏很不规
: 矩。
: 听国内一教钢琴的人说, 4岁孩子3个月应该能弹曲子了, 我们的好像有点慢呀!
: 大家都有什么经验呀?
: 如有好的老师, 请推荐!

avatar
M*e
6
都是他说的?说得很好啊!

【在 c*******7 的大作中提到】
: 【 以下文字转载自 Joke 讨论区 】
: 发信人: linvk (aglee), 信区: Joke
: 标 题: 莫言最伤人的83句话,非常经典
: 发信站: BBS 未名空间站 (Fri Oct 12 01:38:46 2012, 美东)
: 看了不准哭

avatar
W*o
8
在微博上看到了,感觉有些不是他写的吧。
avatar
J*r
9
scatterplot也需要用package么。。。
anyway, 我觉得你的data input有点问题。
如果可以,你可不可以贴投一行of your data file上来
avatar
e*a
10
只有最后一句是他说的。

【在 M**e 的大作中提到】
: 都是他说的?说得很好啊!
avatar
x*m
11
1.you may want to use scatterplot() in R,which is easier to be handle for
new users.
2. Don't use "data" as variable name!!! It's a function in R!!! This mistake
is common in rookie programmers.
3. R complains it cannot find y variable. something wrong with your data or
variable name. post your data here.
avatar
s*a
12
“文图自动生成,不代表官方观点”

【在 e*a 的大作中提到】
: 只有最后一句是他说的。
avatar
h*t
13
谢谢大家的回复。我贴几行data 上来。我目的是想要做一个matrix 的scatter plot。
upper panel包含 correlation coefficient,lower panel 包含数据的scatterplot,
最好含有estimated trendline,线性即可,不用panel。smooth,那个smooth出来对我
的结果没有什么意义。
网上有一个写好的R程序,不需要用package,但是有很多function,我看不太懂,而且
不能换 correlation的method,比如pearson, spearman (或者是我换的地方不对,
出错了)。
以下几行是我的data,第一行是header,我用了header = TRUE,一共是8个column。我
希望能把这8个vector做一个 scatterplot matrix。请大家多给我指点指点。早上看到
bronstein的回复,我又试了试用read。table,还是没能解决问题。
num_rev_z size_tran_z mode_rev_z cir_cb_z t_log2_-z soft_sq_z
hard_sq_z sp_sqrt_z
-0.547192766 -0.868636422 -1.052436002 -0.484344371 0 1.
816662392 -0.373535719 0.095527427
-0.440678884 0.037311032 1.106336502 -0.711983195 0 -0.
284128179 -0.080620804 -0.921497753
avatar
h*g
14
这该多NC 才会相信这是莫言写的

【在 c*******7 的大作中提到】
: 【 以下文字转载自 Joke 讨论区 】
: 发信人: linvk (aglee), 信区: Joke
: 标 题: 莫言最伤人的83句话,非常经典
: 发信站: BBS 未名空间站 (Fri Oct 12 01:38:46 2012, 美东)
: 看了不准哭

avatar
c*g
15
你把这俩行去了试试
##
panel.cor.res(x, y, digits = 2, meth = "spearman", cex.cor=1)
panel.lm(x, y, col =1, bg = NA, pch = 16, cex = 1,
col.line = 2, lty = 1)
##
这俩行不需要的
另外
head(data)
看看你的数据是不是被读进来了

upper.

【在 h******t 的大作中提到】
: 我刚刚开始自学R,作图的时候遇到了一些问题。可能很基本,但是搜索了半天也没有找
: 到答案。这里高手多,请大家帮我看看。
: 我需要做一个 scatterplot 的matrix,找到了 asbio package。 我是这么做的:
: data : library(asbio)
: attach(data)
: panel.cor.res(x, y, digits = 2, meth = "spearman", cex.cor=1)
: panel.lm(x, y, col =1, bg = NA, pch = 16, cex = 1,
: col.line = 2, lty = 1)
: pairs(data, cex.labels=1, cex=.95, gap=.1, lower.panel=panel.cor.res, upper.

avatar
M*d
16
hao

【在 M**e 的大作中提到】
: 都是他说的?说得很好啊!
avatar
h*t
17
我用了head,读出了前6行,是对的么。如果用read之后直接>mydata,是可以把我的
data全部读出的。
不用那两行的话,如何能在upper panel 和 lower panel 输出不同的信息呢。比如
upper 输出数字的correlation,lower输出scattered plot,外带线性的 trendline?
我知道我的问题比较多,多谢多谢。

【在 c*****g 的大作中提到】
: 你把这俩行去了试试
: ##
: panel.cor.res(x, y, digits = 2, meth = "spearman", cex.cor=1)
: panel.lm(x, y, col =1, bg = NA, pch = 16, cex = 1,
: col.line = 2, lty = 1)
: ##
: 这俩行不需要的
: 另外
: head(data)
: 看看你的数据是不是被读进来了

avatar
c*g
18
那么你试了去掉那俩行的结果了吗?你的pairs里面的upper和lower的arguments指明了
用什么function去算了。

【在 h******t 的大作中提到】
: 我用了head,读出了前6行,是对的么。如果用read之后直接>mydata,是可以把我的
: data全部读出的。
: 不用那两行的话,如何能在upper panel 和 lower panel 输出不同的信息呢。比如
: upper 输出数字的correlation,lower输出scattered plot,外带线性的 trendline?
: 我知道我的问题比较多,多谢多谢。

avatar
h*t
19
啊,我太弱了,居然不知道可以直接用这两个function。我试了,可以作图。但是又有
了问题,只用最后一行的话,怎么样选择panel.cor.res 里面的meth呢?默认的是
pearson 吧,如果换成 spearman,我想当然的用了:
pairs(mydata, cex.labels=1, cex=.95, gap=.1, lower.panel=panel.cor.res(
method=”spearman”), upper.panel=panel.lm)
结果行不通。 :(((((

【在 c*****g 的大作中提到】
: 那么你试了去掉那俩行的结果了吗?你的pairs里面的upper和lower的arguments指明了
: 用什么function去算了。

avatar
c*g
20
我猜的,你可以试试,或者你可以尝试修改panel.cor.res
lower.panel = function(x, y, ...) { panel.cor.res(x, y, meth = "spearman") }

【在 h******t 的大作中提到】
: 啊,我太弱了,居然不知道可以直接用这两个function。我试了,可以作图。但是又有
: 了问题,只用最后一行的话,怎么样选择panel.cor.res 里面的meth呢?默认的是
: pearson 吧,如果换成 spearman,我想当然的用了:
: pairs(mydata, cex.labels=1, cex=.95, gap=.1, lower.panel=panel.cor.res(
: method=”spearman”), upper.panel=panel.lm)
: 结果行不通。 :(((((

avatar
h*t
21
成了!太感谢了!有高手点拨好过自己瞎撞啊。
我刚刚才接触 R,觉得很有用。可是没有什么编程的基础,每次都是靠google 或者
youtube,很多command看得不太懂是什么意思。能请你介绍几个适合新手学习的网站或
者资料么,谢谢了。

}

【在 c*****g 的大作中提到】
: 我猜的,你可以试试,或者你可以尝试修改panel.cor.res
: lower.panel = function(x, y, ...) { panel.cor.res(x, y, meth = "spearman") }

avatar
c*g
22
也许有的人会建议你看书,我个人的看法,对于R这种语法不严格的工具,由于在各种
mailing list上问答的太多了,其实 google + 实际问题 + 实践 是最好的学习方法。
当然如果你想看书,google一下R的tutorial,到处都是。
复杂点的问题,google + mailing list 提问会快很多。

【在 h******t 的大作中提到】
: 成了!太感谢了!有高手点拨好过自己瞎撞啊。
: 我刚刚才接触 R,觉得很有用。可是没有什么编程的基础,每次都是靠google 或者
: youtube,很多command看得不太懂是什么意思。能请你介绍几个适合新手学习的网站或
: 者资料么,谢谢了。
:
: }

avatar
h*t
23
非常感谢,看来不是一时之功。要多看多练。
我刚才用你给我改正的 function,也照样改了一下原来的,也可以画一样的图了 :)
,不过就是在 pair 里面不知道为什么不能用cex=。95了,如果删掉就可以。那两个星
号是我后加的,指出cex 的位置。可不可以再麻烦你指教一下。
> pairs.mine )}
> pairs.up cex = 1,
+ col.line = 2, lty = 1 )}
> pairs(zscore, cex.labels=1, gap=0.1, **cex=.95**, lower.panel=pairs.up,
upper.panel=pairs.mine)

法。

【在 c*****g 的大作中提到】
: 也许有的人会建议你看书,我个人的看法,对于R这种语法不严格的工具,由于在各种
: mailing list上问答的太多了,其实 google + 实际问题 + 实践 是最好的学习方法。
: 当然如果你想看书,google一下R的tutorial,到处都是。
: 复杂点的问题,google + mailing list 提问会快很多。

avatar
c*g
24
你的pairs.up里已经有cex了,把它删去就可以了。

pearson"

【在 h******t 的大作中提到】
: 非常感谢,看来不是一时之功。要多看多练。
: 我刚才用你给我改正的 function,也照样改了一下原来的,也可以画一样的图了 :)
: ,不过就是在 pair 里面不知道为什么不能用cex=。95了,如果删掉就可以。那两个星
: 号是我后加的,指出cex 的位置。可不可以再麻烦你指教一下。
: > pairs.mine : )}
: > pairs.up : cex = 1,
: + col.line = 2, lty = 1 )}
: > pairs(zscore, cex.labels=1, gap=0.1, **cex=.95**, lower.panel=pairs.up,

avatar
h*t
25
哦,你这样一说我知道了,cex.label已经定义了 vector 的label大小了,所以cex 本
来是定义pair.lm的,但是在pair.up里定义了,就不需要了,对吧。

【在 c*****g 的大作中提到】
: 你的pairs.up里已经有cex了,把它删去就可以了。
:
: pearson"

avatar
h*t
26
非常感谢你今天的耐心帮助,也谢谢楼上的各位。

【在 c*****g 的大作中提到】
: 你的pairs.up里已经有cex了,把它删去就可以了。
:
: pearson"

相关阅读
logo
联系我们隐私协议©2024 redian.news
Redian新闻
Redian.news刊载任何文章,不代表同意其说法或描述,仅为提供更多信息,也不构成任何建议。文章信息的合法性及真实性由其作者负责,与Redian.news及其运营公司无关。欢迎投稿,如发现稿件侵权,或作者不愿在本网发表文章,请版权拥有者通知本网处理。