小编ska*_*kan的帖子

更改图表轴上刻度线的间距?

如何更改绘图轴上刻度线的间距?

我应该使用哪些参数与基础图或rgl

plot r

73
推荐指数
3
解决办法
19万
查看次数

R:删除向量的最后一个元素

如何删除动物园系列的最后100个元素?

我知道名字[-element]符号,但我不能让它减去一个完整的部分

r vector elements

43
推荐指数
5
解决办法
6万
查看次数

R:从x,y,z绘制3D表面

想象我有一个3列矩阵
x,y,z,其中z是x和y的函数.

我知道如何绘制这些点的"散点图" plot3d(x,y,z)

但是,如果我想要一个表面,我必须使用其他命令,如surface3d问题是它不接受与plot3d相同的输入它似乎需要一个矩阵与

(nº elements of z) = (n of elements of x) * (n of elements of x)
Run Code Online (Sandbox Code Playgroud)

我怎样才能得到这个矩阵?我尝试使用命令interp,就像我需要使用等高线图时一样.

如何在不计算此矩阵的情况下直接从x,y,z绘制曲面?如果我的分数太多,那么这个矩阵就太大了.

干杯

3d plot r matrix rgl

31
推荐指数
4
解决办法
10万
查看次数

数字格式,写入1e-5而不是0.00001

我曾经read.table读过一个包含0.00001等数字的文件

当我把它们写回来时,write.table这些数字显示为1e-5

我怎样才能保留旧格式?

format numbers r

25
推荐指数
3
解决办法
7万
查看次数

从GitHub更新所有包

我知道您可以使用以下语法从CRAN安装软件包:

install.packages(c("Rcpp"), dependencies=TRUE)
Run Code Online (Sandbox Code Playgroud)

您可以使用以下命令从CRAN更新所有这些内容:

update.packages()
Run Code Online (Sandbox Code Playgroud)

另一方面,您可以使用以下命令从GitHub安装软件包(编译它们):

install_github("hadley/tidyr")
Run Code Online (Sandbox Code Playgroud)

如何升级所有GitHub包?

我的意思是不一次重新安装(和编译)它们.像update.packages()github 这样的东西.

r github

23
推荐指数
3
解决办法
3238
查看次数

R:填写时间序列中缺少的日期?

我有一个动物园时间序列,错过了几天.为了填补它并有一个连续的系列我做...

我从头到尾生成一个chron日期时间序列.

我把我的系列与这个合并.

我使用na.locf代替具有las遮挡的NAs.

我删除了syntetic chron序列.

我可以更容易吗?也许有一些与频率相关的指数函数?

r time-series fill zoo

17
推荐指数
2
解决办法
1万
查看次数

斯坦.使用target + =语法

我开始学习斯坦.

任何人都可以解释何时以及如何使用...等语法?

target +=
Run Code Online (Sandbox Code Playgroud)

而不只是:

y ~ normal(mu, sigma)
Run Code Online (Sandbox Code Playgroud)

例如,在Stan手册中,您可以找到以下示例.

model {
real ps[K]; // temp for log component densities
sigma ~ cauchy(0, 2.5);
mu ~ normal(0, 10);
for (n in 1:N) {
for (k in 1:K) {
ps[k] = log(theta[k])
+ normal_lpdf(y[n] | mu[k], sigma[k]);
}
target += log_sum_exp(ps);
}
}
Run Code Online (Sandbox Code Playgroud)

我认为目标线增加了目标值,我认为这是后验密度的对数.
但是后验密度是什么参数?

什么时候更新并初始化?

在Stan完成(并收敛)之后,您如何获取其价值以及我如何使用它?

其他例子:

data {
  int<lower=0> J;          // number of schools 
  real y[J];               // estimated treatment effects
  real<lower=0> sigma[J];  // s.e. of effect estimates 
}
parameters …
Run Code Online (Sandbox Code Playgroud)

bayesian stan

14
推荐指数
2
解决办法
4458
查看次数

美丽的饼图与R

假设我有这个简单的数据:

 mydata <- data.frame(group=c("A", "B", "0", "AB"), FR=c(20, 32, 32, 16))
Run Code Online (Sandbox Code Playgroud)

如果我想从这个数据帧创建一个饼图,我可以这样做:

 with(mydata,pie(FR, labels=paste0(as.character(group), " ", FR, "%"), radius=1))
Run Code Online (Sandbox Code Playgroud)

基本馅饼

这很简单但可以接受.

我怎样才能得到与ggplot2或格子类似的东西?

经过多次试验和错误,我得到了

ggplot(mydata, aes(x = factor(1), y=FR,fill=factor(group)) ) + geom_bar(width = 1,stat="identity")+coord_polar(theta = "y") 
Run Code Online (Sandbox Code Playgroud)

ggplot

它更加复杂和丑陋.难道不是很容易吗?ggplot书籍只提供一些例子,并且不鼓励使用饼图.

莱迪思更糟糕的是,你需要很多线才能让它变得吓人.

请问有人能帮助我获得一个漂亮而简单的饼图吗?例如......

例1

例题

没有任何R包能够轻松完成,没有20行代码吗?

r ggplot2 lattice pie-chart

11
推荐指数
4
解决办法
4万
查看次数

Spark或其他技术中的混合效果模型

是否可以在Spark中运行混合效果回归模型?(正如我们可以在R中使用lme4,或在Julia中使用MixedModels).
任何例子都会很棒.

我读过有一个GLMix函数,但我不知道用户是否可以直接使用它来拟合模型并获得系数和p值,或者它是否只能由机器学习库在内部使用.

我想转向Spark,因为我的数据集比内存大得多.

是否有任何其他通用数据库或框架能够执行类似磁盘流数据的操作?
我只看到一些能够进行简单的线性回归.

问候

regression mixed-models apache-spark

10
推荐指数
1
解决办法
725
查看次数

R,Knitr,Rnw,美丽的科学数字

我希望我的knitr代码生成的所有数字看起来都不像一个老式的计算器.

在此输入图像描述

是否有任何选项可以获得最后一个数字(用10而不是e或E)?

options(scipen = ...)似乎没有那个选项.

我一直在搜索信息,我发现可以使用软件包siunitx直接在LaTex中完成,写下这样的每个数字\num {1e-10}

但是我想knitr自动为所有数字做过,包括表格中的数字.

PD:当我打印东西时,我该如何避免[1]?

PD2:也许是gsub的东西?

PD3:
我回到这个问题了.想象一下,我没有定义自己的表,但我从回归中得到它并使用xtable来生成它.

\documentclass{article}
\usepackage{siunitx}
\usepackage{booktabs}
\sisetup{     group-minimum-digits = {3},    group-separator = {,}, exponent-product = \cdot     }
\begin{document}

<<r, results='asis'>>=

library(xtable)
data(tli)
fm2 <- lm(tlimth ~ sex*ethnicty, data = tli)
xxx <- xtable(fm2)
print(xxx, booktabs = TRUE)

@
\end{document}
Run Code Online (Sandbox Code Playgroud)

但它效果不好.我应该使用哪些选项?

这是印刷的结果 在此输入图像描述

这是print +"booktabs = T"+我的函数beauty()的结果. 在此输入图像描述 问候.

我不知道为什么它产生两个表而不是1.而且数字没有正确对齐.无论如何,我不想依赖我的beauty()函数但只使用suintx,我该怎么办?

latex r scientific-notation digit knitr

9
推荐指数
2
解决办法
402
查看次数