search for: plotbs

Hi Folks, I'm curious for an explanation of the following -- it's a matter of trying to understand how R parses it. I've written sundry little "helper" variants of functions, in particular plot(), to save repetitively typing the same options over and over again. For example: plotb <- function(x,...){plot(x,pch="+",col="blue",...)} This does exactly

Hi, I've read up on readBin() and chapter 6 in the R Data Import/Export manual, but I still can't read a binary file. Here is how the creator of the file described the code that would be needed in Fortran: "Every record has a return in fortran. The length of each record is nx*ny*4. To read you would use the following: nlayx = nx*ny*4 do iz=1,nz,4 read(binary file) var(1:nlayx)

Hi, I need a help. I am new in R and I need to run a nested anova with fixed and random factors (Mixed Model). I have a design with three factors: Day, Area and Plot and the dependent variable is density. The factors Day and Area are fixed while Plot is random, factor Area is nested in factor Day, and factor Plot is nested in Area. I can do it using aov by: mod1<-aov(density~ day +

Dear R-helpers and Klaus, I would like to know how to write out a tree file with bootstrap from phangorn package. That tree file could be in newick format or others. I am new for phylogenetic operation in R. Could you please give me any directions on that? Thanks in advance. Best wishes, Jian-Feng, ######## # as a example # I accomplished 1000 bootstrap simulation on a fit object (a maximum

Please forgive me if this has been asked before but I could not readily find an answer. First, from the example I was able to determine that graphics device commands such as plot can be redirected using 'png' and 'bmp' type commands. I tried this and it works as I understand it. The question is how this graphics device works. If I issue another 'plot' command will it also

Dear R-users, Does anybody know a good way to create and use several graphical devices at the same time in interactive mode? Ideally, I'd like to open 2 to 3 devices and assign names to them. This way, I could make any addition/modification/update to a particular device using its name. I did not see anything like a name argument in ?X11. Is there an alternative? Thanks in advance for your

...lackscholes <- function(s, k, r=.1, t=5, sigma=.9,call=TRUE) { #calculate call/put option d1 <- (log(s/k)+(r+sigma^2/2)*t)/(sigma*sqrt(t)) d2 <- d1 - sigma * sqrt(t) ifelse(call==TRUE,s*pnorm(d1) - k*exp(-r*t)*pnorm(d2),k*exp(-r*t) * pnorm(-d2) - s*pnorm(-d1)) } plotbs <- function(price){ #create s<-seq(0,price,len=price/2) k<-s t<-0:5 sigma<-seq(0,0.9,by=.1) #expand information OptionPrice<-matrix(nrow=length(s),ncol=length(k)) for(i in 1:length(s)) OptionPrice[i,]<-mapply(blackscholes,s[i],k) g...

I am setting up base plots -- one in viewport A and and one in B. This part works fine. But if I go back to A after having done B and add horizontal lines it seems to not use the correct coordinates. How do I tell it to resume using A's coordinates? I am already using par(fig = gridFIG()) but it seems that that's not enough to reestablish them. What happens is that when I go back to

I am setting up base plots -- one in viewport A and and one in B. This part works fine. But if I go back to A after having done B and add horizontal lines it seems to not use the correct coordinates. How do I tell it to resume using A's coordinates? I am already using par(fig = gridFIG()) but it seems that that's not enough to reestablish them. What happens is that when I go back to

Esteemed R Users, Would some one be patient enough to explain the variation i see when creating .pdf file output for my plots? FYI- my goal is produce the highest quality .pdf output from the R 'command line' as opposed to using the menu of the acitve graphics window. Im using 32bit WinXP. Session info at the bottom. I have Ghostscript v8.71 installed and in the Path such that

Version of R: Windows Version 2.0.0 The experimental design contains two plant lines - a control (C) and a mutant (M) - grown out three separate times in plots A, B, C. The design is unbalanced: In plot A, 9 control plants were grown with 29 mutant plants. In plot B, 8 control plants were grown with 20 mutant plants. In plot C, 8 control plants were grown with 22 mutant plants. The