小编mne*_*nel的帖子

我用cast()从reshape包装相当频繁.几乎每次都会弹出这个警告:

Aggregation requires fun.aggregate: length used as default

我试图设定options( warn =-1),但无济于事.如何抑制这些警告？

我想制作条形图,其中可以指定最小条形(很像盒子中的方框和胡须图).barplot可以做到吗？我怀疑ggplot中的答案,但我找不到一个例子.这是一些数据:

                    X  Jan  Feb  Mar  Apr  May  Jun  Jul  Aug  Sep  Oct  Nov  Dec
1      Highest recorded 31.5 31.8 30.3 28.0 24.9 24.4 21.7 20.9 24.5 25.4 26.0 28.7
2  Mean monthly maximum 27.8 28.6 27.0 24.8 22.0 20.0 18.9 18.8 20.4 22.4 23.9 26.8
3    Mean daily maximum 24.2 24.8 23.1 20.9 18.4 16.3 15.5 15.7 16.9 18.3 20.0 22.4
4                  Mean 19.1 19.8 18.1 16.2 13.8 11.9 11.2 11.6 12.7 14.1 15.7 17.7
5    Mean daily minimum 14.0 …

Run Code Online (Sandbox Code Playgroud)

我正在使用如下所示的数据框:

date<-c("2012-02-01", "2012-02-01", "2012-02-03", "2012-02-04", "2012-02-04", "2012-02-05", "2012-02-09", "2012-02-12", "2012-02-12")
var<-c("a","b","c","d","e","f","g","h","i")
df1<-data.frame(date,var)

我想创建第二个数据框,将每天观察的数量制成表格.在该数据帧中,未提及的日期将为零...导致如下所示:

date<-c("2012-02-01","2012-02-02","2012-02-03","2012-02-04","2012-02-05","2012-02-06","2012-02-07","2012-02-08","2012-02-09","2012-02-10","2012-02-11","2012-02-12")
num<-c(2,0,1,2,1,0,0,0,1,0,0,2)
df2<-data.frame(date,num)

我已经尝试了一些聚合函数的东西,但无法弄清楚如何包含没有观察的日期(零).

这个想法是无论是下雨还是下雨,潮湿或精细,用户都必须得到他付出的全部费用.

来自Apple:

Store Kit provides built-in functionality to restore transactions for non-consumable products, auto-renewable subscriptions and free subscriptions

对于这些交易,Apple Store Kit具有良好的内置工具.我想专注于其他类型(特别是消耗品).唯一的交易信息是我们在成功购买后通过Store Kit接收的标识符和收据数据.我们的应用程序使用服务器端模型向其提供产品.但是仍然存在许多丢失购买数据的情况,例如当用户通过App Store进行购买时服务器放下,因此无法向服务器发送收据以完成验证过程.

目前的解决方法是:

1. 服务器返回产品标识符列表
2. 用户选择一个; app将其标识符保存在设备上(通过SQLite或Core Data).Standart Apple Store交易流程紧随其后.
3. 如果成功,应用程序将收据数据与其设备上的标识符一起保存并发送给服务器.如果出现故障或取消,则立即从设备中删除标识符.
4. 如果服务器的响应正常,则应用程序将从设备中删除带有收据数据的标识符.否则,它将定期向服务器发送请求,直到成功响应行为为止.

但这种方法仍然存在泄漏.例如,用户可以从设备中删除不等待交易到服务器的应用程序,因此根本不会有任何关于他的购买的证据.

你的建议？

在R中,我们可以运行时间序列的ACF相关图,并且置信区间带将以浅蓝色绘制.但是当我拉动ACF对象的结构时,我找不到这些值.有谁知道如何提取置信区间带的值？

例如

List of 6
 $ acf   : num [1:27, 1, 1] 1 0.06453 -0.06354 0.00213 -0.01324 ...
 $ type  : chr "correlation"
 $ n.used: int 501
 $ lag   : num [1:27, 1, 1] 0 1 2 3 4 5 6 7 8 9 ...
 $ series: chr "tser[i:(i + 500)]"
 $ snames: NULL
 - attr(*, "class")= chr "acf"

使用 R 中的默认对比度（contr.treatment），线性模型对象的摘要根据级别名称给出行名称。当我将对比度更改为 contr.sum 时，线性模型对象的摘要根据组合数字给出行名称。

对于下面的示例代码，处理对比的行名称是 xa xb xc xd xe，对于总对比，它们是 x1 x2 x3 x4 x5。

除了手动重命名行之外，有没有办法使这些行为相同？

例子：

y <- rnorm(10, 0, 1)
x <- factor(rep(letters[1:5], each = 2))

options(contrasts = c("contr.treatment", "contr.poly"))
summary(lm(y ~ x))

options(contrasts = c("contr.sum", "contr.poly"))
summary(lm(y ~ x))

这是一个简单的问题,但它开始让我烦恼,我找不到解决方案....

我希望能够在使用paste或paste0时将它作为输出保留0.0,所以如果我有以下内容:

y <- c(-1.5,-1.0,-0.5,0.0,0.5,1.0,1.5)
> y
[1] -1.5 -1.0 -0.5  0.0  0.5  1.0  1.5
paste0("x",y,"x")

我明白了:

[1] "x-1.5x" "x-1x"   "x-0.5x" "x0x"    "x0.5x"  "x1x"    "x1.5x" 

但想要:

[1] "x-1.5x" "x-1.0x"   "x-0.5x" "x0.0x"    "x0.5x"  "x1.0x"    "x1.5x" 

这似乎是一项简单的任务,但我在执行它时遇到了困难,并且找到了如何在R中完成它的示例.

我有一个大约75K记录的数据框.一个领域有日期,可以追溯到大约11年.我已经使用以下代码去除小时分钟秒,所以我只有一个月 - 一天

 dat$date=round.POSIXt(dat$date,units="day")

我现在正在尝试根据日期创建一个新的字段"FiscalYear",如果日期是在2008-07-01和2009-06-30之间,那么它将在2009-07-01和2010之间分配给... 06-30被分配到2010财年等......

我一直在尝试的方法是使用ifelse> = && <=来设置日期范围,但这不起作用.有什么想法吗？

我是R的新手,目前正在读一本书"广义的添加模型",一本R by Wood(2006)的介绍,并经历了一些练习,特别是关于空气污染和死亡的部分,这是我感兴趣的领域.使用mgcv包我运行以下模型.

library(gamair) 
library(mgcv) 
data(chicago) 

ap1<-gam(death ~ pm10median + so2median + o3median +s(time,bs="cr",k=200)+ s(tmpd,bs="cr"), data=chicago,family=poisson)

如何提取pm10median的效果估计值和x的95%CI并将输出导出为CSV或任何其他选项？

我不确定是否有可能做我想做的事.我想用两个比例绘制一个x轴,如下图所示.非常感谢,