我想传递一个带矩阵的R txt文件,其中从行中省略拖尾零(除了第一个拖尾零,如果有的话).那些缺失的值被视为零.
例如:
8 7 0
5 4 3 2 1
4 8 9
应该理解为:
8 7 0 0 0
5 4 3 2 1
4 8 9 0 0
在读取矩阵之前,最大行大小(即矩阵列的数量)是未知的.
我的重新编码尝试
df$test[(df$1st==(1:3) & df$2nd <= 4)] <- 1
df$test[(df$1st==(1:3) & df$2nd <= 5)] <- 2
df$test[(df$1st==(1:3) & df$2nd <= 6)] <- 3
导致"longer object length is not a multiple of shorter object length"警告和大量的NA以s df$test,即使一些重新编码正常工作.
我错过了什么？任何帮助赞赏.
DW
我有一个每日动物园(xts),有几十年的数据,格式如下:
head(almorol)
1973-10-02 1973-10-03 1973-10-04 1973-10-05 1973-10-06 1973-10-07
183.9 208.2 153.7 84.8 52.5 35.5
我想绘制冬季数据(12月,1月和2月的整月).我找到了xts的子集,所以我想我可以使用以下方法提取所有的Decembers:
x<-apply.yearly(almorol, FUN=last(almorol, "1 month"))
然后在Jan和Feb做类似的事情,但是我收到以下错误:
Error in get(as.character(FUN), mode = "function", envir = envir) :
object 'FUN' of mode 'function' was not found
我可以单独使用apply.yearly和,last(almorol, "1 month")但是当我将它们组合起来时它不起作用.有没有人知道如何对这3个月的时间序列进行子集化？谢谢你的帮助!
欧元符号€显示为...PDF输出
ggplot() + theme_bw() + geom_line() + scale_y_continuous(formatter = "euro")
我在R中使用了样本函数来生成100个公平硬币翻转的样本,这是我使用的命令.
fair.coin = c("heads" = 0.5, "tails" = 0.5)
然后,
x <- sample(fair.coin, size = 100, replace = TRUE)
> x
tails tails tails heads tails tails tails heads heads tails heads heads tails tails
0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5
heads tails heads tails tails tails tails heads heads heads heads tails heads tails
0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5
tails tails tails …想象一下我有一个2列的数据框
Id Value
12 13
32 3
6022 11
9142 231
12 23
119 312
...
我想获取每个“ Id”的平均值。您知道执行此操作的任何快速方法吗？
我在 *.Rnw.log 中遇到以下错误：
Error: chunk 3 (label=printgraph)
Error in grid.Call.graphics("L_text", as.graphicsAnnot(x$label), x$x, :
Metric information not available for this family/device
仅在 Lyx 中并取决于声明
+ scale_y_log10(ylim=c(0.025,5))
在以下示例代码中：
<<Code, tidy=T, sanitize=TRUE>>=
require(ggplot2)
df.qdf=data.frame(T = c(0 ,1,2 ,0 ,1.5,2.5)
,q.500 = c(0.025,2,1 ,0.025, 4,2 )
,q.025 = 0.5*c(0.025,2,1 ,0.025, 4,2 )
,q.975 = 1.5*c(0.025,2,1 ,0.025, 4,2 )
,Occasion = c( 1,1,1 ,2 , 2,2 )
)
theme_set(theme_bw())
Graph=(qplot(T,q.500,data=df.qdf)
+geom_smooth(aes(ymin=q.025,ymax=q.975)
,data=df.qdf
,stat='identity'
,fill='blue')
+ylab('Response')
+xlab('Time [h]')
+facet_grid(facet=.~Occasion)
) +scale_y_log10(ylim=c(0.025,5))
@
<<printgraph,echo=FALSE,fig=TRUE,width=10,height=7,sanitize=TRUE>>=
print(Graph)
@
在下面 …
假设我希望我的R程序的用户输入正数.
如果他输入负数或输入任何字母(a,b,c,d等),我会想让他继续尝试,直到他输入所需的正数.
如何检查他是否输入了一个角色(a,b,c等)？
例如,如果我有:(事实上,这个网站上有人帮我正确编写了这段代码,但我不明白某些特别的东西,我"#explanation of this line"在代码中指出的第三行
n <- -1
while(is.na(n) | (n < 1) ){
n <- readline("enter a positive integer for the number of simulations: ")
n <- ifelse(grepl("\\D",n),-1,as.integer(n)) #explanation of this line
}
问题:
我知道这is.na(n)意味着如果n为null(不可用),那么如果用户按下Enter则is.na(n)成为TRUE.那是对的吗？
是什么意思:n <- ifelse(grepl("\\D",n),-1,as.integer(n)？
如何检查用户是否输入(a,b,c,d,...)？
感谢那些可以帮助我的人.
Owusu Isaac
我有一个fasta格式文件,其中我只需要保留那些长度小于100的节点.但是,我目前面临的问题是我能够分离节点但是无法放置每个节点的字符在单独的变量中,其长度可以检查并随后将必需的节点与较长的节点分开.所以我的意思是我能够读取标题和单独的节点,但我如何将每个节点中的字符放在变量中.
这是我的数据样本
>NODE_1
GTTGGCCGAGCCCCAGGACGCGTGGTTGTTGAACCAGATCAGGTCCGGGCTCCACTGCAC
GTAGTCCTCGTTGGACAGCAGCGGGGCGTACGAGGCCAGCTTGACCACGTCGGCGTTGCG
CTCGAGCCGGTCATGAACGCGGCCTCGGCGAGGGCGTTCTTCCAGGCGTTGCCCTGGGAA
>NODE_2
CCTCCGGCGGCACCACGGTCGGCGAGGCCCTCAACATCCTGGAGCGCACCGACCTGTCCA
CCGCGGACAAGGCCGGTTACCTGCACCGCTACATCGAGGCCAGCCGCATCGCGTTCGCGG
ACCGCGGGCGCTGGGTCGGCGACCCCGCCTTCGAGGACGTAC
>NODE_3
CCTCCGGCGGCACCACGGTCGGCGAGGCCCTCAACATCCTGGAGCGCACCGACCTGTCCA
CCGCGGACAAGGCCGGTTACCTGCACCGCTACATCGAGGCCAGCCGCATCGCGTTCGCGG
ACCGCGGGCGCTGGGTCGGCGACCCCGCCTTCGAGGACGTACATCATTCCTTAATCTTCC
我的代码:
x <- readLines("1.fa", n = -1L, ok = TRUE, warn = TRUE)
for (i in 1:length(x)) {
if (substr(x[i],1,1)=='>') {
head <- c(head,x[i])
q <- x[i+1]
if (q=!0) {
contig <- c(contig,q)
print(contig)
contig.length <- c(contig.length, nchar(q))
} else {
break
}
} else {
z <- paste(z,x[i], sep=" ")
}
}