Home Comments Thread
New Thread

408 Comments

Page 1 / 5
giscus-bot giscus-bot 2022-12-16 05:41:35
访客 *五脚星* @ 2007-08-11 20:48:07 写道:

😎

giscus-bot giscus-bot 2022-12-16 05:41:36
访客 *deadwind* @ 2007-09-28 10:12:13 写道:

师兄  那个统计之都得文档什么的是不是还不能下?什么时候能下呢?

yihui yihui 2022-12-16 06:34:14

这个……取决于Zhang, X. M.同学的速度……你可以问她……现在已经传了少数一部分上去了

——原帖发布于 2007-09-28 10:31:34

giscus-bot giscus-bot 2022-12-16 05:41:37
访客 *pisces0118* @ 2007-10-13 06:56:16 写道:

注册以后要到注册的邮箱去激活帐号怎么用啊

yihui yihui 2022-12-16 06:34:16

注册的时候填写真实可用的邮箱,随后会收到一封激活邮件,点击链接激活即可。

——原帖发布于 2007-10-13 12:56:25

giscus-bot giscus-bot 2022-12-16 05:41:37
访客 *lanyahbh* @ 2007-10-13 15:10:57 写道:

谢谢您在S-Plus/R中的指导,我想把我的问题再具体的描述一下,我的论文做的是:应用遥感数据及少量的地面调查数据对森林的蓄积量进行预测(预测变量为遥感影像的灰度值和地物类型,目标变量为森林蓄积量;它们都是以二维表的形式组织的),现在我基本上已经完成了应用样本数据拟合好了模型,并在知道了如何应用拟合好的对未知蓄积量的数据(也是二维表的形式)进行预测,现在剩下的问题是:如何把预测完了数据(二维表)转换成图像?
希望您可以抽空帮我看一下,像一个办法,在这里先谢谢您了!

yihui yihui 2022-12-16 06:34:15

我已经回复邮件。

参考image(),例:

x=y=seq(-6,6,.1)
z=outer(x,y,function(x,y)x^2+6*sin(y))
image(x,y,z)
box()

——原帖发布于 2007-10-13 08:34:33

giscus-bot giscus-bot 2022-12-16 05:41:38
访客 *lanyahbh* @ 2007-10-15 17:13:52 写道:

现在又有一个问题了,是这样的:我要用遥感信息(连续型)和地物类(字符型)型建立一个关于预测森林蓄积量(需要是连续型)的神经网络模型。可是在建立神经网络模型时,设置蓄积量为连续型变量,系统提示出错,说什么该变量不能为该类型!但在我的预测中,目标变量需要是连续型的(因为各块林地上的森林蓄积量是连续变化),而且当我把它改成是离散型(只是试验,实践是不应该这样的),预测的精度的各项结果均为0,非常的迷惑!!
抽个时间,帮忙看一下吧!
非常感谢

yihui yihui 2022-12-16 06:34:17

神经网络毫无疑问是可以用于连续型因变量的,你之所以出现这样的问题我猜还是对数据结构不熟悉,检查一下你的参数是否符合函数nnet()的要求,至少我用R的nnet包作连续变量的预测是没问题的。你用S-Plus还是R?(我没有S-Plus,但S-Plus的神经网络与R的函数是相同的名字,你可以参考Venables和Ripley他们那本Modern Applied Statistics with S-Plus)

> x=rnorm(20);y=x+rnorm(20,0,.2)
> m=nnet(y~x,size=2) #隐藏层中单元数为2
# weights:  7
initial  value 21.014472
iter  10 value 17.272225
iter  20 value 11.576956
iter  30 value 10.625425
iter  40 value 10.572321
iter  50 value 10.563150
iter  60 value 10.540803
iter  70 value 10.535614
iter  80 value 10.531502
iter  90 value 10.528969
iter 100 value 10.524533
final  value 10.524533
stopped after 100 iterations
> predict(m) #预测
          [,1]
1  7.192212e-06
2  6.848704e-06
3  9.273907e-06
4  9.999606e-01
5  6.860589e-06
6  5.305207e-01
7  4.216611e-01
8  9.785739e-05
9  2.618520e-01
10 6.847917e-06
11 4.857670e-05
12 6.848141e-06
13 6.847917e-06
14 5.796481e-04
15 6.847917e-06
16 9.999690e-01
17 7.793182e-01
18 9.999652e-01
19 6.847958e-06
20 1.269201e-03

——原帖发布于 2007-10-15 11:28:53

giscus-bot giscus-bot 2022-12-16 05:41:39
访客 *lanyahbh* @ 2007-10-15 22:56:37 写道:
X  Y  NDVI  W_PC_3  PC_1  HIS_3  Broad_laef  conifer  C_B  Bamboo  Nursery  shrub  Volume
20674  3553  -0.002  87.518  91.275  85.117  1  0  0  0  0  0  8.0111
20674.5  3553  -0.027  89.887  103.271  103.893  1  0  0  0  0  0  14.1876
20675  3553  -0.038  92.111  109.917  111.289  1  0  0  0  0  0  5.7136
20673  3552  0.002  87.145  99.851  88.566  1  0  0  0  0  0  6.8663
20673.5  3552  -0.018  82.301  79.805  69.661  1  0  0  0  0  0  10.7293
20674  3552  -0.022  89.515  104.273  99.619  1  0  0  0  0  0  3.9644
20674.5  3552  -0.021  85.004  107.69  92.174  1  0  0  0  0  0  5.127
20675  3552  -0.01  82.276  84.833  72.259  0  1  0  0  0  0  10.4952
20675.5  3552  -0.002  85.815  89.206  84.912  1  0  0  0  0  0  8.0866
20676  3552  -0.003  88.992  122.024  112.106  1  0  0  0  0  0  19.8122
20676.5  3552  -0.016  91.514  112.237  113.429  1  0  0  0  0  0  5.9058
20671.5  3551  -0.014  94.721  132.585  134.009    0  0  1  0  0  2.6946
20672  3551  0.001  86.335  80.149  71.95  1  0  0  0  0  0  4.6919
20672.5  3551  -0.005  84.086  79.04  73.22  1  0  0  0  0  0  8.1268
20673  3551  -0.035  90.996  121.337  125.117  0  0  1  0  0  0  6.3197
20673.5  3551  -0.051  92.234  126.423  120.782  1  0  0  0  0  0  7.9688
20674  3551  -0.026  88.714  119.846  112.791  1  0  0  0  0  0  3.7742
20674.5  3551  -0.05  91.436  115.112  111.937  1  0  0  0  0  0  7.165
20675  3551  -0.038  89.346  115.338  113.41  1  0  0  0  0  0  7.1586
20675.5  3551  -0.052  92.82  102.227  100.321  0  0  1  0  0  0  8.7565
20676  3551  -0.061  91.546  111.542  106.601  0  1  0  0  0  0  4.7667
20676.5  3551  -0.039  89.133  114.427  112.057  0  1  0  0  0  0  5.6961
20677  3551  -0.059  90.797  96.911  96.716  0  0  0  0  0  0  8.1386
20672  3550  -0.023  95.657  132.516  133.636  1  0  1  0  0  0  6.2305
20672.5  3550  -0.05  96.943  113.145  118.109  1  0  0  0  0  0  12.8494
20673.5  3550  -0.03  92.233  129.003  125.508  1  0  0  0  0  0  5.0575
20674.5  3550  -0.045  92.439  121.313  121.931  1  0  0  0  0  0  8.5315
20675  3550  -0.065  93.504  123.078  125.119  1  0  0  0  0  0  18.1014
20675.5  3550  -0.033  93.274  124.418  125.409  1  0  0  0  0  0  11.6885
20676  3550  -0.034  94.531  127.386  129.057  1  0  0  0  0  0  16.0616
20676.5  3550  -0.038  91.399  116.638  117.784  0  0  1  0  0  0  5.8788
20677  3550  -0.047  91.651  123.462  124.469  0  1  0  0  0  0  7.5143
20677.5  3550  -0.036  92.854  116.933  117.664  1  0  0  0  0  0  10.7145
20674  3549  -0.038  94.86  126.249  128.753  1  0  0  0  0  0  13.253
20674.5  3549  -0.045  96.565  123.839  127.273  1  0  0  0  0  0  9.0357
20675  3549  -0.039  97.486  113.069  115.813  0  0  0  0  1  0  1.4214
20675.5  3549  -0.063  96.6  120.225  126.182  1  0  0  0  0  0  18.8498
2674.5  3548  -0.059  98.558  122.468  124.699  1  0  0  0  0  0  12.244

