小编lmo*_*lmo的帖子

我正在使用该quantreg包在R中运行以下分位数回归:

bank <-rq(gekX~laggekVIXclose+laggekliquidityspread+lagdiffthreeMTBILL+
lagdiffslopeyieldcurve+lagdiffcreditspread+laggekSPret, tau=0.99)

并通过提取系数和摘要统计

bank$coefficients
summary(bank)

我得到的结果是

Call: rq(formula = gekX ~ laggekVIXclose + laggekliquidityspread + 
lagdiffthreeMTBILL + lagdiffslopeyieldcurve + lagdiffcreditspread + 
laggekSPret, tau = 0.99)

tau: [1] 0.99

Coefficients:
                       Value    Std. Error t value  Pr(>|t|)
(Intercept)            -0.03005  0.01018   -2.95124  0.00319
laggekVIXclose          0.00471  0.00069    6.81515  0.00000
laggekliquidityspread  -0.01295  0.01619   -0.79976  0.42392
lagdiffthreeMTBILL     -0.12273  0.12016   -1.02136  0.30717
lagdiffslopeyieldcurve -0.13100  0.06457   -2.02876  0.04258
lagdiffcreditspread    -0.21198  0.15659   -1.35377  0.17592
laggekSPret            -0.01205  0.46559   -0.02588  0.97936

但是,我想知道R ^ 2 /调整后的R ^ …

我试图在这个盒子图上添加两组男性和女性的平均年龄标签.到目前为止,我只能按组进行,而不是按性别和组进行.

我的数据框:

Age=c(60, 62, 22, 24, 21, 23) 
Sex=c("f", "m", "f","f","f","m")
Group=c("Old", "Old", "Young", "Young", "Young", "Young")

aging<-data.frame(Age, Sex, Group)

和剧情的命令:

ggplot(data=aging, aes(x=Group, y=Age))+geom_boxplot(aes(fill=Sex))
+geom_text(data =aggregate(Age~Group,aging, mean), 
aes(label =round(Age,1), y = Age + 3), size=6)

我想使用该quantreg包从分位数回归中提取系数以及上限和下限。这是帮助文件中的示例。

data(engel)
attach(engel)
taus <- c(.05,.1,.25,.75,.9,.95)
f <- rq((foodexp)~(income),tau=taus)
sf <- summary(f)
sf[1]
#[[1]]

#Call: rq(formula = (foodexp) ~ (income), tau = taus)

#tau: [1] 0.05

#Coefficients:
#            coefficients lower bd  upper bd 
#(Intercept) 124.88004     98.30212 130.51695
#income        0.34336      0.34333   0.38975

我知道我可以用来coefficients()获取系数。

cf <- t(data.frame(coefficients(f)))    # transpose for better arrangement
cf
#              (Intercept)    income
#tau..0.05   124.88004 0.3433611
#tau..0.10   110.14157 0.4017658
#tau..0.25    95.48354 0.4741032
#tau..0.75    62.39659 0.6440141
#tau..0.90    67.35087 0.6862995
#tau..0.95    64.10396 0.7090685

但我不知道如何获得出现在summary(). 我看了 …

当我试图在我的工作笔记本电脑上安装软件包时,我收到了一个错误.我被告知我拥有笔记本电脑的完全权限,所以应该没有块.

我认为它可能与安装包所在的库的文件路径有关？我正在使用最新版本的R和R Studio,64位,并使用加拿大的U of T CRAN镜像站点.

以下是我尝试安装plyr时出现的内容

install.packages("plyr")

还安装了只有源代码形式的依赖'Rcpp'包,可能需要编译C/C++/Fortran:'Rcpp''plyr'

当我尝试安装时ggplot2,这是一个很长的错误消息:

错误:依赖'stringi'不适用于包'stringr'*删除'C:/RPackages/R-3.2.3/library/stringr'在install.packages中警告:运行命令'"C:/ RPACKA~1/R -32~1.3/bin/x64/R"CMD INSTALL -l"C:\ RPackages\R-3.2.3\library"C:\ Users\MICHAE~1.COL\AppData\Local\Temp\RtmpGs24zA/downloaded_pa​​ckages/stringr_1.0.0.tar.gz'在install.packages中有状态1警告:软件包'stringr'的安装具有非零退出状态错误:依赖'colorspace'不适用于软件包'munsell'*删除'C:/ RPackages /R-3.2.3/library/munsell'在install.packages中警告:运行命令'"C:/RPACKA〜1/R-32〜1.3/bin/x64/R"CMD INSTALL -l"C:\ RPackages\R-3.2.3\library"C:\ Users\MICHAE~1.COL\AppData\Local\Temp\RtmpGs24zA/downloaded_pa​​ckages/munsell_0.4.3.tar.gz"在install.packages中有状态1警告:安装包' munsell'具有非零退出状态错误:依赖'摘要','plyr','reshape2','scale'不适用于包'ggplot2'*删除'C:/ RPack age/R-3.2.3/library/ggplot2'inight.packages中的警告:运行命令'"C:/RPACKA〜1/R-32〜1.3/bin/x64/R"CMD INSTALL -l"C:\ RPackages\R-3.2.3\library"C:\ Users\MICHAE~1.COL\AppData\Local\Temp\RtmpGs24zA/downloaded_pa​​ckages/ggplot2_2.0.0.tar.gz"在install.packages中有状态1警告:安装包'ggplot2'具有非零退出状态下载的源包位于'C:\ Users\michael.colangelo\AppData\Local\Temp\RtmpGs24zA\downloaded_pa​​ckages'中

