我正在从数据库中提取数据,并通过rehsape2运行结果.由于某些原因,这会将POSIXct日期时间戳记变为数字.没问题我想,除了我一小时外,你可以把它们拒之门外.
这是一个最小的例子
foo<-as.POSIXct("2011-04-04 14:18:58")
as.numeric(foo) #gives 130192318
bar<-as.POSIXct(as.numeric(foo),
tz=Sys.timezone(),
origin=as.POSIXct(
strptime("1970-01-01 00:00:00", "%Y-%m-%d %H:%M:%S", tz="UTC")))
as.numeric(bar) #gives 130192318 identical !
foo #Gives "2011-04-04 14:18:58 BST"
bar #Gives "2011-04-04 13:18:58 UTC"
显然foo和bar在数值上是相同的,但是R认为foo需要显示为BST并且显示为UTC.如何将两者显示为BST.这也不起作用;
as.POSIXct(bar, tz="BST") #still gives "2011-04-04 13:18:58 UTC"
作为开发新的分析化学方法的一部分,我需要计算一些数据的运行差异和运行差异.我还需要使用R语言从这些数据中获得置信区间
我想我需要使用anova或其他东西？
我的数据就像
> variance
Run Rep Value
1 1 1 9.85
2 1 2 9.95
3 1 3 10.00
4 2 1 9.90
5 2 2 8.80
6 2 3 9.50
7 3 1 11.20
8 3 2 11.10
9 3 3 9.80
10 4 1 9.70
11 4 2 10.10
12 4 3 10.00
我想要一个包含可变数量表的sweave文档.我认为下面的例子可行,但事实并非如此.我想循环遍历列表foo并打印每个元素作为它自己的表.
%
\documentclass[a4paper]{article}
\usepackage[OT1]{fontenc}
\usepackage{longtable}
\usepackage{geometry}
\usepackage{Sweave}
\geometry{left=1.25in, right=1.25in, top=1in, bottom=1in}
\listfiles
\begin{document}
<<label=start, echo=FALSE, include=FALSE>>=
startt<-proc.time()[3]
library(RODBC)
library(psych)
library(xtable)
library(plyr)
library(ggplot2)
options(width=80)
#Produce some example data, here I'm creating some dummy dataframes and putting them in a list
foo<-list()
foo[[1]]<-data.frame(GRP=c(rep("AA",10), rep("Aa",10), rep("aa",10)), X1=rnorm(30), X2=rnorm(30,5,2))
foo[[2]]<-data.frame(GRP=c(rep("BB",10), rep("bB",10), rep("BB",10)), X1=rnorm(30), X2=rnorm(30,5,2))
foo[[3]]<-data.frame(GRP=c(rep("CC",12), rep("cc",18)), X1=rnorm(30), X2=rnorm(30,5,2))
foo[[4]]<-data.frame(GRP=c(rep("DD",10), rep("Dd",10), rep("dd",10)), X1=rnorm(30), X2=rnorm(30,5,2))
@
\title{Docuemnt to test putting a variable number of tables into a sweave Document}
\author{"Paul Hurley"}
\maketitle
\section{Text}
This document was …我正在尝试使用Sweave生成统计报告,并且我试图在xtable上添加标题,但是如果我包含百分比符号则会破坏它.
一些示例代码
<<label=Analyte2_Results, results=tex, echo=FALSE>>=
print(xtable(result[[2]], caption=paste(levels(vardata$Analyte)[1],
" percent bias and precision estimates with 95 \% confidence intervals",sep="")),
tabular.environment='longtable',
latex.environments=c("center"),
floating=FALSE,
include.rownames=FALSE)
@
如果我使用%或\%那么Sweave过程可以工作,但是Latex没有,如果我尝试\%然后Sweave失败.似乎没有另一种方法可以在Latex中指定百分比符号
有没有办法在标题中实际设置百分比符号？
我有一个数据帧,这是一个运行历史.一些变量包括日期(在POSIXct中)和该运行的值(此处=大小).我想生成各种图表,显示基于特定日期范围的大小列总数的行.理想情况下,我想使用相同的数据集,并从每周,2周,每月季度的总数进行更改.
这是一个示例数据集;
require(ggplot2)
set.seed(666)
seq(Sys.time()-(365*24*60*60), Sys.time(), by="day")
foo<-data.frame(Date=sample(seq(today-(365*24*60*60), today, by="day"),50, replace=FALSE),
value=rnorm(50, mean=100, sd=25),
type=sample(c("Red", "Blue", "Green"), 50, replace=TRUE))
我可以创建这个显示个人价值的情节;
ggplot(data=foo, aes(x=Date, y=value, colour=type))+stat_summary(fun.y=sum, geom="line")
或者我可以这样做以显示每月的金额;
ggplot(data=foo, aes(x=format(Date, "%m %y"), y=value, colour=type))+stat_summary(fun.y=sum, geom="line", aes(group=type))
然而,每季度/ 2周总和会变得更复杂.理想情况下,我喜欢像stat_bin和stat_summary这样的东西,所以我可以指定一个binwidth(或者让ggplot根据范围进行最佳猜测)
我错过了一些明显的东西,或者这是不可能的？
R的帮助系统很棒,我喜欢几乎所有功能的帮助页都包含一个例子,但有时我希望有一种方法可以看到示例代码的结果,而无需复制并粘贴到我当前的会话中.这对图形功能来说是双倍的.
有办法做以下其中一项;
我已经看到了一些关于这个的问题,但是无法解决如何做我想做的事情.
默认情况下,Sweave通过连接Rnw文件的名称和图形对象标签的名称来创建图形.
从这个问题(Make Sweave + RweaveHTML将所有图形放在指定的文件夹中)如果我想将所有图形放在文件夹foo中并命名为bar-graphic我可以使用
\SweaveOpts{prefix.string=foo/bar}
但是如何在文件夹foo中获取图形但命名为rnwfile-graphic？
我有一个很长的脚本来做一些数据分析,它必须连接到几个数据库.其中一个数据库往往比我更频繁地更新我的密码,所以我想为我添加一个弹出框来输入我当前的密码.我找到了这个例子,但是当对话框被解除时,不明白tcltk看到hwo返回我的值.我曾想过将它定义为<<-OnOK函数中的全局变量(),但这看起来很混乱
require(tcltk)
tt<-tktoplevel()
Name <- tclVar("Password")
entry.Name <-tkentry(tt,width="20",textvariable=Name)
tkgrid(tklabel(tt,text="Enter Password"))
tkgrid(entry.Name)
OnOK <- function()
{
NameVal <- tclvalue(Name)
tkdestroy(tt)
}
OK.but <-tkbutton(tt,text=" OK ",command=OnOK)
tkbind(entry.Name, "<Return>",OnOK)
tkgrid(OK.but)
tkfocus(tt)
我有一个数据框,其中包含(除其他外)具有浓度的数字列和具有状态标志的因子列.此状态标志包含NA.
这是一个例子
df<-structure(list(conc = c(101.769, 1.734, 62.944, 92.697, 25.091, 27.377, 24.343, 55.084, 0.335, 23.280), status = structure(c(NA, NA, NA, NA, NA, NA, 2L, NA, 1L, NA), .Label = c("<LLOQ", "NR"), class = "factor")), .Names = c("conc", "status"), row.names = c(NA, -10L), class = "data.frame")
我想用一个字符串替换浓度列以获取标志列的某些值,或者将浓度值格式化为一定数量的有效数字.
当我试试这个
ifelse(df$status=="NR","NR",df$conc)
状态标志中的NA不会触发真或假条件(并返回NA) - 正如文档所暗示的那样.我可以循环遍历行,然后在每个行上使用IF,但这似乎效率低下.
我错过了什么吗？我试过as.character(df $ status)也不行.我的魔力必须变低......
我使用devtools开发了一个小 R 包并将其放在 github ( https://github.com/paulhurleyuk/PaulHurleyMisc ) 上。我用这个结构创建了它
PaulHurleyMisc\
R\
data\
man\
因为它是一个 git 存储库，所以根目录中还有一个 .git 文件夹，当我执行 R CMD Check 时，R 会抱怨可执行文件。
我应该在包的根目录下创建一个子文件夹吗？我还没有尝试使用 devtools 从 github 安装，那么 devtools 会期待什么？