预测变量是NDVI W_PC_3 PC_1 HIS_3 Broad_laef conifer C_B Bamboo Nursery shrub
目标变量是Volume

yihui yihui 2022-12-16 06:34:18

第12行是否含有缺失值?

——原帖发布于 2007-10-20 16:03:08

giscus-bot giscus-bot 2022-12-16 05:41:39
访客 *lanyahbh* @ 2007-10-16 06:45:50 写道:

谢谢你
我还是不懂,我用的是insightful
我下次把数据给你看一下

giscus-bot giscus-bot 2022-12-16 05:41:40
访客 *lanyahbh* @ 2007-10-26 07:56:38 写道:

谢谢!
这几天忙于其他事,都没有时间看您的回帖,我自己也研究了一下,基本上已经解决了

giscus-bot giscus-bot 2022-12-16 05:41:41
访客 *lanyahbh* @ 2007-10-26 08:02:08 写道:

不过有一点就是我的模型不是太稳定,就是说的参数的设置的不同,对起进行评价时,Relative S..  erorr相差很大

yihui yihui 2022-12-16 06:34:20

调整参数,找预测误差最小的

——原帖发布于 2007-10-29 16:46:32

giscus-bot giscus-bot 2022-12-16 05:41:41
Guest *anita_jiu* @ 2007-10-31 15:39:04 originally posted:

Hi, there
Have been extremely busy with the new challenge. Trying to get some sense out of a 'sophisticated' dissertation -  1.8 gender agreed...; 1.54 gender has heard of ... How 'creative'!! Applying mean, median, and mode to analyse norminal data!! ❓

Are u busy over there?

Take care.

A.

yihui yihui 2022-12-16 06:34:21

OMG... It's hard to imagine they can calculate ("cats"+"dogs")/2...
BTW, mode is appropriate for nominal data, but mean & median are absolutely NOT :)

Originally posted on 2007-10-31 18:49:20

giscus-bot giscus-bot 2022-12-16 05:41:42
访客 *xmg* @ 2007-11-10 21:49:13 写道:

看你的blog很有收获,收藏了。

yihui yihui 2022-12-16 06:34:22

不过我倾向于把有用的东西写在英文Blog中,呵呵

——原帖发布于 2007-11-10 22:43:11

giscus-bot giscus-bot 2022-12-16 05:41:43
访客 *等梦拜访* @ 2007-11-13 16:13:17 写道:

谢大人好~
最近偶在看对应分析
R中的MASS包偶怎么找了半天也没有找到啊
rcran上没有发现这个包:(
据吴老的书说corresp函数在这个MASS包中。。

yihui yihui 2022-12-16 06:34:23

CRAN上没有这个包?这玩笑开大了,MASS可是两位超级大牛Venables和Ripley的大作,它在R的众Packages中的Priority是“recommended”,不可能没有。你直接用install.packages()装就可以了,没必要上CRAN网站,你在网站上不能直接找到它是因为V和R两位大人把MASS和class、nnet、spatial四个包一起写成一个Bundle了,你在网站上肯定可以看见VR。

——原帖发布于 2007-11-13 17:53:56

giscus-bot giscus-bot 2022-12-16 05:41:43
访客 *等梦拜访* @ 2007-11-14 05:45:41 写道:

啊~搞定了
谢谢~~~~
原来我都是到cran网上找package。。
原来这样就可以啊

yihui yihui 2022-12-16 06:34:24

如果是用Windows的话,应该自然而然就用菜单安装了。主要是包与包之间有dependency,你若是搞不清楚,只从CRAN上面下载某一个zip文件安装,那么你就得挨个找相关联的包去安装,这是很麻烦的。RGui里面菜单默认就是把关联的包全都安装上,这也对应着install.packages()里面的一个参数,自己找。

——原帖发布于 2007-11-14 10:47:22

giscus-bot giscus-bot 2022-12-16 05:41:44
访客 *一支梅* @ 2007-11-24 05:50:25 写道:

你好!请问有最新版的马克威软件吗?如果有的话,能不能麻烦发到邮箱yizhimeihua2006@126.com

yihui yihui 2022-12-16 06:34:25

一支梅大人,貌似您还不了解我的立场,本人是坚决反对盗版统计软件的。

——原帖发布于 2007-11-24 09:10:09

giscus-bot giscus-bot 2022-12-16 06:36:05
访客 *ryl* @ 2009-07-03 13:04:43 写道:

呵呵,看到版主在FAQ里提的标准留言了

giscus-bot giscus-bot 2022-12-16 05:41:45
访客 *一支梅* @ 2007-11-24 05:51:01 写道:

补充下,谢谢!

giscus-bot giscus-bot 2022-12-16 05:41:45
访客 *dairyman* @ 2007-11-24 08:14:49 写道:

谢先生能否将A HANDBOOK OF STATISTICAL ANALYSES USING SAS给我一份,网上无法下载。多谢!

yihui yihui 2022-12-16 06:34:30

论坛里可以下载了。

——原帖发布于 2007-12-19 13:38:20

giscus-bot giscus-bot 2022-12-16 05:41:46
Guest *dairyman* @ 2007-11-24 08:23:00 originally posted:

email:xwbdairyman@yahoo.com.cn

giscus-bot giscus-bot 2022-12-16 05:41:47
访客 *lanyahbh* @ 2007-12-04 08:56:24 写道:

您好!
问一下:在SPSS里可不可以建立神经网络模型阿?
并且用建立的模型进行预测
在那个模块里?
我用的14.0的版本
谢谢!!!!

yihui yihui 2022-12-16 06:34:25

我基本不用SPSS。据我所知,SPSS本身能做的数据挖掘方法很少,他们把数据挖掘的模块都放在专门的产品Clementine里面做了。

——原帖发布于 2007-12-04 10:31:48

giscus-bot giscus-bot 2022-12-16 05:41:47
访客 *lanyahbh* @ 2007-12-04 11:23:09 写道:


这样啊

giscus-bot giscus-bot 2022-12-16 05:41:48
访客 *wb0800* @ 2007-12-06 22:14:02 写道:

板主,你好!
看到你的统计知识这么牛,本人好生佩服!希望今后能多多指教。
不知你是否有AMOS安装软件及教程,能否给我发过来,本人可能会用到这个软件,感激!
邮箱:wb0800@163.com
"0"是零,不是英文字母。
2007-12-7

yihui yihui 2022-12-16 06:34:26

你好,对不起我没有这个软件。我现在也比较反对SEM的建模方法。

——原帖发布于 2007-12-07 23:56:53

giscus-bot giscus-bot 2022-12-16 05:41:49
访客 *feng* @ 2007-12-13 10:02:41 写道:

你好,自己对统计学比较感觉兴趣,看到“统计之都”后比较喜欢想注册,但注册后多次无法收到,登陆时提示“你的帐号没有激活,请先到您注册的邮箱里激活帐号!”,烦请帮忙解决,多谢!用户名为:hefeng77或hefeng7710

yihui yihui 2022-12-16 06:34:27

两个帐号都已帮你激活。你的邮箱可能把激活邮件当作垃圾邮件了。

——原帖发布于 2007-12-14 09:09:56

giscus-bot giscus-bot 2022-12-16 05:41:50
访客 *zwdai* @ 2007-12-17 21:22:52 写道:

yihui你好,我也在用r,当然是非常初级的菜鸟。有几个概念的问题需要请教:
1)下列概念的区别在什么地方:
linear model,
nonlinearmodel
mixed linear model
general linear model
global linear moedel.
2)在何种情况运用上述不同的模型。
3)在r中,aov()和anov()的区别有哪些?
4)如果要做一个mixed linear model, 在选择fixed effect and random effect 以及covariate的时候有些什么好的建议,以及需要考虑什么问题。
希望指教!

