留言 | /cn/about/

哇哈哈哈，老夫昨天可是费了一整天的功夫寻思怎么把主站重新排版一下啊

——原帖发布于 2009-03-12 09:31:30

此时xlim和ylim由设备窗口大小（宽高）决定，你把窗口拖窄一些就会发现x轴范围变了。也可以在作图之前就规定设备的宽和高，如：

windows(3, 6)  
plot(1:10,11:20,xlim=c(0,10),ylim=c(0,20),asp=1)  
# 或者png()等设备

——原帖发布于 2009-03-15 12:22:58

抱歉，我无法回答，我不知道分数线是多少，也不知道其他人情况怎么样……

——原帖发布于 2009-03-17 07:56:14

多谢相告，不过第一个会议注册费也太贵了，第二个会议也没说是否收费……

——原帖发布于 2009-03-17 15:45:05

安装animation包：install.packages("animation")
然后运行网页中的代码（生成HTML动画网页），或者?sample.strat示例代码
或者安装ImageMagick程序，这样可以把动画转换为GIF格式插入PPT
或者安装SWF Tools程序，把动画转换为Flash插入PPT

——原帖发布于 2009-03-17 19:55:42

您是？……

——原帖发布于 2009-03-22 19:23:27

需要将这些程序的路径加入系统PATH中

——原帖发布于 2009-03-22 19:23:08

dotchart函数和你有同样的问题，我看它是用plot.new()而不是plot(..., type = 'n')解决的，不妨参考一下。

——原帖发布于 2009-03-26 11:05:46

我陷入了统计图形的“局部最优解”：感觉什么图形都可以用点线等元素去做，但事实上肯定有更方便的方式，例如grid、lattice和ggplot2我都没有仔细去学习

——原帖发布于 2009-03-26 19:41:03

先解决数学问题，再解决R问题。
z=f(x,y)=z0，求x和y的函数关系。

——原帖发布于 2009-04-03 14:22:16

Efron: An Introduction to Bootstrap

Originally posted on 2009-04-07 10:32:04

凸包怎么定义？

——原帖发布于 2009-04-07 10:31:42

那就写一个lattice版本吧，你的图中都被圆圈占满了，作图区域中没有地方放图例了，只能放到边界中。

——原帖发布于 2009-04-09 16:18:46

看到了，谢谢啊

——原帖发布于 2009-04-09 09:37:40

指教不敢，互相学习：）

——原帖发布于 2009-04-16 17:16:48

http://cos.name/bbs/thread.php?fid=15

Originally posted on 2009-04-19 19:01:12

麻烦事还有一大堆呢……清理内容、调整格式……最重要的是把以前的链接“无缝”过渡到WP中来，否则那些旧链接都失效了就可惜了

——原帖发布于 2009-05-12 12:05:07

没啥没啥，多提意见，俺好好改进，谢谢 :twisted:

——原帖发布于 2009-05-19 06:50:23

已经发了，希望你清楚怎么处理shape file 🆒

——原帖发布于 2009-05-22 07:48:39

不熟悉，你说的是SQL Server吧？某个分析工具？我从不用数据库工具分析数据。

——原帖发布于 2009-05-22 08:04:46

访客 *唐卡* @ 2009-05-25 05:34:24 写道：

