小编Kar*_*ius的帖子

我想在R中叠加两个散点图,使得每组点具有其自己的(不同的)y轴(即,在图中的位置2和4),但是这些点看起来叠加在同一图上.

有可能这样做plot吗？

编辑显示问题的示例代码

# example code for SO question
y1 <- rnorm(10, 100, 20)
y2 <- rnorm(10, 1, 1)
x <- 1:10
# in this plot y2 is plotted on what is clearly an inappropriate scale
plot(y1 ~ x, ylim = c(-1, 150))
points(y2 ~ x, pch = 2)

我想用R覆盖同一设备上的2个密度图.我该怎么做？我搜索了网络,但我没有找到任何明显的解决方案(我对R来说很新).

我的想法是从文本文件(列)中读取数据然后使用

plot(density(MyData$Column1))
plot(density(MyData$Column2), add=T)

这种精神......

提前致谢

我注意到三星S8,Android 7.0上的问题(更新.这也发生在Android 7.0上:三星S7,Nexus 5x),它告诉(经过几次测试)应用程序扫描过于频繁:

08-14 12:44:20.693 25329-25329/com.my.app D/BluetoothAdapter: startLeScan(): null
08-14 12:44:20.695 25329-25329/com.my.app D/BluetoothAdapter: STATE_ON
08-14 12:44:20.696 25329-25329/com.my.app D/BluetoothAdapter: STATE_ON
08-14 12:44:20.698 25329-25329/com.my.app D/BluetoothLeScanner: Start Scan
08-14 12:44:20.699 25329-25329/com.my.app D/BluetoothAdapter: STATE_ON
08-14 12:44:20.700 25329-25329/com.my.app D/BluetoothAdapter: STATE_ON
08-14 12:44:20.700 25329-25329/com.my.app D/BluetoothAdapter: STATE_ON
08-14 12:44:20.701 25329-25329/com.my.app D/BluetoothAdapter: STATE_ON
08-14 12:44:20.703 4079-4093/? D/BtGatt.GattService: registerClient() - UUID=dbaafee1-caf1-4482-9025-b712f000eeab
08-14 12:44:20.807 4079-4204/? D/BtGatt.GattService: onClientRegistered() - UUID=dbaafee1-caf1-4482-9025-b712f000eeab, clientIf=5, status=0
08-14 12:44:20.808 25329-25342/com.my.app D/BluetoothLeScanner: onClientRegistered() - status=0 clientIf=5 mClientIf=0
08-14 12:44:20.809 4079-7185/? D/BtGatt.GattService: start scan with …

我经常遇到需要检查获得的差异是否高于机器精度的情况。似乎为此目的，R 有一个方便的变量：.Machine$double.eps. 但是，当我转向 R 源代码以获取有关使用此值的指南时，我看到了多种不同的模式。

例子

以下是stats库中的一些示例：

t.test.R

if(stderr < 10 *.Machine$double.eps * abs(mx))

chisq.test.R

if(abs(sum(p)-1) > sqrt(.Machine$double.eps))

积分

rel.tol < max(50*.Machine$double.eps, 0.5e-28)

影响力

e[abs(e) < 100 * .Machine$double.eps * median(abs(e))] <- 0

原理图

if (any(ev[neg] < - 9 * .Machine$double.eps * ev[1L]))

等等。

问题

1. 一个如何理解这些不同背后的理由10 *，100 *，50 *和sqrt()修饰？
2. 是否有关于.Machine$double.eps用于调整因精度问题而导致的差异的指南？

这个R代码抛出警告

# Fit regression model to each cluster
y <- list() 
length(y) <- k
vars <- list() 
length(vars) <- k
f <- list()
length(f) <- k

for (i in 1:k) {
  vars[[i]] <- names(corc[[i]][corc[[i]]!= "1"])
  f[[i]]  <- as.formula(paste("Death ~", paste(vars[[i]], collapse= "+")))
  y[[i]]  <- lm(f[[i]], data=C1[[i]]) #training set
  C1[[i]] <- cbind(C1[[i]], fitted(y[[i]]))
  C2[[i]] <- cbind(C2[[i]], predict(y[[i]], C2[[i]])) #test set
}