zhwdai@yahoo.com

yihui yihui 2022-12-16 06:34:28

这问题真够多的,目前我只有能力回答一小部分。
1)线性模型不说了;非线性模型指的是参数是非线性形式;混合模型的系数有随机变量;一般线性模型和线性模型没什么区别,只是个称谓,但generalized linear model是另外一码事;global那个不懂。
2)问题太大。符合假设条件、效果好就用。
3)aov是普通意义上的方差分析;而anova一般是针对模型而言的比较;
4)我对混合模型并不熟悉,抱歉。

——原帖发布于 2007-12-19 13:34:24

giscus-bot giscus-bot 2022-12-16 05:41:50
访客 *millicent* @ 2007-12-18 20:44:39 写道:

您好!看到您统计研究得如此透彻,很是佩服哦,我是英语专业的硕士,对统计的要求没有那么高深,但是如果统计学得不好,软件用的不确切也许在圈内可以混过一些人,但是最终也会随着时间被揭穿的,而且真正要做好研究还是把该学的学到家。问题是我的导师也不太懂统计(别见笑,呵呵),问别人没问着很懂得呢,我自己钻研了好久SPSS,好不容易感觉搞懂了,可是你这里说SPSS不是很好,我需要建一个结构方程模型,买了些书,也没看明白怎么编程,也不知哪儿能弄到prelis,好不容易知道amos可以建模,也可以下到软件,您说您不太支持这个建模。请您指点一下,对于英语教育研究中变量不很多的统计我该使用什么样的软件分析变量,什么建模的软件比较呢?而且不要笑话,太贵的软件学生买不起,我们做研究没有一分钱经费。收到回帖将不胜感激!邮箱:millicent1@163.com

yihui yihui 2022-12-16 06:34:29

现有问题再有数据最后有模型,这才是做学问的思路。免费统计软件有R啊,比商业软件好多了,只是怕你们用不来,因为它实在是不傻瓜,不懂统计的话几乎没法用。
另:我研究得其实很不透彻……

——原帖发布于 2007-12-19 13:37:52

giscus-bot giscus-bot 2022-12-16 05:41:51
访客 *zwdai* @ 2007-12-19 15:40:07 写道:

呵呵,多谢yihui的回复。
不过看了回复还是感到有些替你担心,因为你说到自己要到什么地方授课!
如果没有搞清楚模型之间的区别,无论什么样的数据到手都一概采用可加线性模型去套,而不去察看是否满足模型的统计学前提,似乎不太可取。
因为前些日子我便吃过这种亏,才努力想去弄清楚它们之间的区别。
前些日子做的统计分析,用了一般线性模型,结果最后发现误差不是正态分布,被老板打回原型,回炉重做。
至于aov及anov之间的区别,我还知道的一点便是,aov必须是balanced的样本,anov则没有这个要求。

不过多些楼主的回复,希望把这里做的更好。以后还会请教。

yihui yihui 2022-12-16 06:34:30

嗯,谢谢你的提醒。其实我自己还是挺注意这些的,我要是自己不明白就说不明白,没搞明白之前我不会出山乱讲的:)
P.S. 你那个应该是“残差”,误差只是理论上存在的东西

——原帖发布于 2007-12-19 16:13:08

giscus-bot giscus-bot 2022-12-16 05:41:52
访客 *zwdai* @ 2007-12-19 15:55:48 写道:

察看了楼主其他部分的内容,多谢分享!

看老外吃粽子,确实让人捧腹,haha

yihui yihui 2022-12-16 06:34:34

你要是看到什么好玩的笑话也可以帮我贡献一些,我平时在网上看Joke,看到好玩的就收集到一个文件中,每次随机读取一条放在这个页面中:)

——原帖发布于 2007-12-20 12:59:30

giscus-bot giscus-bot 2022-12-16 05:41:52
访客 *millicent* @ 2007-12-20 12:35:26 写道:

非常感谢您的回复!我是先有问题再有数据再有模型的哦,我用多元回归分析做的变量之间相关性很大的,我觉得这是一个很有意思的发现,而且很有意义,以前很少有人做,做的也是因为没读懂原著,做得不好。可是现在我就是不会比较专业的建模方法,我觉得还是要靠自己努力啊,周围的人懂得很少,我在北京上硕士哦,大家都不会,也觉得不需要,我倒是对这个很感兴趣,做好研究不会统计还是不行的,压根可能连人家的论文都看不动,我去顶尖级的研究中心学过一阵,老师也不敢恭维,我一定要学好它,别人会我能会,非常感谢回帖,我再看看这个R是什么。不过我没见到我的专业国内外的文献有用这个的,呵呵。看了再决定吧,再次谢谢!

yihui yihui 2022-12-16 06:34:32

嗯,不错不错。结构方程模型现在过于热火朝天,而我很久以来就想写一篇批判文章,我们往往只顾别人模型的所谓“优点”,却忽视了很多严重的问题,典型的就是计算。更有甚者,竟然将结构方程模型称之为统计模型之母,简直是令人愤慨。

——原帖发布于 2007-12-20 12:58:17

giscus-bot giscus-bot 2022-12-16 05:41:53
访客 *millicent* @ 2007-12-24 18:55:42 写道:

又留言了,真是不好意思,我最近被这个结构方程模型快搞晕了,amos学生版只能处理8个变量,我刚好有9个,lisrel我没有prelis啊,弄不到,不知去哪弄啊,你说的r,我去官方网站看了,为什么没有中国的镜像站?试了几个链接没有windows的,再者我现在再去摸索这个不知是不是很浪费时间?还有我们中国外语教育研究中心也不用这个软件(当然这不代表它不好,但一定程度上说明如果我用了我还要和那些人解释它的好处,那些experts是否信任我还不好说),我在准备考博,总感觉时间很紧张,虽然没打算一定考上,但是我还是想赶紧弄完论文,看点书。如果你有的话,能否将amos正式版和prelis传给我,不胜感激!我还有一点不明白就是,R既然是免费,又很好为什么用的人很少呢?是太高深了呢?还是什么别的原因呢?既然很好为什么不做得好用些呢?我不懂,完全外行的问题,别见笑!如能把软件传给我,不胜感激!邮箱:millicent1@163.com,QQ:343577096

yihui yihui 2022-12-16 06:34:35

国内用的人相对少一些,毕竟是英文软件。国外统计界的很多大牛们都用R或S自己写统计程序。这种语言对于统计来说已经是非常方便了,历史上统计软件中只有S获得过ACM软件系统奖,没见那些“好用的”SPSS或者“强大的”SAS获此殊荣,还是有一定道理的。R是基于S语言开发起来的免费软件。中国镜像本来有,但后来那个学校不干了,因为设立镜像不会有人给你钱。
SEM其实并不是特别适合在R里面做,原因之一是R的sem包的功能还相对比较单一,之二是你要是搞清楚了R是怎么做的你就明白了SEM的不合理之处,为了避免把SEM的理论搞清楚,最好还是不要用开源软件,大家都蒙在鼓里用商业软件算好了。我不是卖软件的,无法提供商业软件。

——原帖发布于 2007-12-24 21:09:48

giscus-bot giscus-bot 2022-12-16 05:41:54
访客 *millicent* @ 2007-12-25 07:56:11 写道:

感谢你的回复!虽然没怎么看懂,但是我懂了你和你的这个搏克!

yihui yihui 2022-12-16 06:34:35

我怎么听着这话心里直打鼓……就像是电视剧里面经常说的“我算是看穿你了……” 😱

——原帖发布于 2007-12-25 09:48:27

giscus-bot giscus-bot 2022-12-16 05:41:54
访客 *fairy* @ 2007-12-28 09:54:00 写道:

呵呵,终于找到你这里了~~~
嗯,以后有很多需要跟你请教的啊

yihui yihui 2022-12-16 06:34:36

你你你……懂装不懂……

——原帖发布于 2007-12-28 14:31:36

giscus-bot giscus-bot 2022-12-16 05:41:55
访客 *Nurseshark* @ 2007-12-31 20:11:43 写道:

我很久很久以后终于又找到你的BLOG了
以后敬仰的向你学习画图 😁

