search for: listhttp

Hi R Users, I have produced a simulated scatter plot of y versus x tightly clustered around the 45 degree line through the origin with the following code: x <- seq(1,100) y <- x+rnorm(100,0,10) plot(x,y,col="blue") abline(0,1) Is there some way to generate a 3-dimensional analogue of this? Can I get a similar simulated scatter plot of points in 3 dimensions where the points

Hi, I have a matrix of results of multiple 2x2 chi^2 tests, non- significant tests are marked as TRUE. Is there a function for grouping those variables in a similar way LSD.test from agricolae library does? I reviewed LSD.test's source but it's not helpful for me. This is my matrix: [,1] [,2] [,3] [,4] [,5] [,6] [,7] 1 TRUE FALSE TRUE TRUE TRUE FALSE TRUE 2 FALSE TRUE

dear useRs, i'm working with a mean reverting model of the following specification: y = mu + beta(x - mu) + errorterm, where mu is a constant currently I estimate just y = x (with lm()) to get beta and then calculate mu = estimated intercept / (1-beta). but I'd like to estimate mu and beta together in one regression-step and also get the test-statistics (including parameter

Dear R-users, I have installed the latest version of Tinn R (Version 2.2.0.2). With the older versions it was possible to run Tinn R next to R and then select code that one could send to R. Now it seems that one has to start the R gui in Tinn R and then it is possible to select en run code of R. But if I start R within Tinn R it does NOT start the R-Gui but only the R console. When I start R

Dear list, I have a question and I'm going to give an example of my problem f<- function(d1, d2, d3) { d<- d1*d2/d3 return(d) } v1<- 1 v2<- 2 If I try f(v1, v2, v3) Error in f(v1, v2, v3) : object "v3" not found I obviously got the above error message. I would like to add something to my function to allow me to get a certain value (say zero) if one of the arguments is

Hi, I am looking for a R package to draw 3d plot. But not in a static way like scatterplot3d or stuff like that. I would like to make rotate the plot, to zoom/unzoom etc..., actually to render the graph 'dynamic'... (for the biologist and bioinformatician, a R package which permitts to reproduce what pdb viewer does). Does a Way exist to do something like that in R ? Thanks,

Hi R helpers, If I have a dataset looks like: ID record 1 20 . 30 . 25 2 26 . 15 3 21 4..................... And I want it becomes ID record 1 20 1 30 1 25 2 26 2 15 3 21 4..................... That is, I have to duplicate IDs for those with multiple records. I am wondering it is possible to be

Hi, dear R users I am a newbie in R and I wantto use the method of meximum likelihood to fit a Weibull distribution to my survival data. I use "optim" as follows: optim(c(1, 0.25),weibull.like,mydata=mydata,method="L-BFGS-B",hessian = TRUE) My question is: how do I setup the constraints so that the two parametrs of Weibull to be pisotive? Or should I use other function

I'm importing a lot of text tables of data (from Latent Gold) that includes hashes in some of the column names ("Cluster#1", "Cluster#2", etc.). Is there an easy way to strip the offending hashes out before pushing the text into a table or data frame? I thought I'd use gsub, e.g., but can't figure out how to read in a text file without reading it into a table or

Good day everyone,   I tried to find a multivariate GARCH package and failed to find one. Although when I searched R I found the following link which describes the package:   http://www.r-project.org/user-2006/Slides/Schmidbauer+Tunalioglu.pdf   can any one help me with this issue.   Thank you in advance [[alternative HTML version deleted]]

Hi I am trying to read genbank est flat files. Is there any function in R to do that. I was looking for something similar to SeqIO.parse(input_handle, "genbank") in bioPython that convert the genebank sequence into fasta sequence. (One option is to call biopython from R but I am not familiar with python functions so I want to avoid that and I am keeping it as last resort) Thanks and

Dear useRs, I'm working on multiple comparison design on two factor (2 3 levels) ANOVA. Each of the tests I have tried (Tukey, multcomp package) seem to do only with one factor at a time. fm1 <- aov(breaks ~ wool * tension, data = warpbreaks) tHSD <- TukeyHSD(fm1, "tension", ordered = FALSE) $tension diff lwr upr p adj M-L -10.000000 -19.35342

Hi list, I want to use the p-value from the survdiff function (package survival) to reuse within a function in a Kaplan-Meier plot. The p-value is somehow not a component of the value list ?! Thanks in advance -- A. Goralczyk G?ttingen, Ger.

I would like to Summarize values that are repeated measures at a certain river mile with box plot i.e. The data matrix looks like this 123 124 125 #fiver mile 0.5 0.6 0.7 0.4 0.5 0.6 ... ... ... #values I would like to make a boxplot with the river mile naming the different box plot. How do you suppress the X123? Stephen -- Let's not

Dear list, Can somebody tell me how to invoke user-defined functions from script files during run-time? Basically I have (almost) one function per script file. I am thinking of something like #include in C++ import in Java/Python use in Perl No, I don't need my functions every time R is started. Neither I am thinking about writing building my functions into packages. I just

Howdy Gurus I 'd like to ask a question about how to build TermDocMatrix in tm text mining package. It is not clear about importing a plain text file, and them converting that text file into TermDocMatrix file, etc to me. How can I build a TermDocMatrix of " a plain text document file for text association? Or are there any good manuals? Thank you in advance, -- Kum-Hoe Hwang, Ph.D.

Dear All: I am trying to read the data from the different sheet but in the same Excel file. Does R provide this kind of command to do it? Anyone has experience on this? Thanks! Xin [[alternative HTML version deleted]]

What do you think of this: http://www.microway.com/whisperstation/whisperstation-r.html I'm considering ditching my Windows Vista 2 GB RAM computer for WhisperStation R using Debian 64-bit Linux with 32 GB RAM and setting the whole thing up for R and WinBUGS. I put in a price request, but I know nothing about Linux, or WhisperStation R for that matter, and am really curious what you think?

Hi, I was trying to copy a directory from one location to another. Is there a command to do this (similar to the file.copy command)? thanks! [[alternative HTML version deleted]]

I need to have the maths symbol for >= in the legend, and to substitute threshold variable with its value. Somehow, various attempts weren't successful. Please help. threshold <- 0.5 plot(NA, xlab="", ylab="", main="", axes=F, xlim=c(0,1), ylim=c(0,1), xaxs="i", yaxs="i") legend(x=0, y=1, fill=c("orange", "white",