"非零退出状态"是什么意思？是否正在尝试为安装找到一个库？>代码中提到的appData文件似乎不存在或者我无法访问它.如何将其更改为另一个本地文件夹？

突然，R Studio停止自动进入列表并输入数据框内容，mylist$然后按Tab。我得到的结果消息是

无匹配

我在非常简单的列表或数据框上尝试过：

 simpledataset <- data.frame(a=rep(0,100),b=rep(1,100))   

但是，当我加载数据集Mtcars时，它会正确显示列名。

我%localappdata%\RStudio-Desktop通过重命名来重设，但未进行任何更改。

有什么可能的原因的想法吗？

我想ggplot2只将科学记数法的第一部分写到轴上，然后在轴的顶部添加一个 $x 10^n$ 数量级。有没有一个功能可以做到这一点？

这是一个带有 hack 的 MWE 来说明我的意思：

ggplot(data = data.frame(x = 1:10, y = seq(1, 2, l = 10)*1000), aes(x,y)) + geom_line()

而我想要的是：

ggplot(data = data.frame(x = 1:10, y = seq(1, 2, l = 10)*1000), aes(x,y)) + geom_line() +
  scale_y_continuous(breaks = c(1, 1.25, 1.5, 1.75, 2, 2.05)*1000, label = c(1, 1.25, 1.5, 1.75, 2, "x 10^3"))

作为一个附带问题，我注意到当刻度标签很大时，轴标签会很快接近刻度标签。有没有办法在它们之间设置自动间距？

我正在使用该earlywarnings包，并且想编辑该qda_ews函数中编写的函数之一。我可以fix(...)，但我想编辑的功能由于某种原因在我使用时没有列出fix。

该函数称为generic_RShiny. 这是 github 的链接（ https://github.com/earlywarningtoolbox/earlywarnings-R/blob/master/earlywarnings/R/qda_ews.R）。

如何访问整个qda_ews.R代码以进行我需要的更改？

我正在处理一个常规需要计算t分布密度的问题,而R中的尾部相当远.

例如,使用R的t分布函数,dt(1.424781, 1486, -5)返回[1] 2.75818e-10.我的一些最终输出(使用此密度作为输入)与我的同事在MATLAB中执行的类似计算的参考值不匹配,我认为这可能是由于R中t分布的尾部不精确.

如果我与MATLAB的t分布函数进行比较,则nctpdf(1.424781, 1486, -5)返回ans = 4.3932e-10,这与R的输出有很大不同.

编辑:R打印两条相同的警告消息

In dt(1.424781, 1486, -5) : full precision may not have been achieved
in 'pnt{final}'

这是在Mac,R版本3.3.1上

我想从日期中仅提取季度,例如,从"2003-02-08"日期获得整数1.我一直在尝试这条线

library(mondate)
as.yearqtr(dat$DATE)

"2003 Q1"

as.character(as.yearqtr(dat$DATE))[1]

"2003 Q1"

没有给出我想要的结果.当然我可以写如下条件

library(data.table)
data$DATE = as.Date(data$DATE, format='%d%b%Y')
data$month=month(data$DATE)
setDT(data)[month==1,  quarter:=1]  
    ...

这会起作用,但根本不优雅.这样做有更美妙的方式吗？

谢谢lmo和user2100721!我真的希望我能接受所有的答案!

我试图隐藏我使用openxlsxR 中的包创建的 Excel 工作簿中的列，但我没有任何运气。我可以制作一个预先格式化的 Excel 工作簿并向其中写入数据，但需要隐藏的列数是可变的。我试过使用：

 setColWidths(wb, sheet = "Sheet 1", cols = col_list, widths = 0)

其中“col_list”是需要隐藏的列的列表。当我打开保存的工作簿时，列宽实际上是 0.38。如果有人对我如何处理此事有任何想法，将不胜感激。