yihui yihui 2022-12-16 06:34:37

你应该知道我有实名制的习惯吧,无论Google或Baidu,搜“谢益辉” or "Yihui Xie" or "Yihui"都可以找到这里,不过一般人不会去搜别人名字的:)
我画图功夫可以说是差到极点……我只会用软件…… 😨

——原帖发布于 2007-12-31 20:42:37

giscus-bot giscus-bot 2022-12-16 05:41:56
访客 *fx2006ms* @ 2008-01-01 14:13:17 写道:

谢谢你的新年贺卡,呵呵!
有一天我发现你是YC人,真高兴
我在那儿上的大学,沾亲带故地说我们是HB老乡
祝新年快乐,学业有成!

yihui yihui 2022-12-16 06:34:38

哦?在YC上的大学?莫非是SX大学?呵呵

——原帖发布于 2008-01-01 18:00:59

giscus-bot giscus-bot 2022-12-16 05:41:56
访客 *fx2006ms* @ 2008-01-01 18:35:03 写道:

YC也就有一所SX大学吧,呵呵

yihui yihui 2022-12-16 06:34:39

貌似也是,除了SX大学,别的我也不知道了

——原帖发布于 2008-01-01 19:07:45

giscus-bot giscus-bot 2022-12-16 05:41:57
访客 *fairy* @ 2008-01-02 14:36:26 写道:

新年快乐,可爱的益辉大人同学 😁

yihui yihui 2022-12-16 06:34:39

看样子我不用再回祝你了,都已经这样笑哈哈了 🤪

——原帖发布于 2008-01-02 15:24:22

giscus-bot giscus-bot 2022-12-16 05:41:58
访客 *janey* @ 2008-01-22 06:34:38 写道:

楼主您好,最近因为要想用R来作植物的DNA序列进化树,访问了一些较好的网站,如:统计之都:http://www.cos.name ,统计学http://stat.ruc.edu.cn
发现你是这些网内的牛人。
想请教您两个问题
1。下载的DNA序列中硷基显示不正常的问题
通过R命令从Genbank中下载了一些DNA序列,按照有关的手册提示的步骤进行操作。

a<-read.GenBank("AY748885","AY833051")
str(a)
从返回的结果看,本来应该是atgc字符的序列,都变成了分别用88  18  28  28表示的序列。
说明书手册上的例子并没有出现这样的结果。真不知道是怎么回事?
2。用R从网上下载DNA序列,是否可以添加数据库
choosebank()
的结果中,给出了一堆的数据库。除了"genbank"是我要的外,我的有些序列必须在别的库中找参考序列,而这个库没有在choosebank()命令返回的结果中。是否有办法将不在该范围的数据库加到R中,以进行网上下载序列?

对了,上述的命令要调用四个软件包:

library(ape)
library(apTreeshape)
library(ade4)
library(seqinr)

不知您能否抽点时间指点一二。感谢!
再次感谢人大的同学和老师在相关论坛上提供的学习材料!

👍 👍

yihui yihui 2022-12-16 06:34:40

1、我对生物非常不熟悉,你提到的数据我也不知道是什么格式的,你可以把str(a)贴出来我看看;
2、看帮助文件吧,实在不行就问作者。

——原帖发布于 2008-01-22 09:34:59

giscus-bot giscus-bot 2022-12-16 05:41:58
访客 *janey* @ 2008-01-24 10:09:54 写道:

谢谢指点,我再看看帮助文件!不懂再请教您!

giscus-bot giscus-bot 2022-12-16 05:41:59
访客 *sarah79* @ 2008-01-25 19:39:34 写道:

不得不佩服谢老师的才气,我会再接再厉的向您学习! 👍

giscus-bot giscus-bot 2022-12-16 05:42:00
访客 *yil036* @ 2008-02-15 22:16:06 写道:

楼主你还有The R Book的电子版么?我搜遍网络也没找到。能不能传给我?我的email是yingli036@gmail.com

yihui yihui 2022-12-16 06:34:41

论坛上似乎发过吧,我最近没时间弄这些事。

——原帖发布于 2008-02-23 19:49:14

giscus-bot giscus-bot 2022-12-16 05:42:00
访客 *双江* @ 2008-02-28 05:52:21 写道:

很有水平,技术、内容以及坚持的恒心,堪称标杆!

yihui yihui 2022-12-16 06:34:42

过奖了!:)

——原帖发布于 2008-03-27 08:31:47

giscus-bot giscus-bot 2022-12-16 05:42:01
访客 *en* @ 2008-03-14 15:35:03 写道:

随机过程哪本教材比较好呢????谢谢!!

yihui yihui 2022-12-16 06:34:43

我对随机过程并不熟悉……

——原帖发布于 2008-03-27 08:31:59

giscus-bot giscus-bot 2022-12-16 05:42:02
访客 *christina_yai* @ 2008-04-05 07:07:29 写道:

楼主,我听我们的老师说统计专业在国外可是相当的吃香啊 ,我想如果本科毕业后去国外留学 ,你觉得哪个国家哪个学校的统计学最牛啊?谢谢

yihui yihui 2022-12-16 06:34:44

佛曰:不好说,不好说,一说便是错。建议你(1)知道自己的兴趣所在(2)早点在国外大学找老师,多多交流(3)趁早读一读英文教材,熟悉专业术语(4)如果有师兄师姐出国了,也可以联系他们征求意见

——原帖发布于 2008-04-07 17:10:43

giscus-bot giscus-bot 2022-12-16 05:44:03
访客 *flingzhang* @ 2008-04-08 16:28:22 写道:

老兄,你的新blog没有老版看着舒服!个人感觉!还有你的统计功夫真牛啊!有空好好讨教一下!

yihui yihui 2022-12-16 06:34:45

过奖了,没你想象得那么厉害:)

——原帖发布于 2008-04-11 14:21:52

giscus-bot giscus-bot 2022-12-16 05:44:03
访客 *Yuan* @ 2008-04-10 08:36:51 写道:

用R编译的聚类算法,请问输入的数据有什么样的要求呢?

yihui yihui 2022-12-16 06:34:44

参见相应算法。聚类的数据一般不都是矩阵么?(甚至所有统计分析的数据都是矩阵形式)

——原帖发布于 2008-04-11 14:21:20

giscus-bot giscus-bot 2022-12-16 05:44:04
访客 *老刘* @ 2008-04-16 17:07:26 写道:

如何对apgarch模型的参数进行估计,
谢谢

yihui yihui 2022-12-16 06:34:46

抱歉,不懂……

——原帖发布于 2008-04-20 10:40:47

giscus-bot giscus-bot 2022-12-16 05:44:05
访客 *老刘* @ 2008-04-21 05:06:16 写道:

Thanks all the same,谢大侠。

giscus-bot giscus-bot 2022-12-16 05:44:05
访客 *英波* @ 2008-05-12 20:04:29 写道:

你好,我是学计量经济学的新手,想学习RATS软件,不知道你那里有没有RATS软件,是否可以发给我.我的邮箱是
yinbo0001@sohu.com

yihui yihui 2022-12-16 06:34:47

我早已声明我和这些商业统计软件“划清界限”了,请不要再追问,谢谢。

——原帖发布于 2008-05-12 20:54:32

giscus-bot giscus-bot 2022-12-16 05:44:06
访客 *brenda* @ 2008-06-02 21:40:19 写道:

👍 大牛~

yihui yihui 2022-12-16 06:34:48

where where...

Originally posted on 2008-06-03 07:03:17

giscus-bot giscus-bot 2022-12-16 05:44:07
访客 *青衫* @ 2008-08-03 13:32:46 写道:

👍

giscus-bot giscus-bot 2022-12-16 05:44:07
访客 *青衫* @ 2008-08-03 13:33:09 写道:

谢益辉大哥,你好:小妹,给你发过邮件,不知道你有没有收到。我现在要做一个结构方程的报告,急需amos,能不能麻烦加我msn或qq,传个amos给我啊?拜托拜托,真的是急用。我这几天会一直在线等的。
给你添麻烦了,不好意思啊!
beixinjiaoji@hotmail.com

yihui yihui 2022-12-16 06:34:51

