Re: [R] print multiple plots to jpeg, one lattice and one ggplot2

OK thanks it works christophe 2013/5/3 Felipe Carrillo > Something like this? > library(gridExtra) > grid.arrange(one,two) > > Felipe D. Carrillo > Supervisory Fishery Biologist > Department of the Interior > US Fish & Wildlife Service > California, USA > http://www.fws.gov/redbluff/rbdd_jsmp.

Re: [R] BIOENV

This is a PRIMER question, not an R question. Please take it up with the PRIMER support people. R-helpers cannot be expected to know how other, closed-source, software works. Cheers, Simon. On 06/05/13 12:34, Gilson Carvalho wrote: Dear all, Does anyone knows why the results of a BIOENV (PR

[R] BIOENV

Dear all, Does anyone knows why the results of a BIOENV (PRIMER v. 6.1.15) are diferent of the bioenv() + mantel() in vegan? Not the spearman correlation, indeed the pseudo-p value. I know that the approach bioenv() + mantel() is biased. So, how the BIOENV (PRIMER) ends with larger p values (perm

Re: [R] Open graphics device within sourced script?

Believe it or not, this is FAQ 7.16 (consider also FAQ 7.22). --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#. ##.#. Live Go...

[R] Markov Chain Discrete

Hi, I want to simulate from Markov Chain Discrete Algorithm and my transition  matrix is: The probability transition matrix trans = matrix(c(0.8,0.05,0.05,0,                  0.05,0.80,0.055,                  0.05,0.05,0.81, 0, 0.05, 0.095, 0, 0.81), ncol=4, byrow=TRUE); The Initiate Values are

Re: [R] Open graphics device within sourced script?

Cross-posting to r-help is not appropriate. (This is especially important when the answer was to be found in the R-FAQ.) Please wait for an answer on SO if you choose to post there first. (Or if you have chosen to ignore this request, then at least have the courtesy to post a SOLVED notice here

[R] Open graphics device within sourced script?

Hello all, I am having the weirdest issue that I don't know how to get around. I am using R, and I am sourcing another script within one master script. Within the sourced script, I have code that looks like the following: pdf("Figs/bar_gni.pdf") m1table$Country1 <- reorder(m1table$Country, m1tabl

[R] compare two quadratic regression using bootstrap

Hello, I want to know if two quadratic regressions are significantly different. I was advised to make the test using step 1 bootstrapping both quadratic regressions and get their slope coefficients. (Let's call the slope coefficient *â*^1 and *â*^2) step 2 use the slope difference *â*^1-

Re: [R] Error in plot.window(...) : need finite 'xlim' values

On 05/05/2013 10:20 PM, shaqifahsani wrote: Hi, I am very new to R. I am planning to plot a dendrogram from a binary file but received error after running it. I retrieved the code from the Internet and customized it accordingly... Error: Error in plot.window(...) : need finite 'xlim' values Ca

Re: [R] slope coefficient of a quadratic regression bootstrap

Hello all, I found the following code worked. code require(boot) function.fit2 <- function(data,i){ d <- data[i,] # select obs. in bootstrap sample fit2 <-lm(logDIS~logBMN+I(logBMN^2), data=dataN) coefficients(fit2) # return coefficient vector

[R] slope coefficient of a quadratic regression bootstrap

Hello, I want to know if two quadratic regressions are significantly different. I was advised to make the test using step 1 bootstrapping both quadratic regressions and get their slope coefficients. (Let's call the slope coefficient *â*^1 and *â*^2) step 2 use the slope difference *â*^1-

Re: [R] Package survey: singularities in linear regression models

Ok, that's more helpful. The problem is with replicate-weight designs, and it's because svyglm() uses the fitted coefficients from the point estimate as starting values for fitting the replicates. And even if that is changed, the computation of the replicate variance doesn't like all the replicat

Re: [R] Error in plot.window(...) : need finite 'xlim' values

On 05.05.2013 14:20, shaqifahsani wrote: Hi, I am very new to R. I am planning to plot a dendrogram from a binary file but received error after running it. I retrieved the code from the Internet and customized it accordingly... Error: Error in plot.window(...) : need finite 'xlim' values Can

Re: [R] Lasso Regression error

On May 4, 2013, at 10:26 PM, Preetam Pal wrote: > Thanks David for the paper, I understand the theory. > > But my question is about R only: the vector of coefficients that R outputs in > lars(), does it apply against the original variable y or against (y-y_bar). I > have put in intercept=T as

[R] negative correlation with corSymm in nlme or another way?

Hi, I have a nonlinear model with residuals that are negatively autocorrelated at short distances. I can find no spatial correlation structures in nlme that can easily handle negative autocorrelation as most have bounds on parameter values so that correlation is from 0 to 1. First, is there so

Re: [R] custom function that plots other functions- problem

This request is clearly off-topic according to the Posting Guide. Please use the study assistance provided in the context of your assignment. --- Jeff NewmillerThe . . Go Live... DCN:

Re: [R] custom function that plots other functions- problem

Hello, Try the following. Note that argument 'how_many' is not needed, it can be derived from the length of the functions vector. plotter<- function(f, range, quality){ x <- seq(range[1], range[2], quality) y <- sapply(seq_along(f), function(.i) f[[.i]](x)) plot(x, y[, 1]) fo

[R] custom function that plots other functions- problem

Hi everyone! Could you please help me with the following assignment? My aim is to write a custom function that draws the plots of functions submitted by the user. The catch is, that we do not know the exact number of functions the user will select to draw. The program is relatively easy with one

[R] Error in plot.window(...) : need finite 'xlim' values

Hi, I am very new to R. I am planning to plot a dendrogram from a binary file but received error after running it. I retrieved the code from the Internet and customized it accordingly... Error: Error in plot.window(...) : need finite 'xlim' values Can you explain why and what should I do... Than

[R] ArimaLike NaN return

Hi, > sessionInfo() R version 3.0.0 (2013-04-03) Platform: x86_64-unknown-linux-gnu (64-bit) I am experiencing following problem. I am trying to fit arima(2,0,4) model, and optim complaint that armafn provides it with infinite value. I found that it is because sumlog (res[2L]) (see stats/R/ari

Re: [R] Size of a refClass instance

Yes, I agree. How does one conceptually achieve polymorphic behavior without instantiating 100,000s of instances? Perhaps one way around this is to represent the data in an efficient R way -- i.e. a data.frame -- and create a set of re-usable singleton instances of different node types. To pe

Re: [R] Size of label axis

On 05/05/2013 08:24 AM, Jose Narillos de Santos wrote: Hi all, I want to know if the is a possibility to increase the size of the default numbers that appears on axis on a plot, I don´t mean a labbel I mean the number on axis and yaxis. You see when in excel you select a plot and select all and

Re: [R] R 3 and Debian Testing

This might be a question best posed to the r-sig-Debian list. Cheers, MW On May 5, 2013, at 9:55, Lorenzo Isella wrote: > Dear All, > I am using Debian testing on multiple machine machines at home. > This is my source list > > deb http://ftp.ch.debian.org/debian/ testing main contrib non-free

[R] R 3 and Debian Testing

Dear All, I am using Debian testing on multiple machine machines at home. This is my source list deb http://ftp.ch.debian.org/debian/ testing main contrib non-free deb-src http://ftp.ch.debian.org/debian/ testing main non-free contrib deb http://security.debian.org/ testing/updates main contrib