接受各种R语言难题

cbw1243

42655

504

收藏 2016-06-26

暑假漫长，生活没有挑战。在此征集各类R语言疑难杂症，欢迎各位骚扰。
提出给力问题的，将予以100论坛币奖励。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

jiangqing001

2016-6-26 10:37:09

霸气啊~~~

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

万人往LVR

2016-6-26 11:39:12

想用R调用qq给某个好友发消息，请问楼主怎么搞

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

yangkongsong

2016-6-26 12:26:21

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243

2016-6-26 22:25:32

万人往LVR 发表于 2016-6-26 11:39
想用R调用qq给某个好友发消息，请问楼主怎么搞

可以利用 system 指令。比如： system('open -a QQ') 这个适用于mac
open 是cmd指令，a 是application， QQ就是你想要玩的。

不过coding的时候最好不要开QQ，容易走神写出BUG。😄

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

elickychan

2016-6-26 22:50:03

cbw1243 发表于 2016-6-26 09:15
暑假漫长，生活没有挑战。在此征集各类R语言疑难杂症，欢迎各位骚扰。
提出给力问题的，将予以100论坛币奖 ...

想用R做空间面板分位数和空间面板门限求大神赐教

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

ivyvivi

2016-6-26 22:53:26

R，有点统计方面的问题，可以问么？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

elickychan

2016-6-26 22:54:10

cbw1243 发表于 2016-6-26 09:15
暑假漫长，生活没有挑战。在此征集各类R语言疑难杂症，欢迎各位骚扰。
提出给力问题的，将予以100论坛币奖 ...

还有就是R怎么做空间面板联立方程模型分析实现GS3SLS估计？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

太阳之子追梦者

2016-6-26 23:02:05

大神，果然任性！！！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

MILK$很好喝

2016-6-26 23:31:25

cbw1243 发表于 2016-6-26 09:15
暑假漫长，生活没有挑战。在此征集各类R语言疑难杂症，欢迎各位骚扰。
提出给力问题的，将予以100论坛币奖 ...

楼主能写一下关于词云怎么做的吗，没找到，感谢

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243

2016-6-27 00:38:34

elickychan 发表于 2016-6-26 22:54
还有就是R怎么做空间面板联立方程模型分析实现GS3SLS估计？

一类模型里可能出现很多具体的问题，能不能把你的问题具体一些。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243

2016-6-27 00:41:11

MILK$很好喝发表于 2016-6-26 23:31
楼主能写一下关于词云怎么做的吗，没找到，感谢

我没有做过，但是推荐你看一下 “wordcloud” 这个包。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243

2016-6-27 00:41:53

太阳之子追梦者发表于 2016-6-26 23:02
大神，果然任性！！！

大神就免了，只为共同进步。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243

2016-6-27 00:42:28

ivyvivi 发表于 2016-6-26 22:53
R，有点统计方面的问题，可以问么？

Very welcome!

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

nkunku

2016-6-27 02:09:45

您好！R做实验设计，都有哪些好用的包或者书籍，可否推荐一下？谢谢！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243

2016-6-27 02:32:44

nkunku 发表于 2016-6-27 02:09
您好！R做实验设计，都有哪些好用的包或者书籍，可否推荐一下？谢谢！

我个人比较喜欢 <The Art of R Programming> 这本书，作者是Norman Matloff.
至于package，没有一个绝对好用的，主要看你要处理的问题。
我自己也写过包。刚发表出来觉得很好用，一个月以后觉得好垃圾。。。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

sophy499

2016-6-27 07:04:26

用hclust的时候只能用distance matrix，有没有办法用distance做cluster但是在y轴显示similarity也就是1-distance？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

pkl

2016-6-27 10:03:20

牛人，顶一下，向你学习！大家一起共同努力、共同进步！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243

2016-6-27 10:06:44

sophy499 发表于 2016-6-27 07:04
用hclust的时候只能用distance matrix，有没有办法用distance做cluster但是在y轴显示similarity也就是1-dis ...

我尝试用mtcars数据，是显示distance的。
plot(hclust(dist(mtcars)), hang = -1)
还可以对坐标轴进行customize，比如：
axis(side = 2, at = seq(0, 500, 100), col = "#F38630", labels = FALSE, lwd = 2)

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

太阳之子追梦者

2016-6-27 11:16:46

cbw1243 发表于 2016-6-27 00:41
大神就免了，只为共同进步。

刚好，这几个月，我都专注在R领域的爬虫、文本挖掘、算法

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

MILK$很好喝

2016-6-27 16:02:13

cbw1243 发表于 2016-6-27 00:41
我没有做过，但是推荐你看一下 “wordcloud” 这个包。

嗯嗯，谢谢

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ivyvivi

2016-6-27 16:50:12