参见SPSS 15.0

——原帖发布于 2008-08-07 09:42:54

giscus-bot giscus-bot 2022-12-16 05:44:08
访客 *adelle* @ 2008-08-07 13:54:41 写道:

你好,在查找数据的标准化的时候浏览到你的主页
首先是深深感谢
学习过程中,要处理一些数据
很多不能理解
3q

yihui yihui 2022-12-16 06:34:52

标准化很容易的,有问题尽管留言:)

——原帖发布于 2008-08-07 14:31:02

giscus-bot giscus-bot 2022-12-16 05:44:09
访客 *zwdbordeaux* @ 2008-08-12 11:59:51 写道:

为什么最近几天我在这边没有办法打开COS论坛?是不是楼主离开之后,服务器发生了什么问题?
还是我这边的问题?

yihui yihui 2022-12-16 06:34:53

没办法,因为服务商在奥孕期间需要避运……

——原帖发布于 2008-08-20 18:22:47

giscus-bot giscus-bot 2022-12-16 05:44:09
访客 *zwdbordeaux* @ 2008-08-20 13:45:54 写道:

避到什么时候是个头?不要从此都没有办法看才好,上了统计之都,依然没有办法进一步点入!急躁。
今天早上发现小木虫也避运了,这事弄得!

yihui yihui 2022-12-16 06:34:54

等运期过了就好了,唉

——原帖发布于 2008-08-20 18:23:59

giscus-bot giscus-bot 2022-12-16 05:44:10
访客 *macross* @ 2008-09-15 05:18:54 写道:

你好,我用R语言绘出了一些统计图,但是PDF格式的,那么该图如何保存呢?又如何转化成其他图片格式呢?谢谢。

yihui yihui 2022-12-16 06:34:54

参见grDevices包,R有很多图形设备,PDF只是其中一种,还有Postscript、PNG、JPEG、TIFF、BMP等等。
?pdf
?png
?bitmap
……

——原帖发布于 2008-09-15 13:28:25

giscus-bot giscus-bot 2022-12-16 05:44:11
访客 *曼妮* @ 2008-09-29 10:42:15 写道:

google零假设,误打误撞不成想就进了楼主大人的blog了。楼主这里建设的那是上有形式下有内容、外有层次内有深度,那是果然……相当的深啊~~
四下一看,忽然就惭愧的感到我脖子上生了个外星脑袋了,要不除了笑话我怎么什么都不好理解呢?

p.s. 本人学医的,这个……职业差别,职业差别 😓 。
p.s.2 楼主很强,在我最菜的地方,强的像是装甲的,佩服,献花! 😳

yihui yihui 2022-12-16 06:34:55

过奖了,术业有专攻,没有高下之分 🐷

——原帖发布于 2008-09-30 04:06:41

giscus-bot giscus-bot 2022-12-16 05:44:11
访客 *麦克斯韦小妖* @ 2008-10-14 06:00:27 写道:

楼主,请问在R语言中如何在同一个坐标下画多条曲线,类似于MATLAB下的plot(x1,y1,x2,y2),并用不同的标志表示出?谢谢。

yihui yihui 2022-12-16 06:34:56

先plot(),再用lines()逐条线添加
或者用matplot(), matlines(), ...
详情看帮助。

——原帖发布于 2008-10-14 07:41:46

giscus-bot giscus-bot 2022-12-16 05:44:12
访客 *zhoulvjun* @ 2008-10-15 13:18:43 写道:

您门户首页上的“Fisher随便想出来的”所指的网址已经变动了,现在是“http://www.jerrydallal.com/LHSP/p05.htm”了...提醒一下,莫怪!

yihui yihui 2022-12-16 06:34:57

感谢还来不及呢,咋会怪呢:)

——原帖发布于 2008-10-15 15:48:36

giscus-bot giscus-bot 2022-12-16 05:44:13
访客 *lisle* @ 2008-10-22 18:12:40 写道:

谢老师,AMOS中,做路径分析,是否可以设置控制变量?
如何设置?

诚恳求教。

我的邮箱:soliloquist@163.com

yihui yihui 2022-12-16 06:34:58

很抱歉,我不懂……

——原帖发布于 2008-10-22 19:34:44

giscus-bot giscus-bot 2022-12-16 05:44:14
访客 *lora* @ 2008-10-23 13:44:46 写道:

您好,看到您有amos的软件,所以想从您那下载,您能不能把gmail的账号和密码给我?
我的信息在网上不便留,您方便的话发邮件给我,我在邮件里都会说明.多谢
我的email:lorashen@126.com

yihui yihui 2022-12-16 06:35:00

自己找SPSS,我不是卖软件的,谢谢

——原帖发布于 2008-10-24 09:12:58

giscus-bot giscus-bot 2022-12-16 05:44:14
访客 *ananjingjing* @ 2008-10-23 16:06:43 写道:

您好:
   我想参加这次的R会议,因为是新手,所以只是想听老师们讲。不知道在哪里报名。

yihui yihui 2022-12-16 06:34:59

会议通知第七条:http://stat.ruc.edu.cn/cn/notice/58785.html

——原帖发布于 2008-10-24 09:08:30

giscus-bot giscus-bot 2022-12-16 05:44:15
访客 *wf* @ 2008-11-01 14:01:34 写道:

我是新手,对贵论坛,也对R语言,不知有资格请教吗?因为学生物多样性,所以想学R语言,不知从何开始。谢谢指点

yihui yihui 2022-12-16 06:35:00

啥叫“资格”?问问题逮着人问就是了。打开R,输入1+1,得到2,你就入门了。然后浏览一下R-intro,熟悉数据结构之后再看统计部分,不断看帮助。所有事情的前提是,你最好有点英语底子以及足够的耐心。

——原帖发布于 2008-11-02 09:36:50

giscus-bot giscus-bot 2022-12-16 05:44:16
访客 *zhoulvjun* @ 2008-11-03 17:18:44 写道:

确实是按您说的那样去做的....

不管我怎么改swftools的赋值,都是类似的信息,无法正确解析目录名,这也是为什么我把png2swf.exe文件复制到c盘根目录的原因。

giscus-bot giscus-bot 2022-12-16 05:44:16
访客 *zhoulvjun* @ 2008-11-03 19:18:59 写道:

将函数修改如下,即去掉dQuote()就能显式地正常使用参数swftools来调用了:
   tool = ifelse(is.null(swftools), paste(dev, "2swf", sep = ""), file.path(paste(swftools, paste(dev, "2swf", sep = ""),sep="")))

如果不修改函数,只支持环境变量中指定路径的方式引用,因为问题的关键部分是在dQuote(file.path(swftools, paste(dev, "2swf", sep = ""))))这部分上,而在使用环境变量的情况下是不会执行此部分的。

yihui yihui 2022-12-16 06:35:04

谢谢,这是一个小Bug,应该用shQuote()而不是dQuote(),下一版本会修正过来。本来以前是没有引号的,后来发现有人把SWF Tools安装在带有空格的目录下,如果不加引号就无法正确执行命令。

——原帖发布于 2008-11-04 10:34:56

giscus-bot giscus-bot 2022-12-16 05:44:17
访客 *zhoulvjun* @ 2008-11-03 22:22:19 写道:

今天在R 2.8.0装了animation包,并把swftools中的png2swf.exe复制到c盘根目录,调用saveSWF()函数时指定swftools="c:/",但运行后总提示swftools not found; please install swftools first: http://www.swftools.org

后来看saveSWF.r中的代码,逐步调试发现是在函数体中执行version <- system(tool, intern = TRUE)后,version值为:"c:"; 不是内部或外部命令,也不是可运行的程序或批处理文件。而此时tool值为:"c://png2swf"

不管我怎么改swftools的赋值,都是类似的信息,无法正确解析目录名,这也是为什么我把png2swf.exe文件复制到c盘根目录的原因,即避免目录复杂。但如果我在console中直接赋值给tool="c://png2swf",并执行system(tool, intern = TRUE)能返回16行的值。

如果我把函数中的判断version的语句删掉,那么结果是Flash will be created at:
Executing: "c://png2swf" Rplot*.png -o movie.swf -r 0.666666666666667
Warning message:
In system(convert) : "c://png2swf" not found