我有一个训练数据集(C1)和一个测试数据集(C2).每个都有129个变量.我做了k意味着对C1进行聚类分析,然后基于聚类成员分割我的数据集,并创建了不同聚类的列表(C1 [[1]],C1 [[2]],...,C1 [[k] ]).我还为C2中的每个案例分配了一个集群成员资格,并创建了C2 [[1]],...,C2 [[k]].然后我对C1中的每个簇进行线性回归.我的因变量是"死亡".我的预测变量在每个群集中都不同,变量[[i]](i = 1,...,k)显示了预测变量名称列表.我想为测试数据集中的每个案例预测死亡(C2 [[1]],...,C2 [[k]).当我运行以下代码时,对于某些集群.

我收到了这个警告:

In predict.lm(y[[i]], C2[[i]]) :
prediction from a rank-deficient …

我有以下内容,以条形图数据框.

c1 <- c(10, 20, 40)
c2 <- c(3, 5, 7)
c3 <- c(1, 1, 1)
df <- data.frame(c1, c2, c3)
ggplot(data=df, aes(x=c1+c2/2, y=c3)) +
  geom_bar(stat="identity", width=c2) +
  scale_fill_manual(values=c("#FF6666"))

我最终只有灰色条纹:条形图的灰色条纹

我想改变酒吧的颜色.我已经尝试过来自http://www.cookbook-r.com/Graphs/Colors_(ggplot2)/的不同scale_fill_manual, 但仍然有灰色条.

谢谢您的帮助.

我正在尝试生成一个条形图,但是y轴刻度太短.这是我的代码:

barplot(as.matrix(dat), log="y", ylim=c(10000,100000000), beside=TRUE,
        ylab = "Number of reads", col = c("gray","black","white"))

它离开了轴的空间(根据ylim),但没有填充实际的轴.我经历过?barplot并尝试了一些事情(从谷歌搜索我认为xpd = F, yaxs = c(10000,10000000,5)应该工作,但它没有).

我知道这是一件小事,但这正是我长期坚持的问题,而不是实际工作,所以任何帮助都会非常感激!

编辑:为输入人员干杯!

我最初没有使用ylim进行绘图,但最终会出现一个更奇怪的轴(同样的问题); 我实际上选择了我的ylim值来给它一个更好的间隔轴.

这是数据:

dat <- read.table(text="D2,D3n,D3m,D4n,D4m
21234722,34262282,31920464,25486357,20712943
35343,64403,22537,39934,46547
126646,312286,101105,81537,76944", header=TRUE, sep=",")

编辑2:@DWin没错 - 我更新了我的R,现在情节很好 - 谢谢大家!

可以帮我改变这些图的标题文字大小.即使它们变大？

脚本

ggplot(NMPSCMOR, aes(Length,  fill=Year)) + 
  geom_histogram(position="dodge", binwidth=60,  colour="black") + xlim(0, 600) +
  scale_fill_grey(start = 1, end = 0)+ 


  geom_vline(data=ddply(NMPSCMOR, Year~Morphology~Sector2, numcolwise(mean)), 
             mapping=aes(xintercept=Length,color=Year),  linetype=2, size=1) + 
  scale_color_grey(start=1,end=0)+

  xlab("Length Class") +
  ylab(expression(paste("Total Count"))) + #( ", m^2, ")", sep = 
  facet_wrap( ~ Morphology + Sector2, ncol=3, scales = "free") +
  theme(

    panel.grid.minor = element_blank(),               #removes minor grid lines
  panel.grid.major = element_blank())

我正在用corrplot绘制相关图.我想绘制相关系数:

require(corrplot)
test <- matrix(data = rnorm(400), nrow=20, ncol=20)
corrplot(cor(test), method = "color", addCoef.col="grey", order = "AOE")

但它们在情节中太大了:

有没有办法让coefficent的字体更小？我一直在看,?corrplot但只有参数可以更改图例和轴字体大小(cl.cex和tl.cex).pch.cex也不起作用.

在图中,是否可以略微抖动状态缩写标签,以便它们不重叠？如果我使用check_overlap = TRUE,那么它会删除一些重叠的观察,我不希望这样.我也不想要geom_label_repel,因为它有标签伸出并穿过我所包含的45度线(我不想发生)

这是我的代码的相关部分供参考:

ggplot(df, aes(x = huff_margin_dem, y = margin16dem_state, label = abbrev)) +
  geom_abline(intercept = 0) +
  geom_text(fontface = "bold")