大神，我想用GLM模型来拟合几个变量的关系，Y~X1+X2+X3+X4+X5. 同伴叫我用negative binomial，因为研究的东西是疾病或者是社交传播，查了一下，好像一般都是用negative binomial，不过看解释都是说传播数量符合这个分布。
而我这里边Y 是个efficiency，是传播的效率，取值范围是[0,1], 而自变量，除了X1是整数，其他的也都是[0,1]区间的。
我这里到底该用什么样的线性或者非线性模型来拟合呢？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243

2016-6-27 22:33:01

ivyvivi 发表于 2016-6-27 16:50
大神，我想用GLM模型来拟合几个变量的关系，Y~X1+X2+X3+X4+X5. 同伴叫我用negative binomial，因为研究的东 ...

首先，GLM拟合是很简单的。以下是我random的数据和GLM回归
y <- runif(100)
x1 <- rpois(100,2)
x2 <- rnorm(100)
test <- MASS::glm.nb(y ~ x1 + x2) （会出现错误，因为y是continuous。但是你可以忽略，因为得到的参数估计是无偏的）
summary(test)
其次，在任何不用OLS的情况下，都要认真思考一下为什么不用OLS的原因。你已经做了功课，literature里说传播数量符合negative binomial分布。同时，你也意识到你可能不能依赖这个，因为你研究的变量不同。negative binomial的特质就是分布的方差大于mean，俗称over-dispersion. 你的模型里的Y的条件分布未必是这个。效率的取值应该是0到1，除非有特殊的原因导致你有很多的0效率观测值。不然，你是没有理由选择glm的。
最后，如果你的样本不够大，另外over-dispersion并不明显，你可以选择quasi－poisson。当然，从模型选择角度，你也可以用AIC来比较各种模型。效率评估里面有个DEA模型挺好用的，有时间可以注意一下。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

fcfc2013

2016-6-27 23:30:30

library(xml2)
library(rvest)

fun <- function(i,url){
  url <- paste0(url,i)
  web <- read_html(url)
  # 小区名称
  id <- web %>% html_nodes("div.col-1 h2 a") %>% html_text()
  # 每平米均价
  price <- web %>% html_nodes("div.average span.num") %>% html_text() %>% as.numeric()
  # 小区地址
  where <- web %>% html_nodes("div.where span.region") %>% html_text()
  # 房屋面积
  area <- web %>% html_nodes("div.area span") %>% html_text()
  lianjia <- data.frame(id, price, where, area)
}

lianjia <- data.frame()
url <- ("http://wh.fang.lianjia.com/loupan/pg")

for (i in 2:47){
  fun(i,url)
  lianjia <- rbind(lianjia,fun(i,url))
  }
我用 R的rvest包函数爬链家网房产价格，循环采集制定的数据，想循环爬取47个页面，（每个页面有10条数据）
写了个fun（自定义爬取内容函数）  用for循环 47次发现每次运行就采集第2页10条数据；  这个怎么回事？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243

2016-6-28 01:04:51

fcfc2013 发表于 2016-6-27 23:30
library(xml2)
library(rvest)

我刚试运行了一下。你的代码没有问题。
这样做你会发现原因：
i <- 1
然后把fun里面的程序一个一个的运行。你会发现 area这个vector的length只有9。其他的都是10. 所以不能合并成dataframe，
你可以继续检查一下这是为什么。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

lww1993

2016-6-28 19:00:08

大牛，如何读一篇中文文档，然后按照不同的字不同的标点去分类统计频数，然后写到EXCEL表格中？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

fcfc2013

2016-6-28 22:09:26

cbw1243 发表于 2016-6-28 01:04
我刚试运行了一下。你的代码没有问题。
这样做你会发现原因：
i

弄了半天还是没搞定“area”这个列为什么是“9”；我去掉这个。
for循环可以正常一次遍历完47页面；
R提倡隐式循环，如用“apply” 如何写，达到 “for” 同样的效果了？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243

2016-6-28 22:52:17

fcfc2013 发表于 2016-6-28 22:09
弄了半天还是没搞定“area”这个列为什么是“9”；我去掉这个。
for循环可以正常一次遍历完47页面；
R ...

apply并不能增加运转速度。你可以用lapply

result <- lapply(1:47, function(i) fun(i, url) #return a list
result <- do.call(rbind, result) #combine all results

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243

2016-6-28 22:58:13

lww1993 发表于 2016-6-28 19:00
大牛，如何读一篇中文文档，然后按照不同的字不同的标点去分类统计频数，然后写到EXCEL表格中？

R的文本抓取是很方便的，比如grep， grepl， strsplit
写到excel就更简单了，比如 write.csv。文本读取可以用 read.table
当然，具体怎么做看你具体的问题

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

太阳之子追梦者

2016-6-29 12:30:01

字符串处理：
[[9995]]
[1] "送货快" "物流" "就是" "好"
[[9996]]
[1] "不错" "值得" "推荐"
[[9997]]
[1] "价格便宜" "味道不错"

如何转为变为
[9995]  "送货快" "物流" "就是" "好"
[9996]  "不错" "值得" "推荐"
[9997]  "价格便宜" "味道不错"

然后按行导出

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群