我该怎么办呢?谢谢!

yihui yihui 2022-12-16 06:35:01

装好SWF Tools之后告诉saveSWF()的swftools参数你的软件安装路径即可,保证需要的程序在那个路径下能找到,没有必要把某个程序单独复制出来。看saveSWF()的帮助。

如果不想每次都给swftools参数赋值,那么就把SWF Tools的路径添加到系统的环境变量PATH中。

——原帖发布于 2008-11-03 22:31:49

giscus-bot giscus-bot 2022-12-16 05:44:18
访客 *在旅途* @ 2008-11-04 03:55:00 写道:

最近焦头烂额啃统计,几次搜到您的大名了,今天找偏最小二乘回归,没想到又跑到这里了。不敢探讨,实在是外行,请教或者救助可以吗?我有一批数据用了4因素饱和D-最优设计,测了几个方面的不同指标,啃来啃去,啃到偏最小二乘回归了,可是用什么软件来还是一头雾水,Spss简单一些的分析还行,SAS直接外行,其他统计软件更不会了。能请您百忙之中帮忙指点一下吗?不胜感激啊。
我QQ:553911749
  E-mail: sunxiaoqin74@163.com

yihui yihui 2022-12-16 06:35:03

抱歉,你的情况我可能无能为力。我一般只用R,估计你就更不会了。

——原帖发布于 2008-11-04 10:25:31

giscus-bot giscus-bot 2022-12-16 05:44:19
访客 *zhoulvjun* @ 2008-11-05 20:23:16 写道:

在看您的论文《基于R软件rpart包的分类与回归树应用》,但其中的数据http://www.stanford.edu/class/stats202/DATA/stagec.data目前已经无法下载,可否email我一份?我的邮箱:319336688@qq.com,谢谢!
另外,您在文中说“相应复杂性参量cp必须大于0.016 462(本例中选为0.017)”,请问这是通过怎样的计算过程得到的呢?
请指教,谢谢!

yihui yihui 2022-12-16 06:35:05

刚检查了一下,确实链接已经失效,我Google了一下,发现这里还有那批数据:http://www.ms.unimelb.edu.au/~odj/Teaching/dm/index.html(查找stagec.csv)
cp是根据1-SE规则算出范围再取的。

——原帖发布于 2008-11-06 09:23:27

giscus-bot giscus-bot 2022-12-16 05:44:20
访客 *zhoulvjun* @ 2008-11-06 02:58:37 写道:

谢谢您提供的数据!
我用数据试算的结果与您在文中print(fit)的结果完全一致,说明两套数据一致。但:
1、fit$cptable的结果中,最后两列不一致。为全面说明情况,我贴出 printcp(fit)的结果:

Survival regression tree:
rpart(formula = Surv(pgtime, pgstat) ~ age + eet + g2 + grade +
   gleason + ploidy, data = stagec, method = "exp")

Variables actually used in tree construction:
[1] age     g2      gleason grade  

Root node error: 195.41/146 = 1.3384

n= 146

       CP nsplit rel error  xerror     xstd
1 0.128256      0   1.00000 1.00593 0.073366
2 0.041437      1   0.87174 0.88799 0.076052
3 0.028913      2   0.83031 0.93030 0.083594
4 0.017707      3   0.80139 0.96093 0.086121
5 0.016462      4   0.78369 0.99496 0.088869
6 0.011390      5   0.76723 1.00928 0.088951
7 0.010000      6   0.75584 1.01495 0.091569
> 0.88799+0.076052
[1] 0.964042

按照1-SE规则,应是3、4行入选,此后依据最小cp原则,选取行4,即最大节点数为3。虽然结果与您的一致,但判断过程显然不同(如果按您原文中的数据,在4、5行中选择行4显然是选择了可选的最大cp,这也是我不明白“相应复杂性参量cp必须大于0.016 462(本例中选为0.017)”的原因,因为它与前面说的最小cp原则直接矛盾)。

2、我选择留言为“隐藏”,为何在页面的左侧“最新留言”处能见到我的加密留言呢?

我是初学者,不懂之处很多,实在是因为求知心切,有冒犯之处还要请您谅解!

yihui yihui 2022-12-16 06:35:06

1、误差是交叉验证计算出来的,而交叉验证的过程涉及到随机抽样,每次抽样结果不同,计算的误差也就不同,这是正常的。
2、显示留言的插件有Bug,这些留言不用加密吧?

——原帖发布于 2008-11-06 16:33:40

giscus-bot giscus-bot 2022-12-16 05:44:20
访客 *zhoulvjun* @ 2008-11-07 09:52:58 写道:

在plot.rpart绘的树中,第一层节点与第二层节点之间的距离实在是太夸张的长了,以致于最后一层节点的相应标签都显示不全,好像没有合适的参数可以用来调整。请您指点,谢谢!

yihui yihui 2022-12-16 06:35:06

这些距离是有特殊意义的。不要用“好像”,?plot.rpart,第二个参数便是

——原帖发布于 2008-11-07 11:32:53

giscus-bot giscus-bot 2022-12-16 05:44:21
访客 *zhoulvjun* @ 2008-11-07 12:33:43 写道:

再问一个rpart的问题,真是不好意思,老是麻烦您些小问题!

这批数据我运行出这样一个结果:

fit$cptable
         CP nsplit rel error    xerror       xstd
1 0.12825561   0 1.0000000 1.0159632 0.07328409
2 0.04143718  1 0.8717444 0.9137376 0.07548916
3 0.02891302  2 0.8303072 0.9448839 0.08067524
4 0.01770700  3 0.8013942 0.9886030 0.08655039
5 0.01646215   4 0.7836872 0.9718825 0.08479855
6 0.01138994  5 0.7672250 0.9741466 0.08478884
7 0.01000000  6 0.7558351 0.9902167 0.08427587

0.9137376+0.07548916
[1] 0.9892268  按照1-SE标准,第3至第6行的xerror都是小于这个值的,如果选择最小的cp,岂不是选择6个节点?但我觉得4个节点仍旧比较合适,但它的cp和xerror都比5、6个节点大啊!怎么判断才比较好呢!

yihui yihui 2022-12-16 06:35:08

为什么你觉得4个节点合适?

——原帖发布于 2008-11-09 19:01:47

giscus-bot giscus-bot 2022-12-16 05:44:21
访客 *关菁菁* @ 2008-11-09 12:03:54 写道:

突然就进来了,又让我开始思考

yihui yihui 2022-12-16 06:35:07

思考啥?:)

——原帖发布于 2008-11-09 18:51:34

giscus-bot giscus-bot 2022-12-16 05:44:22
访客 *zhoulvjun* @ 2008-11-09 20:22:13 写道:

呵呵,选4个节点是综合了其它几次运行结果后的感觉,也算是“民主投票”方式吧!此外,由于Variables actually used in tree construction:
[1] age     g2      gleason grade
,如果最后的树能够充分使用这四个变量的话,我想也是比较合理的,即5个节点。

giscus-bot giscus-bot 2022-12-16 05:44:23
访客 *wg* @ 2008-11-15 13:32:44 写道:

问一下使用R是否可以对jpg格式的数码照片进行处理,输入图片得到图片中一个个椭圆的轴长、面积等?谢谢!

yihui yihui 2022-12-16 06:35:09

rimage包中有读取JPEG格式图片的函数,读进来之后是一个三维数组,分别表示三原色,你要考虑的是怎样从颜色识别封闭椭圆的位置。

——原帖发布于 2008-11-16 17:12:02

giscus-bot giscus-bot 2022-12-16 05:44:24
访客 *随风* @ 2008-11-17 18:08:01 写道:

😳 很有意思,对统计学知之甚少,但未来一年却要频繁用到,看了这里,忽然对它感兴趣!

yihui yihui 2022-12-16 06:35:10

统计本来就是很有意思的学科 😁

——原帖发布于 2008-11-18 16:05:58

giscus-bot giscus-bot 2022-12-16 05:44:24
访客 *struggle12* @ 2008-11-25 11:47:43 写道:

非常好的论坛!期待有好的成果!

