小编Yan*_*ang的帖子

我想在图例中写一个像“R^2=0.00575”这样的方程，数字 0.00575 可以自动嵌入到图例中。这是一个例子。

set.seed(100)
x=rnorm(100)
y=1:100
fit=lm(y~x)
R_squared=format(summary(fit)$r.squared,digits = 3)
plot(x,y,type="l")
legend("topleft",legend =expression(R^{2}~"="~R_squared),bty = "n")

如图所示，变量"R_squared"未嵌入方程中。有什么解决办法吗？谢谢。

我正在尝试使用 R 包构建不同的回归模型caret。对于数据来说，它既包括数值，也包括因子。

问题 1： 在回归模型中同时包含数值和因子的正确方法是什么caret？

问题2：回归模型通常需要进行数据预处理（中心和尺度），那么因子的预处理如何进行？

library(caret)
data("mtcars")
mydata = mtcars[, -c(8,9)]
set.seed(100)
mydata$dir = sample(x=c("N", "E", "S", "W"), size = 32, replace = T)
mydata$dir = as.factor(mydata$dir)
class(mydata$dir)    # Factor with four levels

MyControl = trainControl(
  method = "repeatedcv", 
  number = 5, 
  repeats = 2, 
  verboseIter = TRUE, 
  savePredictions = "final"
)

model_glm <- train(
  hp ~ ., 
  data = mydata, 
  method = "glm", 
  metric = "RMSE", 
  preProcess = c('center', 'scale'), 
  trControl …

我现在正试图绘制加拿大的地图,ggplot2我发现经度和纬度没有在地图中正确显示.有什么解决方案吗？非常感谢.

arcgis shapfile是从https://www.arcgis.com/home/item.html?id=dcbcdf86939548af81efbd2d732336db下载的

library(ggplot2)
library(rgdal)
countries<-readOGR("Canada.shp", layer="Canada")
ggplot()+geom_polygon(data=countries,aes(x=long,y=lat,group=group),fill='white',color = "black")

地图中的经度应该是110W,100W,90W.并且地图中的纬度应该像50N,60N,70N.但是,现在情况并非如此.

我正在使用ecdf()函数从一些随机样本中计算经验累积密度函数（ECDF）：

set.seed(0)
X = rnorm(100)
P = ecdf(X)

现在P给出ECDF，我们可以绘制它：

plot(P)
abline(h = 0.6, lty = 3)

我的问题是：如何找到样本值x，例如P(x) = 0.6ECDF的0.6位数，或ECDF与的交点的x坐标h = 0.6？

对于R中的一般情节,legend用于对数字进行编号.

set.seed(100)
Mydata=rnorm(65)
Year=1950:2014
plot(x=Year,y=Mydata,type = "l")
legend("topleft","(a)",bty = "n")

我想知道我们如何使用ggplot2做同样的事情.谢谢.

我有一些数据,我想找到连续值<= 2,连续值的长度应> 3.我的目标是(1)找到不同连续值组的长度,以及(2)找到每组中第一个值的位置.我试过以下代码:

set.seed(100)
pre = sample(x=1:5, size = 90, replace = T)
which(pre<=2)

这将产生以下结果:

[1]  1  2  4  8 10 13 14 17 18 19 26 30 33 37 40 41 49 50 51 52 53 54 56 57 58 60 66 69 72 80 85 88 89

因此,保守的价值观群体包括:(1)17,18,19; (2)40,41; (3)49,59,51,52,53,54; (4)56,57,58; (5)88,89.

但是,由于我只需要长度> = 3的连续值,因此应从结果中排除组(2)和(5).我想知道我怎么能在R中做到这一点？谢谢你的帮助.

我有一个小标题，我正在尝试计算多个指标。

library(tidymodels)
price = 1:50
prediction = price * 0.9
My_tibble = tibble(price=price, prediction=prediction)

# The following code can calculate the rmse
My_tibble %>% 
  rmse(truth = price, estimate = prediction)

# Is it possible to calculate `rmse` and `rsq` at the same time?
# The following code reports an error: object '.pred' not found
My_tibble %>%
  rmse(truth = price, estimate = prediction ) %>% 
  rsq(truth = price, estimate = prediction )

把问题延伸一点，是否可以同时计算rmse和？cor

My_tibble %>%
      rmse(truth = price, …

我是Python新手，我对以下代码感到很困惑：

示例1：

n = 1
m = n
n = 2
print(n)   # 2
print(m)   # 1

示例2：

names = ["a", "b", "c"]
visitor = names
names.pop()
print(names)      # ['a', 'b']
print(visitor)    # ['a', 'b']

示例 1 显示n是 1 并且m是 1。然而，示例 2 显示names是 ['a', 'b']，并且visitor也是 ['a', 'b']。

对我来说，示例 1 和示例 2 很相似，所以我想知道为什么结果如此不同？谢谢。

我有一些数据，并且尝试了一个filled.contour看起来不错的情节。然而，图例很难控制，所以我正在考虑使用ggplo2。但我不知道如何绘制filled.contourusing ggplot2。

数据包含 840 行（代表日期）和 12 列（代表 12 个时间尺度）。这是一个例子

set.seed(66)
Mydata <- sample(x=(-3:3),size = 840*12,replace = T)
Mydata <- matrix(data=Mydata,nrow=840,ncol=12)
Dates <- seq(from=1948+1/24, to= 2018,by=1/12)
data.breaks <- c(-3.5,-2.5,-1.5,0,1.5,2.5,3.5)
filled.contour(Dates,seq(1:12),Mydata,col=cols(11),xlab="",ylab="time-scale",levels=data.breaks)

正如我们所看到的，图例间隔不是我想要的。我想-3.5,-2.5,-1.5,0,1.5,2.5,3.5在图例上展示，我相信用 来做到这一点要容易得多ggplot2。谢谢你的帮助。

我有两个概率密度函数,我想知道它们的分布是否相似.我知道R中的KS测试可以做到这一点,但是当我运行代码时,会发生错误.谢谢你的帮助.

set.seed(100)
a=density(sample(x=1:30,size = 30,replace = T))
b=density(sample(x=1:40,size = 35,replace = T))
plot(a)
lines(b)

ks.test(a,b)
Error in ks.test(a, b) : 
 'y' must be numeric or a function or a string naming a valid function

我试图在R中创建一个箱线图，但是，我发现该图的y轴刻度值错误。

.rdata可从https://www.dropbox.com/s/vbgf3mhgd2mjx8o/Mydata2.rdata?dl=0获得

load("Mydata2.rdata",.GlobalEnv)
boxplot(Value~Type+Level, data=Mydata2)

如图所示，y轴标记为“ 0、50、100”，但是我的数据范围为-36.9到133.7。我想知道如何解决这个问题？