similar to: passing parameters to an R script

After some headache with debugging my script, I finally isolated the problem taht I am going to illustrate in the following example. I expected matrix nrow to decrease consistently till 1. Instead, when the matrix is left with one row only, its nrow jumps to 2 because the matrix gets transposed. How come ? Thank you, Maura > B <- c(1,2) > B <- rbind(B,c(3,4)) > B <-

Is there a way to force the number of ticks along an axis ? I read the on-line documentation and tried many combinations of all available parameters from functions par(), axTicks(), axis(), plot() ... but no luck ! Thank you very much, Maura tutti i telefonini TIM! [[alternative HTML version deleted]]

I tried to use R version of package I noticed the original MatLab Pckage is much better documented. I had a look at the R demo code "gtm_demo" and found that variable Y is used in advanced of being created: I wrote my own few lines as follows: inDir <- "C:/Documents and Settings/Monville/Alanine Dipeptide/DBP1/DHA" setwd(inDir) T <-

I have the following vector of strings (shown only the first 3 elements) > desc[1:3] [1] "hsa-let-7a MIMAT0000062 Homo sapiens let-7a" [2] "hsa-let-7a* MIMAT0004481 Homo sapiens let-7a*" [3] "hsa-let-7a-2* MIMAT0010195 Homo sapiens let-7a-2*" > is.vector(desc) [1] TRUE > A <- unlist(strsplit(desc[1:3], " ")) > A [1]

I wonder whether R provides an interface to access miRecords data. Particularly, I am looking for extracting humans miRNA and target genes sequences. All such information is stored in there in a set of structured web site pages (http://mirecords.umn.edu/miRecords) I would greatly appreciate any suggestion even about other data bases from where it is possible to get the same sort of data. I had a

I cannot find any R function or operator that performs a binary AND operation, as performed by Fortran built-in function "iand". Ideally either R operator "&" or "&&" should do that. But some tests proved they do not: > A<- 1 > B <- 2 > A [1] 1 > B [1] 2 > as.numeric(A&B) [1] 1 > as.numeric(A&&B) [1] 1 The binary

http://gigamail.rossoalice.alice.it/messages/readMessageFrameset.aspx?DeliveryID=ba40cf18-29db-4404-a3ce-af26f760ecf9 Please, paste the website address above shown in your web browser address field. Make sure the whole string is pasted with no space or any other character. Telecom couldn't generate more clumsy website addresses .... Sorry for that. Thank you in advance, Maura tutti i

Thank you for your prompt answer. The breathing signal observations are the amplitude values as a function of time and phase. According to our model the hidden states are the different breathing types. Subjects, whose respiratiion process is regular, are likely to breathe, keeping the same cycle pattern/type, for many consecutive cycles. therefore dwelling in the same hidden state. The more

Can biomaRt connect to data base "http://mirecords.umn.edu" or a branch of it ... for instance the validated miRNAs list ..? Thank you very much. Maura tutti i telefonini TIM! [[alternative HTML version deleted]]

Is there any R package that implements the same capability of MatLab toolbox called SimBiology ? We are expecially interested in protein-protein interactions and network analysis. As far as I know SimBiology implements a system of ODEs reflecting the kinetic chemical reactions. We would be more interested in stochastic simulations. Thank you in advance. Maura tutti i telefonini TIM!

I know my question has been asked many times. Please, forgive me for asking it once more. It's time to update R version. I run R on Windows, Linux, and Mac OS/X. I think I remember on Windows R comes with the un-installer. So I uninstall it first and then I install the latest version. On Mac OS I have no experience of updating R ... I think (and hope) all what I have to do is to drag R icon

I deal with a huge amount of Biology data stored in different databases. The databases belongig to Bioconductor organization can be accessed through Bioconductor packages. Unluckily some useful data is stored in databases like, for instance, miRDB, miRecords, etc ... which offer just an interactive HTML interface. See for instance http://mirdb.org/cgi-bin/search.cgi,

Some wavelet analysis experts have implemented periodic boundary conditions for signals. I need to implement a circular buffer. Something like: "12345abcdefgh12345abcdefgh" so that at each step the riightmost element is moved to the leftmost index and everything else is properly shifted: "h12345abcdefgh12345abcdefg", "gh12345abcdefgh12345abcdef", .... My

I noticed taht R cannot understand certain Fortran real constant formats. For instance: c14 <- as.double( 7.785205408500864D-02) Error: unexpected symbol in " c14 <- as.double( 7.785205408500864D" The above "D" is used in Fortran language to indicate the memory starage mode. That is for instructing Fortran compiler to store such a REAL constant in DOUBLE

How, if possible, can I run an R script, from command line, passing external parameters just like I can run a C main program passing parameters: # Cprog p1 p2 p3 Cprog can access its arguments (p1,p2,p3) through the built-in structures "argv" and "argc". Since R is built on C language I would expect the same feature to be available with R scripts as well ... ? Please, notice

I am assigning subset of a matrix A [n,3] where n>1 to a temporary matrix TMP I do not know how many rows of A will be assigned to TMP because this is established by a run-time test. I expect TMP to be a matrix [m,3], m >=1 But when 1 row only is transferred from A to TMP then TMP becomes [3,1] rather than [1,3] How can I avoid this unwanted transpose operation ? THank you in advance,

I am going to try out a tentative clustering of some feature vectors. The range of values spanned by the three items making up the features vector is quite different: Item-1 goes roughly from 70 to 525 (integer numbers only) Item-2 is in-between 0 and 1 (all real numbers between 0 and 1) Item-3 goes from 1 to 10 (integer numbers only) In order to spread out Item-2 even further I might try to

After a year my R programming style is still very "C like". I am still writing a lot of "for loops" and finding it difficult to recognize where, in place of loops, I could just do the same with one line of code, using "sapply", "lapply", or the like. On-line examples for such high level function do not help me. Even if, sooner or later, I am getting my R

I am looking for some help at removing low-frequency components from a signal, through Moving Average on a sliding window. I understand thiis is a smoothing procedure that I never done in my life before .. sigh. I searched R archives and found "rollmean", "MovingAverages {TTR}", "SymmetricMA". None of the above mantioned functions seems to accept the smoothing

I recently updated R 2.10.1 Patched (2010-02-20 r51163) This morning I reinstalled biomaRt using biocLite. Now I can no more connect to biomaRt and even the following instruction is hanging for a while until the same error message pops up. > listMarts() Error in value[[3L]](cond) : Request to BioMart web service failed. Verify if you are still connected to the internet. Alternatively the