yihui yihui 2022-12-16 06:35:11

谢谢

——原帖发布于 2008-11-25 12:15:40

giscus-bot giscus-bot 2022-12-16 05:44:25
访客 *sh* @ 2008-11-26 06:37:21 写道:

我问问题的邮件没有回复~~是不是应该来这里问啊~~

yihui yihui 2022-12-16 06:35:12

你问了我也不知道你是谁啊……

——原帖发布于 2008-11-28 09:10:49

giscus-bot giscus-bot 2022-12-16 05:44:26
访客 *小鸡快跑* @ 2008-11-27 12:39:20 写道:

顶一个~
我同样学统计~同样对计算机感兴趣~偶尔做站~
要是早发现有这个论坛和站就好了。。。。

yihui yihui 2022-12-16 06:35:13

这位同仁是?一起加入统计之都吧 http://cos.name

——原帖发布于 2008-11-28 09:11:50

giscus-bot giscus-bot 2022-12-16 05:44:27
访客 *yanlinlin82* @ 2008-11-28 14:07:27 写道:

无意间发现这里,不错,支持!我是学生物的,希望能逐步掌握好统计学。

@yanlinlin82

yihui yihui 2022-12-16 06:35:14

加油!

——原帖发布于 2008-11-30 11:12:41

giscus-bot giscus-bot 2022-12-16 05:44:27
访客 *gong* @ 2008-12-02 23:05:28 写道:

谢先生, 你好。

看了你的那个R包,感觉不错,首先祝贺你的成功。

不过发现了一点小小缺陷(其实也不算缺陷),就是无法用你的包名调用help文件,例如,help(animation)?animation。我想你在包里加一个这种函数也很容易。

例如,picante这个包一样。

picante-package(picante) R Documentation 
Picante: Phylocom Integration Community Analyses Null-models Traits and Evolution in R
Description
Collection of functions for integrating ecology and evolutionary biology 

Details
Package: picante 
Type: Package 
Version: 0.2-0 
Date: 2008-06-17 
License: GPL-2 

Author(s)
Steve Kembel <skembel@berkeley.edu>, 
David Ackerly <dackerly@berkeley.edu>, 
Simon Blomberg <s.blomberg1@uq.edu.au>, 
Peter Cowan <pdc@berkeley.edu>, 
Matthew Helmus <mrhelmus@wisc.edu>, 
Cam Webb <cwebb@oeb.harvard.edu> 

Maintainer: Steve Kembel <skembel@berkeley.edu>
yihui yihui 2022-12-16 06:35:14

谢谢您的关注,不过我不知道您要看的是什么信息,是help(package = 'animation')呢还是library(animation); ?'animation-package'
一般来说对于一个函数包,其名称并不是一个函数,所以不能直接查看帮助(例如?MASS就得不到MASS的介绍),不过您要是觉得?animation这种方式比较方便的话,我会在下一个版本中加入这个功能。目前我已经在R-Forge上更新:https://r-forge.r-project.org/projects/animation/

——原帖发布于 2008-12-02 18:42:16

giscus-bot giscus-bot 2022-12-16 05:44:28
访客 *zhoulvjun* @ 2008-12-04 07:23:59 写道:

想对一批记录做综合评价,但变量是定类、定序尺度的,如果因为不是定距变量因而不能用主成分或因子分析,那么该怎么办呢?
请指教,谢谢!

yihui yihui 2022-12-16 06:35:15

1、对于没有顺序的定类变量你怎么评价?
2、因子分析我不懂,不过主成分分析哪个条件说了要求变量必须是连续变量的?
3、你认为主成分和因子分析和综合评价有什么实质联系?(我是问认为,不是问别人怎么说)

——原帖发布于 2008-12-04 10:40:16

giscus-bot giscus-bot 2022-12-16 05:44:29
访客 *zhoulvjun* @ 2008-12-04 14:48:20 写道:

谢谢您这么快就回答了
1、嗯,我还真是糊涂了,如果只是类别的话,还真是没法做出什么评价。那二分变量应该可以吧,比如说同意或不同意,回答同意的人应该可以评价他更好或差一点。
2、如果定序变量可以的话,那是把它视为定距变量来使用么?可我的数据主要是社会态度方面的定序变量,对每个人而言,其满意与非常满意之间的距离并不与其它人的相同,因此不能将其“近似”视为定距变量,虽然现在很多论文里是这样做。如果主成分分析本来就可以用定序变量的话,这与我所知道的分析范例不同,它们都是用定距变量做的。可能是我误解了吧,因为觉得主成分既然用各变量来线性表出,那么,用定序变量就是在将它“近似”成定距变量了。
3、主要是因为变量太多了,而自己了解的降维方法就这两种,而它们又能计算出综合的主成分或因子得分,我就想用这两种得分来做评价。

yihui yihui 2022-12-16 06:35:17

2、主成分分析需要的是数据的协方差阵(或相关系数阵),所以关键在于变量之间的协方差或相关系数怎么算(或估计),这和定序还是定距抑或定类没有必然关系。
3、“得分”这个名词很能忽悠人,我希望哪一天谁能真正帮我解释清楚什么叫得分,至今为止我个人没有遇到这样的人。你也记下这个问题吧,得分和评价有什么实质联系?(别玩儿数学,玩儿实际的东西)

——原帖发布于 2008-12-05 13:14:49

giscus-bot giscus-bot 2022-12-16 05:44:29
访客 *大树* @ 2008-12-05 08:12:29 写道:

谢老师,你好:

本月13日在人大举行的会议若只是旁听,也需要注册吗?

谢谢。

大树

2008年12月5日11:12:24

yihui yihui 2022-12-16 06:35:16

不好意思,已经没有座位了,如果你不介意站着的话当然可以来:)

——原帖发布于 2008-12-05 13:10:26

giscus-bot giscus-bot 2022-12-16 05:46:30
访客 *缰* @ 2008-12-07 19:53:19 写道:

今天下午去听了“一二·九”合唱,晚上上网的时候偶然看到师兄在05年写的一篇博客,恰巧也是关于“一二·九”的。顿时想到了很多东西,而概括起来,不外乎就是“中国青年”四个字。
感觉已经很长时间没有仔细思考过这四个字的含义了,本来想大发感慨几句,却发现自己已经没有发言权了。中国青年在做些什么,我又再做些什么。
05年的时候师兄是大三吧,现在大三的自己,看到那些文章真是有种时空穿梭的感觉。看着师兄在大三时留下的文字,既像是在翻一本泛黄的相册,又像是在对比镜子中的自己。
也许我要开始做些什么了……

@yixuan

yihui yihui 2022-12-16 06:35:18

你肯定会比我成长的快,子曾经曰过,后生可畏啊
另:今年很奇怪,天地人大上面怎么没有因为129打架?现在都变成和谐129快乐129了?

——原帖发布于 2008-12-08 20:40:04

giscus-bot giscus-bot 2022-12-16 05:46:31
访客 *缰* @ 2008-12-10 08:12:46 写道:

相互比较竞争肯定都会有的吧,只不过貌似现在更关注“软实力”了…… 😎

@yixuan

yihui yihui 2022-12-16 06:35:20

原来如此,大不如前了,想当年,那家伙,号称顶上十大一周不沉呐!

——原帖发布于 2008-12-10 12:07:09

giscus-bot giscus-bot 2022-12-16 05:46:32
访客 *huihui* @ 2008-12-10 11:11:54 写道:

还有位置吗?我是武汉的学生,刚好出差到北京.

yihui yihui 2022-12-16 06:35:19

可以,来吧,请把报名表发给我

——原帖发布于 2008-12-10 12:06:04

giscus-bot giscus-bot 2022-12-16 05:46:32
访客 *缰* @ 2008-12-10 18:35:59 写道:

弱问一句,师兄当时是不是唱了超高难度的《蓝色多瑙河》啊? 😁

@yixuan

yihui yihui 2022-12-16 06:35:21

正是……春天来鸟(哈哈吼吼),大地在欢笑(哈哈吼吼),蜜蜂嗡嗡叫(哈哈吼吼),风吹动树梢(哈哈吼吼)……
当年这个129真是烦死人了,为此跟某些学生干部吵架不知吵过多少回

