科学网—leymus的博文

用R进行聚类分析时R2的简便计算方法

牛海山 2015-2-20 15:00

R在进行聚类分析时，没提供R^2。网上有人提供了自己编制的代码。觉得还应该有更简洁的计算方法,我于是借鉴了kmeans函数的帮助，写了下面三行代码计算R2。使用了scale() 和 by() 两个函数，大大减少了代码量。这两个函数都在base包里，不需要额外加载包。使用时需要注意一点：只保留分类变量列，把名字列去掉。 tss ...

个人分类: R使用|6693 次阅读|没有评论

超级计算机上用PBS提交R计算任务详解

热度 1 牛海山 2014-12-22 22:07

最近想用R做点大活儿，可能得用超级计算机。把我能接触的超级计算机研究了一下，试着用PBS提交了几次计算任务。中间产生了一些小问题，都一一解决了。现把心得共享。下面是依葫芦画瓢写的一个用于提交计算任务的脚本。PBS是Altair公司出品的一个任务管理系统，好像超级计算机上用得很普遍。比如下面这个脚本的名字叫“Rjo ...

个人分类: R使用|12331 次阅读|1 个评论热度 1

R做图时无法转换unicode字符的临时解决方案

热度 1 牛海山 2014-12-2 20:52

不知道什么原因，一些情况下unicode字符在R出图时会转换失败。报错信息诸如： 'mbcsToSbcs'里转换'④'出错：91代替了dot 或者 conversion failure on '≤' in 'mbcsToSbcs' : dot substituted for e2 ...

个人分类: R使用|8190 次阅读|3 个评论热度 1

用R画中国地图时如何把南海做成小图

热度 2 牛海山 2014-12-1 14:08

平常的中国地图，南海部分是一个独立的小图，放在整张地图的一角的。不知道专业的地图软件是不是很容易实现这样的效果，反正我是费了整整一个工作日的时间才在R里做到这一点。赶紧上来宣传一下，一个是表表功，再一个是记录下来、留作后用。其实，我以前想都没想过为什么要把中国地图折成这样的两截。也许审美观随着 ...

个人分类: R使用|8609 次阅读|3 个评论热度 2

能把R下多图合并

热度 1 牛海山 2014-11-3 20:01

在R下多图合并，通常用par()或者layout()。无论用哪一个命令，在微调时候，只有对plot的内核有深入理解，才能做好。能把这内核讲深、讲透的文章不多，是使用者大为苦恼之事。下面提供两个我认为很好的网文。只有吃透了，才能讲清楚。对参数oma（边际外区）与mai（边际区）的意义，是首先需要了解的；可以参考这篇网文 ...

个人分类: R使用|4842 次阅读|3 个评论热度 1

R中如何去掉/删除字符串中指定长度的字符

牛海山 2012-10-18 17:21

如果想从左侧删除4个字符用sub('^....','',str) 如果想从右侧删除4个字符用sub('....$','',str) 例如: str-"I have a dream" sub('^....','',str) "ve a dream" sub('....$','',str) "I have a d"

个人分类: R使用|25185 次阅读|没有评论

R graph：如何自定义坐标轴刻度标示(tick label)

牛海山 2012-6-14 10:23

这里要说的不是轴标。轴标用xlab或者ylab赋值。我要说的是轴上的刻度标志（tickmark）相对应的那些标识（labels)如何自定义。首先要在plot里面压制自动生成坐标轴；因为你想自定义，当然需要压制自动生成的了。例如用xaxt="n"就压制了x轴刻度标志以及刻度标识。其次使用axis()命令自定义刻度标识。重要参数有四个：side( ...

个人分类: R使用|18158 次阅读|没有评论

R数据分析中interaction.plot很有用处

牛海山 2012-3-1 00:58

可以很方便地展示各个因素之间相互作用关系。如果ANOVA显示某些互作项显著，可以用这个命令迅速查看互作的形式。关键是方便。最基本语法是：interaction.plot(x.factor,trace.factor,response,...) 其中，x.factor是形成x轴的因素（互作的一个构成因素）；trace.factor是互作的另一个构成因素，它的每个水平将会形成互作图 ...

个人分类: R使用|17372 次阅读|没有评论

R做图中插入特殊符号

热度 3 牛海山 2012-2-20 14:59

R做图中经常需要输入一些特殊符号（例如希腊字母），如果这些符号无法手工输入，那么可以以"uXXXX"代码形式输入。例如可以用“u0394”输入希腊字母DELTA，用“u2030”输入千分号。更多符号可上网查unicode字符集。

个人分类: R使用|16274 次阅读|4 个评论热度 3

帐号		自动登录	找回密码
密码			注册

关闭 安全验证

leymus

关闭安全验证