> seq(1,4,1)
[1] 1 2 3 4
> seq(1,4,-1)
Error in seq.default(1, 4, -1) : wrong sign in 'by' argument
>
On Tue, Feb 21, 2012 at 1:53 PM, bioinformatics <uccstud...@live.ie> wrote:
> Here is my code
>
>
> slidingwindowplotATGC = function(windowsize, inputseq)
> {
>
> starts = seq(1, length(inputseq)-windowsize, by = windowsize)
> n = length(starts)
> chunkGs = numeric(n)
> chunkAs = numeric(n)
> chunkTs = numeric(n)
> chunkCs = numeric(n)
> for (i in 1:n) {
> chunk = windowsize[starts[i]:(starts[i]+9999)]
> chunkG = sum("g" == chunk)/length(chunk)
> chunkA = sum("a" == chunk)/length(chunk)
> chunkT = sum("t" == chunk)/length(chunk)
> chunkC = sum("c" == chunk)/length(chunk)
> chunkGs[i] = chunkG
> chunkAs[i] = chunkA
> chunkTs[i] = chunkT
> chunkCs[i] = chunkC
> }
>
>
> plot(starts,chunkGs,type="b",ylim=c(min(min(chunkAs),min(chunkTs),min(chunkCs),min(chunkGs)),max(max(chunkAs),max(chunkTs),max(chunkCs),max(chunkGs))),col
> = "red")
> points(starts,chunkTs,col = "blue")
> points(starts,chunkAs,col = "green")
> points(starts,chunkCs)
> }
>
>
> Im getting the following error message,
>
> Error in seq.default(1, length(inputseq) - windowsize, by = windowsize)
> :
> wrong sign in 'by' argument
>
>
> which I never got before when running codes of this sort, infact I re ran
> old code that worked perfectly before, except this time Im getting this
> error message which doesn't seem to make any sense at all! I need help with
> this before I go completely insane... Maybe Im just bad at this program,
> but
> it seems to me that it has a mind of its own... I was also getting an error
> message before regarding the ylim function, stating that it needed to be a
> finite value, which is what I was giving it? HELP!!!
>
> --
> View this message in context:
> http://r.789695.n4.nabble.com/In-R-2-14-1-what-does-wrong-sign-in-by-argument-mean-tp4407926p4407926.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html<http://www.r-project.org/posting-guide.html>
> and provide commented, minimal, self-contained, reproducible code.
>
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