——原帖发布于 2008-12-11 10:28:08

giscus-bot giscus-bot 2022-12-16 05:46:33
访客 *大树* @ 2008-12-11 18:40:30 写道:

谢,

你收到我发给你的短信了吗?

周六我们所去包括我在内的至少2人,或许是3人。

大树

2008年12月11日21:40:13

yihui yihui 2022-12-16 06:35:22

可以来,会场可能比较挤,而且后面报名的都没有安排午餐

——原帖发布于 2008-12-11 19:27:14

giscus-bot giscus-bot 2022-12-16 05:46:34
访客 *大树* @ 2008-12-11 18:43:14 写道:

谢,

代一个学数学的MM问一个弱智问题:

因子分析变量多的时候为什么不出现因子载荷阵 ?

大树

2008年12月11日21:42:59

giscus-bot giscus-bot 2022-12-16 05:46:34
访客 *庞庞* @ 2008-12-12 06:19:45 写道:

我对大概120个数据进行因子分析,spss默认会输出特征值大于1的主成分,最后的 计算结果有30个特征值大于1,但是最后的结果不显示因子载荷阵,只有一个标题,前天我做的时候有19个主成分都显示了,我改为输出特征值大于2的只能输出10个,再多的时候就不再显示了

我现在急需要这个,然后再进行下一步的分析

yihui yihui 2022-12-16 06:35:24

这是SPSS的问题,我不是客服,谢谢

——原帖发布于 2008-12-13 18:34:53

giscus-bot giscus-bot 2022-12-16 05:46:35
访客 *anjing* @ 2008-12-13 17:34:43 写道:

请问:今天R统计会议的文章,在什么地方可以下载

yihui yihui 2022-12-16 06:35:23

全文不会挂在网上,摘要和幻灯片征得作者同意之后会挂出来,COS论坛和统计学院网站都会放

——原帖发布于 2008-12-13 18:31:18

giscus-bot giscus-bot 2022-12-16 05:46:36
访客 *sky* @ 2008-12-29 22:00:44 写道:

请教个问题:在看MASS,感觉Robust Regression和Resistant Regression非常相似,而它们的中译名称(稳健回归、耐抗回归)几乎没啥区别,稳健和耐抗两词差不多是同义词了吧。Robust Regression回归选好psi函数的话应该有不错的耐抗性,而Resistant Regression本来就比较稳健。陈希孺的《近代回归分析》中的稳健性的定义是“具备一定的抗干扰性”,就是耐抗性的意思了吧。这么一来,这两种回归方法的命名方式令我非常困惑,难道其中还有玄机?

当然,Robust Regression和Resistant Regression的区别也是明显的,Robust Regression是一种特殊的加权方式,而Resistant Regression是通过稳健估计量(如中位数对应lms,切尾和trimmed sum对应lts)来估计参数的,它们计算难度也很有区别。

不知老大对这两种方法有何高见?在下竖起耳朵来听……

yihui yihui 2022-12-16 06:35:25

字面上看Resistant是专指对离群点的耐抗,而Robust回归不一定有这个效果。MASS一书中对于稳健回归主要是从M-估计的角度说的,而耐抗回归主要是LMS、LMT和S-估计。

——原帖发布于 2008-12-31 20:17:19

giscus-bot giscus-bot 2022-12-16 05:46:36
访客 *petrel* @ 2008-12-30 17:04:35 写道:

您好,我请教一个R软件上的问题,自己编写了一个函数,想应用lines,画在图上。但不成功,给出提示:错误于as.double(y) :
 cannot coerce type 'closure' to vector of type 'double'
。应用plot却可以画出,请问这是什么原因?迫切,谢谢,
也可邮箱联系:liuhaiyen(at)126.com

yihui yihui 2022-12-16 06:35:26

str(y)

Originally posted on 2009-01-01 19:05:50

giscus-bot giscus-bot 2022-12-16 05:46:37
访客 *zhoulvjun* @ 2008-12-31 17:58:17 写道:

元旦快乐!
牛人在牛年更牛!

Happy 牛 year!

yihui yihui 2022-12-16 06:35:24

这句话很牛……

——原帖发布于 2008-12-31 19:40:20

giscus-bot giscus-bot 2022-12-16 05:46:38
访客 *缰* @ 2008-12-31 21:09:18 写道:

师兄黑皮牛也!

@yixuan

giscus-bot giscus-bot 2022-12-16 05:46:38
访客 *sunshine* @ 2009-01-03 18:27:06 写道:

我发现R很多的程序包功能非常强大,但目前发现无论是刊物还是网络上用中文介绍的非常的少,请问你有没有对程序包做过深入的研究...你是否有关于这方面的信息。..谢谢

yihui yihui 2022-12-16 06:35:27

先看CRAN Task Views,牛包一般都会被收录在那里面,然后就是读R News,编辑看得上的论文也表明那个包有发表的价值,然后就是道听途说,最后当你发现实在找不到你要的东西时就挥笔自己写一个包吧。你这个问题问得太泛了,1600多个包,我自然不可能挨个研究……

——原帖发布于 2009-01-03 19:47:35

giscus-bot giscus-bot 2022-12-16 05:46:39
访客 *sunshine* @ 2009-01-12 21:18:28 写道:

现在市面上的R书都是主要结合统计分析来编,要是你能开拓性编一本详细介绍R中其他的强大功能的书就好咯,R实在太有趣了..非常期待!

yihui yihui 2022-12-16 06:35:28

这种书在出版社恐怕难以过关……

——原帖发布于 2009-01-13 13:29:03

giscus-bot giscus-bot 2022-12-16 05:46:40
访客 *小黎* @ 2009-01-30 15:29:20 写道:

新年快乐,打个招呼,呵呵

yihui yihui 2022-12-16 06:35:29

新年好!

——原帖发布于 2009-02-10 17:11:39

giscus-bot giscus-bot 2022-12-16 05:46:40
访客 *gtly* @ 2009-02-09 18:08:07 写道:

一直逛统计之都和这个博客,非常喜欢,学到很多方法论,受益匪浅!牛年希望能继续看谢老师大作! 😳

yihui yihui 2022-12-16 06:35:29

多谢支持!

——原帖发布于 2009-02-10 17:11:29

giscus-bot giscus-bot 2022-12-16 05:46:41
访客 *StimetCross* @ 2009-02-18 10:48:03 写道:

好久不见。最近过的好吗?来拜个晚年
我换日志了 一年多没写了。呵呵 有空来看看:)

giscus-bot giscus-bot 2022-12-16 05:46:42
访客 *lynn* @ 2009-03-02 08:33:51 写道:

谢老师,在网上查到您能为我们提供一些统计软件的令人高兴的好消息,只可惜不知如何与您取得联系,谢老师,我是武汉大学研究生,目前正在寻找Amos5软件,期望能得到老师的帮助,先行谢谢啦~祝您一切顺利!我的邮箱是fuzhiyun3087@126.com

yihui yihui 2022-12-16 06:35:30

请问你在哪里看见我本人说可以提供软件的?

——原帖发布于 2009-03-02 10:26:00

giscus-bot giscus-bot 2022-12-16 05:46:42
访客 *chenwen* @ 2009-03-04 07:54:17 写道:

您好,请问,您可以给我传一个AMOS的软件么?不是不想自己去弄,实在是下不了.拜托啊 .请您加我吧:QQ:962118512

yihui yihui 2022-12-16 06:35:32

http://yihui.name/cn/post/238.htm

Originally posted on 2009-03-04 12:40:03

giscus-bot giscus-bot 2022-12-16 06:35:46
访客 *jiaojiaojing* @ 2009-05-12 11:56:06 写道:

新版面着实让我吃了一惊,先研究了研究怎么留言。

giscus-bot giscus-bot 2022-12-16 05:46:43
访客 *chenwen* @ 2009-03-04 07:55:18 写道:

发完了,才发现别人都叫您谢老师,不好意思啊,之前发的也没什么礼貌.

yihui yihui 2022-12-16 06:35:31

我不是老师……

——原帖发布于 2009-03-04 12:38:37

Sign in to join the discussion

Sign in with GitHub