search for: inrae

Hi ! I need a plot for data extending over several orders of magnitude on the y axis. The following command generates a nice looking semi-log plot for my data: plot(x,y,log="y",type="l",lty=3, ylim=c(0.01,2),yaxp=c(0.01,1,1),las=1) I would appreciate having also minor tick marks in-between the 3 major ticks obtained with the above command. The "minor.tick" function

Dear all, I am a SAS user, who's trying R. I am a little bit lost for graphics. What is the simplest way for plotting y as a function of x with one symbol (line or dot) for each level of a class variable z (for SAS langage= plot x*y=z). Can I add vertical bars for standard deviation stored in a separate variable (say ystd). Finally, is it possible to do all this with xyplot function

Dear Colleagues, I have x, y data (pollen and seed dispersal from oaks !) that I would like to fit with a function which look like this: p(a,b,x,y)=b/(2*pi*a?gamma(2/b))*exp(-(square_root(x?+y?)/a)power(b)) I am looking for a and b values that fit my data at best. Can someone give me hints to perform such an analysis with R ? Thanks a lot Sophie Sophie Gerber

Dear all, Does any know how to solve PDE with R? The archive list refers to the use of ODE if PDE are parabolic. I am not a mathematician and this does not mean anything for me! help would be very appreciated. Many thanks -- ___________________________________________________ Christophe NGUYEN UMR 1220 INRA-ENITAB Transfert sol-plante et cycle des ?l?ments min?raux dans les ?cosyst?mes

I would like to import "tif" images in R and I do not find any function that can do that. In Matlab there exists the function "imread" that can read the most known images format. Does a similar function exist for R ? Thanks in advance -- Herve CARDOT ____________________________________________________________ Unite Biometrie et Intelligence Artificielle, INRA Toulouse BP

I am not yet a R user but I will be soon. I am looking for the R command and syntax to compute the roots of Bessel function i.e. computing the z values that lead to Jnu(z)=0 where J is a Bessel function or order nu. May You help me ? thanks in advance. Dr Catherine COUTAND Institut National de la Recherche Agronomique (INRA) umr Physiologie Int?grative de l'Arbre Fruitier et Forestier

Hello My question is the following. I have tried to find a similar subject in archives but not found (perhaps bad search!) and I tried R-help and I was advised that R-devel would perhaps better for my question So I wonder if it is possible to create a dll from a R code to use it in another language/program. Indeed, one of my colleague wrote a R code combining some R functions to create a

Dear developer, in my version of R (2.4.0) as weel as in a more recent version (2.6.0) on different computers, we found this problem : > a<-(58/40-1) > a [1] 0.45 > b<-(18/40) > b [1] 0.45 > a<b [1] TRUE > a==b [1] FALSE > Something seems wrong here. but if we do > c<-0.45 > d<-0.45 > c<d [1] FALSE then everything is ok. If we use 59

Hi all, By default, the heatmap function gives an image with a dendrogram added to the left side and to the top. Is it possible to only add the dendrogram to the left side and let the order of the columns unchanged ? I tried heatmap(mat, col=rbg,Rowv=res.hclust$order,Colv=1:dim(mat)[[2]]). In this case, the order of the columns are unchanged but a dendrogram is added to the top. How can I

I'm looking for an R function that simply recodes a quantitative variable into a number of classes according to specified break-points. Obviously I can do this using nested ifelse() commands, but I want to write it into a function where I can't pre-specify the number of classes. Is there an obvious way to do this? An example to clarify: how to convert c(0,10,5,1,9,6) to

Hi all, I'm working out some Fortran code for which I want to compute the Choleski decomposition of a covariance matrix in Fortran. I tried to do it by two methods : 1) Calling the lapack function DPOTRF. I can see the source code and check that my call is correct, but it does not compile with: system("R CMD SHLIB ~/main.f") dyn.load("~/main.so") I get: Error in

Hello R users, I found the odfWeave package to create an odf document. It seems to be a very nice tool. So i tried to used it to create a report with multiple plot: I create an odt file with some code inside: I connect to a mysql database I get a list of projects foreach project I would like to make a plot (a map exactly) then in a R console I use the odfweave (inFile, outFile) function.

I found a problem with scale function, while using center=FALSE and scale=TRUE: all column are not divided by standard error, but divided by sqrt (1/(n-1) sum Xi^2 ) Example: > l<- c(1,2,3) > scale(l,F,T) [,1] [1,] 0.3779645 [2,] 0.7559289 [3,] 1.1338934 attr(,"scaled:scale") [1] 2.645751 2.645751 = sqrt( 1/2 * (1+4+9) ) Antoine Lucas & Aymeric Labourdette

Full_Name: Antoine Lucas Version: 1.7.0 (2003-04-16) OS: Linux Submission from: (NULL) (193.51.197.253) I have a problem using read.table: If in a dataframe, we have a string containing a dot, write.table will not write any file while using option "dec=','". Example > m <- "1.5" > write.table(m,dec=',') Error in if (n%%nrowv == 0) value <-

Hello, I am plotting a 3D function using persp and I would like to use greek symbols in the axes labels. I have found examples like this one on the web: plot(0,0,xlab=expression(kappa[lambda]),ylab=substitute(paste(phi,"=",true,sigma),list(true=5))) this works well with plot but not with persp: with the command persp(M,theta = -20,phi =

Hi all, Let Names a vector of chatacters. For example, > Names [1] "g 604 be-0 -p1 (602 matches)" "g 606 Phli-0 -p2 (517 matches)" [3] "g 608 alu-0 (659 matches)" I try to use gsub or grep functions for two problems : 1. First, I would like to delete all the characters between parentheses. [1] "g 604 be-0 -p1" "g 606 be-0 -p2" [3] "g

Hi, I'm plotting a bar chart like this: ggplot() + geom_bar(data=res,aes(fill=f1,x=f2,y=y),stat="identity",position="dodge") f1 contains quite a few levels and the plot is really quite difficult to read when the order of bars on the graph and on the legend does not match. This problem has been discussed recently here:

Dear R helpers, I have the following model model_1<-glm(y~A+B+C+E+A:D,contrasts=list(D=contrasts_D),data=mydata,na.action=na.omit) with: options(contrasts=c("contr.sum", "contr.poly")) A,B and E are 2-levels factor, C is covariate, D is 20 levels factor with 10 in relation with the first levels of factor A, the other in relation with the second levels of factor A

How can I insert mathematical expressions for variable names in a lattice parallel plot? I tried to implement mathematical expressions in varnames, however, without success. For example, neither parallel(~iris[1:4] | Species, iris, varnames=c("P[Width]", "Petal[length]", "alpha[Width]", "Sepal[Length]")) nor parallel(~iris[1:4] | Species,

Hi, I have some trouble with the result of a fortran function calling the R binomial RNG: The C function rbinom is wrapped as in the file attached. My main fortran program starts as call rndstart() and ends as call rndend() I happen to call the binomial RNG within a loop as b = ggrbinom(1.d0,0.5d0) write(*,*) 'b=',b In certain cases, after a few iterations in the loop, b get