那我只好继续想办法~`

抓紧最后的两个月学习时间吧，上了研究生就没空学习啦，我看硕士童鞋们一天到晚都在写作业，忙倒是忙得很，学到什么就不好说了。这个现象很奇怪，不知是谁的问题。要说推荐书嘛，因为我不太了解你的背景，所以也不好推荐，数理统计方面找本中文书看就差不多了，看个梗概思想就足够，我可以以过来人的身份武断地说，你现在想全面补数理统计的基础基本上已经来不及了。数学方面若重看高代数分必然太慢，不如直接找本多元的书看附录（如张尧庭、方开泰《多元统计分析引论》的第一章和附录），熟悉对统计有用的矩阵操作；（如果你不是纯理论研究者的话）统计方面离不开实践，尽看书不如不看书，不做点事情的话，将来拿到数据两眼一抹黑，肯定也不行，而实践又离不开软件，说到软件的话在我这里就没什么别的软件可选择了 😁 因此归根结底，我建议你啃一啃MASS这本书（Modern Applied Statistics with S）或者其它跟R语言有关的书，统计之都下载站的Statistical Software目录下面有大量的电子书，你可以在暑假里好好学一门软件，在玩数据的过程中练模型和方法。祝你顺利 :lol:

——原帖发布于 2009-05-23 13:13:31

访客 *kathy* @ 2009-05-23 16:02:43 写道：

谢谢师兄~我觉得软件还是得好好学一下
fighting～～～

嗯，加油加油！ 😉

——原帖发布于 2009-05-23 17:16:39

谢谢，可是花了俺好几天的功夫啊，界面大为简化了，但担心读者找不到日志的全文链接…… :twisted:

——原帖发布于 2009-05-24 18:57:11

研究一个问题，首先要给概念下一个清楚的定义，比如你的问题中，什么是“好”的定义？如果“排名越靠前就越好”的话，那么你根本不需要管什么长尾分布，把所有歌手的数据拿来，看你关注的歌手在卖碟和卖书方面分别排名第几，马上就知道答案了。

“标准分”的出处我从考证，可能源自 $(x_i-\bar{x})/\sqrt{Var(x)}$，个人觉得这个概念并没什么用，前面我写过一篇日志“这个年代该怎样教统计（或不该怎样教）”，今天想想这个概念，觉得也应该扔进垃圾堆。在计算分布概率不方便的年代，标准分可以当作概率的替代工具，因为分位数和概率基本上是等价的概念，如今计算机这么方便，还计算标准分干什么？要作比较的话，直接比较概率则一目了然，概率的尺子比分位数要易懂得多。

当总体分布未知的时候，常用经验分布代替，此时概率和秩（即排名）基本上等价，所以绕了一大圈，最终仍然是看排名就够了。注意这里说的概率一般情况下是尾概率，当然， $P(X < x)$ 和 $P(X > x)$ 也是等价的。

——原帖发布于 2009-05-26 08:50:40

访客 *bnudahu* @ 2009-05-26 11:39:21 写道：

呃，算你回答了吧 😢
至少作为一个非统计学出身的外行，我基本上明白你的意思了 💋
谢谢

看来我文不对题了 :?: 对不起呀……

——原帖发布于 2009-05-26 15:52:23

因为相关部门要强行备案，COS宁死不从，现在搬到国外去了 😁

——原帖发布于 2009-05-31 15:17:03

我想在COS设立一个奖项，奖品就是网站空间，COS现在有能力提供足量的博客空间了 :lol: 不过还没想好组织什么竞赛或活动…… :mrgreen:

——原帖发布于 2009-05-31 15:49:22

访客 *魏太云* @ 2009-05-31 16:09:01 写道：

啊，COS这么强大了 😄

我也有个主意，就是用不着为了资金问题在COS上面放广告（虽然我很少看见，主要是firefox下的ABP过滤了），COS靠大家捐助也不是不光彩的事情，买空间的费用也不太贵，大家捐点都是心甘情愿的，我下学期得了奖学金一定抽一份子出来报效COS 😄

@taiyun

COS的广告会尽快撤掉，只保留下载站的广告，论坛和主站的广告即将消失。将来我想COS的费用应该主要来源于通过COS本身得到的各种统计项目，而不要依靠捐赠。学生如果从COS项目中获得收益，可以考虑一定比例捐赠，如果是奖学金之类的钱就算了，这个事情我说了很多次了，不提倡学生捐赠。

——原帖发布于 2009-05-31 18:15:13

又來了一位UoA的學生啊，歡迎歡迎，哈哈
我早已經把推廣R作為我的常規任務了 :lol:

——原帖发布于 2009-06-03 19:18:22

邮件、留言、评论任选一种方式告诉我即可，请不要三种都用。谢谢。

——原帖发布于 2009-06-07 11:23:39

请附上数学公式，谢谢。从公式翻译到代码容易，从代码翻译回公式很要命。如果你会用LaTeX，请用**$latex** latex formula**$**写公式。

——原帖发布于 2009-06-22 19:38:16

访客 *Samson* @ 2009-06-25 05:23:41 写道：

非常感谢您能回复!

这两天我在试的时候还看到BB里的dfsane()也能求解非线性方程组,不过奇怪的是我要求解的30次方程组中,有10来次是可以收敛的,另外的那些即使我多次把求解结果代回去做为初始估计值重新求解也无法收敛.

所以我现在的问题是:optim(), nleqslv(), dfsane()哪个可以解决我的问题?应该怎么用呢?
如果您还知道哪里能找到这方面的具体介绍也请指点一下,网上实在很难找到这方面的详细介绍.
麻烦您了!

下面是我要求解的方程组 (在我的例子里 $\tau$ 为1,所以省去了):

%下面是第一个方程:
%期权定价公式
$S_0=V_0\Phi(d_1)-De^{-r\tau}\Phi(d_2)$

其中：
$d_1=\frac{\ln(V_0/D)+(r+\frac{1}{2}\sigma_V^2)\tau}{\sigma_V\sqrt{\tau}}=\frac{\ln(V_0/De^{-r\tau}+\frac{1}{2}\tau\sigma_V^2)}{\sigma_V\sqrt{\tau}}$

$d_2=d_1-\sigma_V\sqrt{\tau}$

%下面是第二个方程:
$\sigma_S=\Phi(d_1)\frac{V}{S}\sigma_V$

在已知 $S$ 、 $\sigma_S$ 、 $D$ 、 $r$ 和 $\tau$ 的情况下，联立方程一和二求解公司价值 $V$ 及公司价值收益率的标准差 $\sigma_V$ 。

optim()是用来做优化（求极值）的，你这里是解方程组，不一样。我没用过nleqslv()函数，刚看了一下帮助文档，和那些经典的优化函数也差不多，无法收敛可能有很多种原因，你可以试试（1）换参数初始值；（2）换求解方法；（3）换搜索方法；（4）换迭代次数；（5）换步长；……等等等等

——原帖发布于 2009-06-25 13:58:24

懒得分页了，一页看起来方便，省得翻 😁

——原帖发布于 2009-07-22 14:19:55

我马上要去Iowa State University读PhD了

——原帖发布于 2009-07-22 14:20:33

对Q没有研究，抱歉无法评论……

——原帖发布于 2009-08-02 00:43:17

过奖了，哈哈

——原帖发布于 2009-08-11 01:58:45

已发。

——原帖发布于 2009-08-18 06:55:47

Sounds like an pleasant trip... And I'm glad to know your successful application.

I've sent the document you requested by email. Thanks for pointing out the problem!

Originally posted on 2009-09-14 08:40:08

我不记得我曾经做过SAS的链接……

——原帖发布于 2009-09-24 07:57:29

标题都在右边……呃，我这个设计确实不够易用，不过懒得改了，谢谢 :lol:

——原帖发布于 2009-10-04 20:29:19

访客 *华晨* @ 2009-10-04 21:10:34 写道：

哦，呵呵，原来这样呀！不说还真发现不了。
不过这也算是创新了，不过就是不太实用，呵呵。挺不错的

谢谢，我们有自己的服务器，现在不缺硬件资源。

——原帖发布于 2009-10-16 00:10:54

为了保证时间的完整性，我几乎不用QQ，也不喜欢即时聊天。办网站如你所说并不容易坚持，统计之都也坚持了几年了，暂时没有合并的需求。谢谢你的关注！

——原帖发布于 2009-10-16 02:18:01

用R把剩下的号码组合写进文件，调用system()执行md5sum算一下就知道了。不过其实要是真想知道我的QQ号，动动脑子就找到了，放md5sum只是个幌子 😁

——原帖发布于 2009-10-19 06:46:03

邮箱是对的，邮件已收到。

——原帖发布于 2009-10-29 20:51:07

我在Iowa State Univ统计系，来吧来吧，一定盛情招待，虽然我们这个农村什么都没有 😁

——原帖发布于 2009-11-03 01:25:33

不知道……没用过你说的那个软件，我都是用Google Reader……

——原帖发布于 2009-11-03 07:23:17

速度够快的啊，我才刚刚开始整容，正在踌躇中用哪个模板呢……前一个模板很多观众朋友们都抱怨说找不到文章标题在哪儿或者两篇文章之间不知道从哪里断开的，这两天感恩节放假，正好整一整

——原帖发布于 2009-11-24 06:56:04

我一般用Firefox，没有发现太占用资源的情况；你的CPU貌似确实老了点啊。不分页也是为了浏览方便，免得一页一页翻。

统计入门的书我也不好推荐，我自己是从概率论数理统计一步步过来的，不知道你想学到什么程度。

——原帖发布于 2009-11-27 10:05:06

你在R的主页上搜一下邮件列表就知道了，我刚看了一眼，至少有HiddenMarkov、msm、RHmm等包都可以做，质量如何自己把握。

——原帖发布于 2009-11-29 20:57:01

今天无意瞅了一眼后台被过滤的垃圾评论，发现你的这条评论被误判为垃圾评论了……看来Akismet不完全可靠啊

——原帖发布于 2009-12-05 00:43:11

参见楼上留言148。md5是个结识人的过滤器 :lol:

——原帖发布于 2009-12-06 07:22:25

访客 *views63* @ 2009-12-06 07:50:23 写道：

没看前面的留言，原来 R 也可以计算 md5sum
对 $epsilon$ 表示$forall epsilon >0$的学科感兴趣，我的 QQ 也很少用，讨论问题不如论坛、邮件（列表），真的要说什么不如打电话。

（抱着试试的心态，自己能力破解用了别人写的工具。）

md5sum可以用R调用系统命令md5sum生成（system("md5sum somefile")），也可以用tools包中的md5sum函数生成，开源软件大多都能相互支持

你是学数学的啊？以后遇见数学题就找你了啊！

——原帖发布于 2009-12-06 08:08:36

访客 *views63* @ 2009-12-06 14:47:08 写道：

应用数学，只是学校的课程数学类的还不到一半……
数学学得很差（但很感兴趣），人也一点不优秀。
如果有数学问题，很乐意一起讨论。

当然不是我自己写的啊，WP的模板千千万，没有必要自己写，我也没那个精力。你看一下网页源代码中CSS文件的路径就知道模板的名字了。我一般不用那些被WP推荐的模板或者最流行的模板，因为用的人太多了。

你的嵌套貌似不是真的嵌套，而是一种古老的回复方式，用@表示对谁回复。WP自从2.5还是哪个版本开始就自身支持嵌套回复了，不需要用插件或者模板实现。我用的就是WP自己的嵌套，在后台可以打开嵌套回复的选项，也可以设置嵌套深度。这一点陈丽云肯定深有体会，哈哈 :lol:

另外，你在本站的URL写错了，http:后面少了一个/。

——原帖发布于 2009-12-08 23:15:09

访客 *魏太云* @ 2009-12-09 07:51:46 写道：

丽云说嵌套深了就不好管理了，她还建议我自己写模板呢。

@taiyun

我自己从来都是改模板，好几年前做网站才自己写模板，后来觉得这里面技术太繁杂了，还是让那些疯狂的程序员和设计师去干这种事情吧，等他们写完拿来改改就好了。不过写模板会让你学到各种网页旁门左道，你每时每刻都要发愁IE和FF以及其它浏览器的兼容性、考虑用户的屏幕分辨率、字体、操作系统，甚至是否支持JS以及Flash的版本。设计模板的人肯定都有自己的烦恼，但我想每个设计过模板的人一定有个共同点，那就是极度痛恨那白痴IE6在世上祸害这么多年仍然阴魂不散（逼得别的浏览器必须用<!--[if IE 6]>这样的语句专门服侍IE6）。

嵌套不要太深就行，我设置的是4层，看起来还凑合。

——原帖发布于 2009-12-09 08:08:06

访客 *魏太云* @ 2009-12-09 08:53:40 写道：

我换了一个模板，嵌套可以了；不过这个模板似乎商业味道很浓。嵌套也是4层:)

@taiyun

访客 *魏太云* @ 2009-12-23 07:34:28 写道：

wordpress真能折腾人，现在选取了一个比较满意的模板，小地方的细节也修改了不少。不过不幸发现了一个大毛病，就是点击分类目录下的元素(比如R)时，生成的 http://taiyun.cos.name/category/r/ 页面中的格式都丢了，包括图片、各级标题等，不知怎么解决？

@taiyun

这个页面的生成一般是靠模板中的archive.php控制的。从你的症状来看，肯定是输出页面列表的时候在什么地方用了PHP函数strip_tags()函数。你对照single.php或者index.php看一看就知道输出原理了，一般都是用WP函数the_content()输出。

——原帖发布于 2009-12-23 07:43:28

访客 *liyun* @ 2010-09-30 03:19:57 写道：

我，我，我怎么时隔近一年之后才看到……
我，之所以设置一层嵌套只是因为我懒得写多层嵌套的CSS，检讨去了……

@cloudly

这就是懂技术的人的痛苦啊，什么都可以自定义，但经常走到一个不能通用的角落，等到核心技术更新了，自己更改的东西是放弃也不好，不放弃也不好。

——原帖发布于 2010-09-30 07:19:18

不如直接找个程序员谈吧，我只是业余玩网站而已，对正式的网站工作兴趣不大。我个人的建议是你可以看看 wordpress.org 上的那些案例。

——原帖发布于 2009-12-15 21:50:22

我没用过nleqslv，刚才看了一下帮助文档，它没有设定x范围的参数。对于有限制的优化问题，如果问题相对简单的话，可以参考?constrOptim或者?nlminb。

——原帖发布于 2009-12-23 02:24:03

你的问题问得太不清楚了（比如我问你究竟什么叫bagging里面变量的权重？）。参见Ihavenothing在论坛给你的回答。另外，论坛发帖和给我留言这两种途径请任选其一，因为我两边都会看，谢谢。

——原帖发布于 2010-01-01 12:11:15

访客 *hanhui* @ 2010-01-02 00:43:47 写道：

师兄，其实问题是这个样子的，有1到50个同学，用X1-X4这4个variable来描述，根据这4个variable可以把50个学生分成ABCDE5组，我要做的是，设立一个benchmark，根据这个benchmark而给每个人打分，从而能够用ABCDE来划分这些学生。我本来想用classification做这个，因为R里bagging这个pacakage有importance这一项，但又来想想也不对。。。一直很纠结这个score怎么弄出来，麻烦师兄帮我解答一下，不胜感激

跑题一下：呃，请问你是不是出国了……这中英混杂打字难道不嫌切换输入法麻烦么……纯属个人偏好：我是个土人，很不习惯看到或者听到别人把一些常见的中文词语整成英文夹在中文里。跑题结束。

窃以为这问题还是没说清楚：

“根据这4个variable可以把50个学生分成ABCDE5组”：“可以”是什么意思？怎么分组？

“设立一个benchmark，根据这个benchmark而给每个人打分，从而能够用ABCDE来划分这些学生。”：benchmark具体是什么？打分又如何“打”？前面根据变量把学生分了ABCDE五组，后面怎么又用ABCDE来分学生？

“因为R里bagging这个pacakage有 importance这一项，但又来想想也不对。。。一直很纠结这个score怎么弄出来”：R里没有bagging这个包（我是该装糊涂不知道呢还是该纠正为“R的adabag``包里的``bagging``函数”）。总之你说了半天，我仍然不知道你要的“score”到底是个神马东西，抱歉，阅读能力有限。

——原帖发布于 2010-01-02 02:53:57

访客 *hanhui* @ 2010-01-02 06:41:31 写道：

谢谢师兄的回复，我重述一遍吧：
有1到50个同学，用X1-X4这4个变量来描述，根据这4个变量，题设已经把50个学生分成了ABCDE5组，其中组A要好于组B要好于组C。。。我的目的就是要找出，为什么要这么分组：要设立一个评分标准，根据每个学生不同的变量，从而能够算出50个同学不同的得分，利用这个得分把50个学生重新进行划分，使得这个划分与题设已给划分相同
数据例子：
姓名 组别 年龄 身高 体重 臂展
甲 A 22 176 120 100
乙 B 30 165 111 87
丙 A 26 200 190 120
.
.
.
我想了好久不知道该用什么方法来解决这个问题，麻烦师兄帮我看看，万分感谢

访客 *yixuan* @ 2010-01-02 06:58:13 写道：

好吧，我来凑个热闹。首先你要判断一下那个事先给定的分组是不是真的有依据，搞不好是别人随便说说呢。如果你只是想验证或怀疑一下那个分组是不是合理，你可以用聚类的思路去做，也就是你把这50个学生聚成4类，然后看看聚出的类与给定的类是不是吻合就行了。

@yixuan

访客 *yixuan* @ 2010-01-02 06:59:15 写道：

呃，没注意，应该是5类。

@yixuan

访客 *yixuan* @ 2010-01-02 07:21:31 写道：

聚类就是一种探索数据分组的方法，如果你通过已有的分类建立了一个指标，然后又去用这个指标去验证分类的合理性，这不是因果循环吗？

@yixuan

访客 *hanhui* @ 2010-01-02 07:32:49 写道：

请问，那如何才能给出一个标准，证明我的聚类方法能够验证原先分类的正确性。 聚类只能得到一个样本之间的距离，并不能给每一个变量一个权重，从而使每个观察值得到一个综合评分，这是我对聚类的理解，如果你能用聚类做，麻烦告诉我一下，不胜感激

访客 *yixuan* @ 2010-01-02 08:12:40 写道：

要不然这样，把验证分组的过程看成两步，第一步是考虑无序的分组是否合理，第二步再考虑组与组之间的大小关系，如果第一步就推翻了那有序的就更不用说了。第一步可以先用聚类来做，第二步可以用主成分来得到每个观测值的得分，然后比较每个类得分的均值。当然我觉得单纯从数据来得到得分不是特别靠谱，毕竟方法本身是不知道实际的意义的。
（我们是不是快成版聊了……）

@yixuan

访客 *hanhui* @ 2010-01-02 10:04:20 写道：

效果还是非常差。。。唉，郁闷了。。。等待益辉师兄

为啥要等我呢，我又不比Yixuan高明，你可不要低估这个$pi$。关于这个出题人的目的，我只能接受一种解释，那就是他故意跟你玩躲猫猫。除此之外，他就是在浪费时间。

即使他要玩躲猫猫，你也未必要纠缠在什么得分上（抱歉你又撞在我不喜欢的东西上了，我很讨厌统计里面”得分“这个东西，尤其是什么加权平均），为什么要有得分才能分组呢？因变量虽然是有序的分类变量，但有序分类也是分类啊，先用无序的方法试试，如果无序的方法能正确分组了，那就不用管这因变量到底是有序还是无序了。（我怎么跟说绕口令似的）

只要能找到办法分开这些学生，也就是有了规则把谁分到哪一组中，你说此时需要什么得分么，一定要把变量线性组合么？

最后，”效果非常差“要说清楚了（我三番五次请你把话说具体一些，给别人尽量多的信息，不要等人来挤牙膏），用了什么方法、得到了什么效果。按你给的三个数据，肉眼能看出来：年龄小则A（否则B）、身高高则A、体重大A、臂展长A……

——原帖发布于 2010-01-03 11:48:50

你这么一提醒，貌似俺该更新简历了。人都是在不nb的时候拼命装nb，nb的时候拼命装不nb……所以不要看那华丽丽的列表，那里面一半是扯淡，另一半是瞎扯淡。

照片嘛，某小子几乎从不把正面清晰无码的活人照片放网上，但似乎无照片不成简历，所以就放个轮廓好了。

——原帖发布于 2010-01-25 20:53:57

我也很少煮面，所以也没啥高方啦，只是觉得煮面时加几滴油并不能防止溢出而已……我若要吃面条就自己擀，还没关注过牌子……

——原帖发布于 2010-01-26 20:45:57

访客 *焦静* @ 2010-03-04 13:45:11 写道：

啊？自己擀面啊，你牛！我是很爱吃面，但不会做，每天都喜欢去食堂牛肉拉面窗口蹭蹭。

煮面条时看它快溢出来了就加一小碗凉水，然后等它再煮开，快溢出来了再加凉水，反复加四五次面条也就差不多好了。（最近连续煮了几天面条 :-) ）

——原帖发布于 2010-03-07 21:10:49

1、因为中文经过编码了，这样存储和访问更“安全”，有些服务器未必直接支持中文文件名。
2、这个系统本身没有限制，但上传文件肯定受PHP的配置限制，一般PHP会限制文件大小为2M，如果你没有修改php.ini的权限，那么就只能按默认大小上传了。
3、不知道什么意思。插件未必带工具栏。

——原帖发布于 2010-03-07 21:08:35

访客 *sunshine* @ 2010-03-08 11:41:43 写道：

对于回答1：中文经过编码后，如何识别在系统中的中文文件夹和文件呢？
对于回答2：如何修改php.ini达到增大上传文件大小的目的？目前正在组建一个用于小型办公的dokuwiki。
对于回答3：明白您的意思了，就是说插件未必有工具栏按钮，如果希望有，是否要编程？如何编写？

新问题4：为何indexmenu插件时而出现在工具栏中，又时而没有呢？
谢谢。。。

访客 *sunshine* @ 2010-03-08 19:01:40 写道：

问题2已解决，谢谢!

1、如果你熟悉PHP的话，参见函数urlencode；至于如何识别中文文件名，我没有100%有效的办法，也许可以urldecode之后再判断，或者判断文件名中是否有%（我不确定）。
3、参见Dokuwiki的开发手册，我不熟。
4、我不熟悉这个插件……

——原帖发布于 2010-03-10 21:26:31

这问题很好，但我几乎无法回答，因为我还没这么深的资历。我在论坛上给版主们发了一个问题，和你这个问题有相通之处，你不妨先看看：http://cos.name/cn/topic/101086

我的感觉是，你不必担心统计人的出路。只是玩玩软件输出的人不可能取代我们的。

——原帖发布于 2010-03-10 21:53:24

抱歉我对生存分析不了解。

——原帖发布于 2010-03-24 18:43:43

访客 *lily* @ 2010-03-25 04:49:24 写道：

哦 还是谢谢了

我真是佩服你这嗅觉……

这个包的第一个版本还没达到我的预期，所以暂时没吱声，昨晚花了几个小时补充了一下，现在总算差不多了

——原帖发布于 2010-03-25 20:33:38

这是我的荣幸啊，谢谢

——原帖发布于 2010-03-28 07:59:05

我当时就给你回信了，没收到？

library(maps)
library(mapdata)
map('china')
# 例：线
lines(c(95, 100,120), c(35, 30, 45), lwd = 2)

——原帖发布于 2010-03-30 17:07:01

抱歉我不了解……

——原帖发布于 2010-04-01 17:38:16

呃，这问题问我其实不太合适。因为如果要我介绍经验的话，我就问一句：你能4年时间几乎每天都开着R么？

无他，唯手熟尔。也许你没这么大兴趣，也许你没这么多时间。我个人没有看完过任何一本书，大部分时间里我都在看帮助。要说起步，如果仅仅是编程起步，那么了解R的几种数据结构就够了，如果是要做统计，那么还有漫长的路要走，因为统计没搞懂的时候那些R函数几乎处处会遇到难题。

——原帖发布于 2010-04-10 10:00:25

你考虑一下QQ号的限制条件之后可能也没那么难（数字、不超过10位），暴力方法都能解出来了。比如： http://cos.name/cn/topic/101321#post-202299

——原帖发布于 2010-04-11 20:26:58

访客 *Linlin Yan* @ 2010-04-12 03:05:45 写道：

奇怪，md5是用私钥或公钥进行加密的么？难道是我记错了？

@yanlinlin82

这种情况我建议直接访谈，不要用什么logistic回归之类的统计模型（如果你的目的是找出挂科因素并提出预防方案而不是仅仅为了发论文的话）。原因我一会儿另发文解释。

——原帖发布于 2010-04-13 20:25:16

访客 *流浪荡云儿* @ 2010-04-13 20:51:09 写道：

好的，顺便说下压根没敢想过发论文呢，只是老师布置的期末作业----------选择合适的题目，收集相应的数据，建立统计模型，进行统计分析，直到你认为满意为止（咱们老师的原话）。是不是我的题目选的不好啊？

访客 *流浪荡云儿* @ 2010-04-14 20:50:10 写道：

为什么是直接访谈呢？

理由见于：
http://yihui.name/cn/2010/04/knowledge-and-common-sense/

——原帖发布于 2010-04-15 09:11:21

我在原文中没说清楚么……数据生成过程和画图的代码全都在原文中： http://yihui.name/cn/2008/10/least-median-of-squares/

——原帖发布于 2010-04-15 10:05:08

访客 *张志坚* @ 2010-04-15 10:38:53 写道：

谢益辉博士,我是在网页上看到的，没有看到原文。另外能否将常见的稳健回归方法中M估计法的人工计算过程发一些资料和实例。

访客 *张志坚* @ 2010-04-15 10:58:49 写道：

数据生成代码看不懂，我没有用过统计软件，但想搞懂这些统计工具。

不熟悉R的话，解释起来就又费劲了……我还是自己替你生成数据附在那篇文章下面好了。

你说的“人工计算过程”是什么？

——原帖发布于 2010-04-15 11:41:57

那么请你指出你所看的材料的出处吧，我不知道你在说什么……

——原帖发布于 2010-04-15 18:39:46

参见148、157、177楼。我知道那个解密的事情，如前面所说，我放个MD5在那里的真正目的不是让人去解密，而是设置了一道“过滤器”，如果有人能加我QQ，我起码知道这个人满足条件（1）真的懂破解MD5，或者（2）遇到问题能仔细想想是不是真的无解，或者（3）会用搜索。

——原帖发布于 2010-04-23 22:04:31

我见你在很多地方问过同样的问题，你说的这两种检验都太含糊了，在普通的线性模型中，t检验和F检验都是精确的，只有在非线性或者混合效应模型或其它类似非线性的模型中，才会出现近似的情况。我所知道的是R的很多模型建好之后都可以用anova()去做F检验（通常是比较嵌套模型）。如：

set.seed(123)
z = x = (1:100)/10
y = 100 + 2/x + z + rnorm(x)/10
mod1 = nls(y ～ Const + A/x + Z*z, start = list(Const = 50, A = 2, Z = 1))
mod2 = nls(y ～ Const + A/x, start = list(Const = 50, A = 2))
anova(mod2, mod1)
## 输出
Analysis of Variance Table

Model 1: y ～ Const + A/x
Model 2: y ～ Const + A/x + Z * z
  Res.Df Res.Sum Sq Df Sum Sq F value    Pr(>F)    
1     98     645.44                                
2     97       0.82  1 644.62   76292 < 2.2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

——原帖发布于 2010-05-01 21:15:41

谢谢。这年头柯南都四处打酱油了……

——原帖发布于 2010-05-03 07:24:51

这名字好，Evan Zhou，莫非江湖人称“一碗粥”……闪……

——原帖发布于 2010-05-11 06:04:17

访客 *流浪荡云儿* @ 2010-05-15 19:48:47 写道：

从理论上说，两组相互独立数据的pearson相关系数应等于0。能设计一个实验进行验证吗

这个……本人谢绝调戏……谢谢！

P.S. 最近大家怎么都这么淘气呢？

——原帖发布于 2010-05-16 05:13:38

生成相关系数绝对为零的两组随机数是很难的，同时也是毫无必要的，你只需要知道理论上相关系数为零就够了。随机数只是一种理论实现。

> set.seed(123)
> cor(runif(1000000), runif(1000000))
[1] 2.255649e-05
# 相关系数够小么？

——原帖发布于 2010-05-11 21:28:43

这老师真调皮……好，那你请你老师生成一个服从伯努利分布的随机数先，并让他/她证明生成的确实是个随机数。

——原帖发布于 2010-05-12 03:05:13

它在我的TODO list上，但近期内我不会深入研究。目前为止我了解也很有限。

——原帖发布于 2010-05-14 13:01:08

访客 *江户川柯南* @ 2010-05-15 05:42:40 写道：

QQQ=3Q,谢啦～

（只是我死活没找到你的TODO list～
=_=!!!，给个链接吧~不甚感激.....）

鼠标一拖到顶，看FAQ3。我没用过这个包，你问问别人吧。

——原帖发布于 2010-05-14 13:09:55

居然是你……更没想到你去nz了……一晃都五年了啊

——原帖发布于 2010-05-14 15:32:42

你说的这个PCO我不熟悉，搜了半天也没发现什么好结果，结果上维基百科一看，原来是Multidimensional scaling（MDS）的前身……那就简单了，MDS就是把高维数据在低维坐标下展示的一种方法，使得在低维坐标中，点与点的距离尽量跟高维坐标中的距离相似。

http://en.wikipedia.org/wiki/Multidimensional_scaling

——原帖发布于 2010-05-16 10:57:56

访客 *冰裹冻* @ 2010-05-16 12:35:35 写道：

嗯，谢谢，茅塞顿开，这个网站真的很有用。
我在用R写有关MDS 作业，被卡住了，可以指点下么？我把详细内容和我写的代码发到你的邮箱里了。拜托拜托。

我的时序只是已经差不多忘光了，而且我现在对时序也很不感冒。印象中，自相关图和偏自相关图分别用来判断MA和AR模型的阶数（也许要反过来，不记得了），R中，acf()和pacf()函数分别可以画这两种图，其中有95%的置信区间（默认蓝色虚线），落在这个区间内的可算是截尾了。

若不是跟我个人有关的R问题，请到论坛提问，见上面FAQ3。

——原帖发布于 2010-05-17 12:48:32

访客 *冰儿* @ 2010-05-17 13:02:02 写道：

恩，是的，应该去论坛。
具体准则教材中很清楚，但是具体样本不会出现理论中的完美情况，有很强的主观意思，不好定夺~·

我不知道你到底是要写论文还是写作业，我觉得没人写论文会写整个“非参数统计”。而变系数模型（如果我理解正确）似乎也并不是非参的内容，“非参用于经济”听起来又像是计量经济学，所以我彻底糊涂了。

变系数应该是混合效应模型的一种吧？如果是这样，lme4包应该是最合适的选择了。我没做过相关的例子。若你不是自谦的话（菜鸟），我觉得把lme4用明白还是有一定难度的。

——原帖发布于 2010-05-19 11:04:58

生成相关系数为0的数字当然是很简单，比如

> cor(c(-1,-1,1,1), c(1,-1,1,-1))
[1] 0

但你能说服别人这是两组随机数么？

你老师的原题是

2、从理论上说，两组相互独立数据的pearson相关系数应等于0。试设计一个实验进行验证和检验并通过ppt在课堂上进行讲解（提前把ppt发到我的邮箱并通过审查），你可直接获得《统计软件使用》课程的免修机会（60分），或者取得该门课程10分期末考试加分机会。

这和你传来的话似乎是不一样的吧？他没让你生成相关系数为0的两组数据，只是让你验证和检验这个相关系数为0。既然这题事关重大，我也没有你老师的授权说可以提供帮助，所以就不发表意见了，省得他找我麻烦。

——原帖发布于 2010-05-20 08:06:12

我刚发现你的老师是COS论坛会员！所以你赶紧卷铺盖逃吧，没准儿他过两天就看到这里来了！ :lol:

——原帖发布于 2010-05-20 08:21:00

访客 *流浪荡云儿* @ 2010-05-22 20:18:25 写道：

老大，我前两天给你发了老师的原题，邮件名称叫——老师调戏学生的统计题，收信人xie@yihui.name

请不要回帖不看帖！麻烦抬头看看

——原帖发布于 2010-05-23 09:24:23

http://cos.name/chinar/chinar-2010/

Originally posted on 2010-05-21 13:40:53