小编ale*_*ulo的帖子

我将使用两个变量计算时间序列的一阶导数(dpH/dtime),time和pH.

在R中是否有任何类型的函数可以执行此操作,还是应该计算额外的函数来执行此操作？

我有这个数据集:

x <- c(0, 40, 80, 120, 160, 200)
y <- c(6.52, 5.10, 4.43, 3.99, 3.75, 3.60)

我使用lm()以下方法计算了线性模型

model <- lm(y ~ x)

我想知道x如果我有新y值的预测值,例如ynew <- c(5.5, 4.5, 3.5),但如果我使用该predict()函数,它只计算新y值.

x如果我有新y值,如何预测新值？

我用这段代码制作了下面的叠加条形图:

library(ggplot2)
library(dplyr)

Year      <- c(rep(c("2006", "2007", "2008", "2009"), each = 4))
Category  <- c(rep(c("A", "B", "C", "D"), times = 4))
Frequency <- c(168, 259, 226, 340, 216, 431, 319, 368, 423, 645, 234, 685, 166, 467, 274, 251)
Data      <- data.frame(Year, Category, Frequency)

Data2 <-Data%>%
  group_by(Year, Category)%>%
    summarise(Sum_grp = sum(Frequency))

Data3 <-transform(Data2, Pos = ave(Frequency, Year, FUN = cumsum) - Frequency / 2)

ggplot(Data3, aes(Year, Frequency, group=Category,fill = Category))+
  geom_bar(stat="identity")+
  geom_text(aes(label = Frequency,y=Pos), size = 3) 

现在,我想在每个栏的顶部添加每个组的总和,但我不知道如何.

有人可以帮帮我吗？

非常感谢!!!!!

我使用lda()MASS包中的函数进行了线性判别分析.现在我将尝试绘制一个像ade4包(forLDA)中的双标图.你知道我怎么能这样做吗？

如果我尝试使用该biplot()功能它不起作用.例如,如果我使用Iris数据并生成LDA:

dis2 <- lda(as.matrix(iris[, 1:4]), iris$Species)

然后我可以使用该函数绘制它plot(),但如果我使用该函数biplot()它不起作用:

biplot(dis2)
Error in nrow(y) : argument "y" is missing, with no default

如何绘制变量的箭头？

我正在尝试将sd函数应用于我的数据框但它不起作用:

sdsd <- by(nowna[, 1:16], nowna$stamm, sd)
Error in FUN(X[[1L]], ...) : could not find function "FUN"

你知道为什么吗？

非常感谢.

如果我在 Ipython notebook 中输入此代码，它会起作用：

plt.scatter(x,y)

但是如果在 python IDLE 中输入它，那么我会收到这个错误：

<matplotlib.collections.PathCollection object at 0x0000000007957908>

你知道为什么吗？

非常感谢。