小编ECI*_*CII的帖子

除了base和之外,如何删除所有已安装的软件包recommended？

我想将本地目录同步到Google云端存储中的存储分区.我想复制远程不存在的本地文件,跳过已经存在的远程和本地文件.使用GSUtil可以做到这一点吗？我似乎无法为GSUtil找到"同步"选项或"不要覆盖".可以编写这个脚本吗？

我在Linux上(Ubuntu 12.04)？

sudo apt install awscli安装当前版本 1.16.113-1

官方网站（https://docs.aws.amazon.com/cli/latest/userguide/install-cliv2-linux.html）支持aarch64但不支持 armhf。还有其他方法或我缺少的东西吗？v1 和 v2 之间的差异不是很大，但例如 v1 不支持--storage-class DEEP_ARCHIVE

我有

plot(rnorm(120), rnorm(120), col="darkblue", pch=16, xlim=c(-3,3), ylim=c(-4,4))
points(rnorm(120,-1,1), rnorm(120,2,1), col="darkred", pch=16)
points(c(-1,-1.5,-3), c(4,2,0), pch=3, cex=3)

我想通过绘制一条通过一组点的平滑曲线来描绘图形的一部分.我可以定义3-4组点但我无法定义一个函数.我希望在R(而不是GIMP)中这样做,因为我想提交为SVG.我想要实现的目标如下

这可能吗？我知道这不是一个复杂的图形问题,但任何基础R解决方案都可以.

我有2个数据帧,每个数据帧的列数不同.一些列在2个数据帧之间是通用的.如何才能将两个数据帧的公共列重新绑定到新的数据帧？

我试过,library(plyr);rbind.fill(A,B)但它在不匹配的列中设置NA值,这对我没有帮助.

非常感谢EC

该muhaz软件包使用内核平滑方法从正确的删失数据中估计危险函数.我的问题是,有没有办法获得计算的危险函数的置信区间？muhaz

options(scipen=999)
library(muhaz)
data(ovarian, package="survival")
attach(ovarian)
fit1 <- muhaz(futime, fustat)
plot(fit1, lwd=3, ylim=c(0,0.002))

在上面的例子中muhaz.object fit有一些条目fit1$msemin,fit1$var.min,fit1$haz.est然而,它们的长度是一半fit1$haz.est.

有没有想法可以提取危险函数的置信区间？

编辑:我根据@ user20650的建议尝试了以下内容

options(scipen=999)
library(muhaz)
data(ovarian, package="survival")
fit1 <- muhaz(ovarian$futime, ovarian$fustat,min.time=0, max.time=744)


h.df<-data.frame(est=fit1$est.grid, h.orig=fit1$haz.est)

for (i in 1:10000){
d.s.onarian<-ovarian[sample(1:nrow(ovarian), nrow(ovarian), replace = T),]
d.s.muhaz<-muhaz(d.s.onarian$futime, d.s.onarian$fustat, min.time=0, max.time=744 )
h.df<-cbind(h.df, d.s.muhaz$haz.est)
}


h.df$upper.ci<-apply(h.df[,c(-1,-2)], 1,  FUN=function(x) quantile(x, probs = 0.975))
h.df$lower.ci<-apply(h.df[,c(-1,-2)], 1,  FUN=function(x) quantile(x, probs = 0.025))
plot(h.df$est, h.df$h.orig, type="l", …

我有以下查找表:

lkp <- data.frame(
         x=c(0,0.2,0.65,0.658,1.3,1.76,2.7), 
         y=c(1,1,1,0.942,0.942, 0.92, 0.89)
       )

我想得到给定X值的Y值.

如果表中存在X值,则应返回表的确切Y. 如果X值不存在,那么Y值应该作为2个最近邻居(仅2个最近邻居)的线性插值返回.我不想让模型适合整体数据.

对于上表

for X=0.2 Y=1 (exact lookup) 
for X=2 Y=0.91 (linear interpolation between the last 2 rows of the data frame)

有没有准备好的功能呢？

我有情节1

curve(exp(x), from=1, to=5, lwd=5)
curve(150-exp(x), from=1, to=5, lwd=5, col="darkblue",add=T)

在里面我想补充下面的情节2

par(mar=c(7,7,1,1))
curve(exp(x), from=1, to=5, lwd=7, xlab="chi", ylab="exp(x)", cex.lab=4,axes=F)
axis(1, labels=NA,at=c(0,5))
axis(2, labels=NA,at=c(0,150))
text(1,120,"Alpha",adj=c(0,0),cex=3)
text(3.5,10,"Beta",adj=c(0,0),cex=3)

为了获得以下内容

我还想让情节2变得透明,这样如果在情节2后面有一些情节1的元素,它们仍会显示(就像蓝线一样).同样重要的是图2的较大标签以及其轴上没有标签和刻度.

这可能吗？请仅基于R解决方案(没有ggplot2 /无格)

例如,我生成了一些高分辨率的出版质量图

library(plot3D)
Volcano<-volcano
zf=10 #zoom factor
tiff("Volcano.tif", width=1800*zf, height=900*zf, res=175*zf, compression="lzw")
image2D(z = Volcano, clab = "height, m",colkey = list(dist = -0.20, shift = 0.15,side = 3, length = 0.5, width = 0.5,cex.clab = 1.2, col.clab = "white", line.clab = 2,col.axis = "white", col.ticks = "white", cex.axis = 0.8))
dev.off()

该文件是22 MB.

现在我用GIMP打开文件而不做任何其他事情我将它导出为"Volcano gimp.tif"(不要改变分辨率,或做任何其他事情).GIMP生成一个1.9 MB的文件("Volcano gimp.tif").

imagemagick 报告相似图片统计:

$ identify Volcano.tif
Volcano.tif TIFF 18000x9000 18000x9000+0+0 8-bit DirectClass 22.37MB 0.000u 0:00.000
$ identify "Volcano gimp.tif"
Volcano gimp.tif TIFF 18000x9000 …

我想在两组上显示倾向得分匹配统计数据(BW中无法比拟,颜色匹配),并希望使用像下面这样的镜像直方图

是否可以在基础R中叠加4个不同的直方图？有没有提供此功能的包？

x轴是[0,1]有界的(它是一个概率),BW列总是大于或等于彩色列(即彩色列后面不能有BW列).

图片来自http://www.ncbi.nlm.nih.gov/pubmed/22244556