小编Mat*_*ath的帖子

我有这个向量:

x = c(1,1,1,1,1,0,1,0,0,0,1,1)

而且我只想为正数做累积求和.我应该有以下向量作为回报:

xc = (1,2,3,4,5,0,1,0,0,0,1,2)

我怎么能这样做？

我已经尝试过了,cumsum(x)但是它会对所有值进行累积求和,并给出:

cumsum(x)
[1] 1 2 3 4 5 5 6 6 6 6 7 8

我有这么简单 data.frame

 lat<-c(1,2,3,10,11,12,20,21,22,23)
 lon<-c(5,6,7,30,31,32,50,51,52,53)
 data=data.frame(lat,lon)

想法是基于距离找到空间聚类

首先,我绘制地图(lon,lat):

plot(data$lon,data$lat)

很明显我有三个基于点位置之间距离的聚类.

为此目的,我在R中尝试了这个代码:

d= as.matrix(dist(cbind(data$lon,data$lat))) #Creat distance matrix
d=ifelse(d<5,d,0) #keep only distance < 5
d=as.dist(d)
hc<-hclust(d) # hierarchical clustering
plot(hc)
data$clust <- cutree(hc,k=3) # cut the dendrogram to generate 3 clusters

这给出了:

现在我尝试使用簇中的颜色绘制相同的点

plot(data$x,data$y, col=c("red","blue","green")[data$clust],pch=19)

结果如下

这不是我想要的.

实际上,我想找到像这样的情节

谢谢你的帮助.

我有

X=c(20 ,18, 34, 45, 30, 51, 63, 52, 29, 36, 27, 24)

随着boxplot,我试图绘制quantile(X,0.25),quantile(X,0.75) 但这不是真正的R中的boxplot中的上下四分位数

boxplot(X)
abline(h=quantile(X,0.25),col="red",lty=2)
abline(h=quantile(X,0.75),col="red",lty=2)

你知道为什么吗？

为什么 pandas 有两个 Boxplot 函数：pandas.DataFrame.plot.box()和pandas.DataFrame.boxplot()？

df = pd.DataFrame(np.random.rand(10, 5), columns=['A', 'B', 'C', 'D', 'E'])
df.plot.box()

df.boxplot()

我想知道是否有一个函数可以随机添加或减法运算符:

x +- y

我想生成一个序列,如:

c(rep(1,7),rep(2,7),rep(3,7),rep(4,7),rep(5,7),rep(6,7),rep(7,7)) 

但对于大数而不是7,说100.

我怎么能这样做而不重复到100？

谢谢

我的例子:

x = c(1:15)
y = c(-4:10)
plot(x, y)

R xlim自动给出:2, 4, 6, ..., 14

我怎样才能改变它: 1, 2, 3,..., 15

我试过了 :

my.limits = as.numeric(seq(1, 15, by = 1))
x = c(1:15)
y = c(-4:10)
plot(x, y, ylim = c(-4,10), xlim = my.limits)

但有错误:

Erreur dans plot.window(...) : valeur 'xlim' incorrecte

我有一个10×10的矩形散点图(说)我想在这个散点图中的每个点周围绘制一个半径= 3的圆.你有什么想法吗？

谢谢你们.

我正在尝试使用 R 进行多项式回归lm，我有

Time=c(1980:2016)
y= rnorm(length(Time))

我用了 ：

reg=lm(y~poly(Time,2))
round(reg$coefficients,3)

这使 ：

 (Intercept) poly(Time, 2)1 poly(Time, 2)2 
        -0.110         -1.298          0.172

和 ：

Time2=Time^2
reg2=lm(y~Time+Time2)
round(reg2$coefficients,3)

它给

(Intercept)        Time       Time2 
   1146.590      -1.128       0.000 

哪里有问题？