如何删除动物园系列的最后100个元素?
我知道名字[-element]符号,但我不能让它减去一个完整的部分
想象我有一个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绘制曲面?如果我的分数太多,那么这个矩阵就太大了.
干杯
我曾经read.table读过一个包含0.00001等数字的文件
当我把它们写回来时,write.table这些数字显示为1e-5
我怎样才能保留旧格式?
我知道您可以使用以下语法从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 这样的东西.
我有一个动物园时间序列,错过了几天.为了填补它并有一个连续的系列我做...
我从头到尾生成一个chron日期时间序列.
我把我的系列与这个合并.
我使用na.locf代替具有las遮挡的NAs.
我删除了syntetic chron序列.
我可以更容易吗?也许有一些与频率相关的指数函数?
我开始学习斯坦.
任何人都可以解释何时以及如何使用...等语法?
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) 假设我有这个简单的数据:
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书籍只提供一些例子,并且不鼓励使用饼图.
莱迪思更糟糕的是,你需要很多线才能让它变得吓人.
请问有人能帮助我获得一个漂亮而简单的饼图吗?例如......
没有任何R包能够轻松完成,没有20行代码吗?
是否可以在Spark中运行混合效果回归模型?(正如我们可以在R中使用lme4,或在Julia中使用MixedModels).
任何例子都会很棒.
我读过有一个GLMix函数,但我不知道用户是否可以直接使用它来拟合模型并获得系数和p值,或者它是否只能由机器学习库在内部使用.
我想转向Spark,因为我的数据集比内存大得多.
是否有任何其他通用数据库或框架能够执行类似磁盘流数据的操作?
我只看到一些能够进行简单的线性回归.
问候
我希望我的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,我该怎么办?