similar to: From table to data.frame

Greetings all! I'm wanting to plot contours of a function, and I don't want to evaluate it at a dense grid of points (i.e. I don't want a huge array of values). Say I have a vector of x-values such as x <- 0.1*c(1:10), and the same for y <- 0.1*(0:10). I then evaluate a 10*10 matrix z of values of z = f(x,y). so I can then do CL <- contourLines(x,y,z) and get a list of

Hello R experts, I am trying to plot ordiellipse for my data but the col according to factors. Metabolites_raw= read.csv(file.choose(), head = TRUE) #file 21Metabolites.csv Metabolites_t=t(Metabolites_raw[,2:82]) ord <- metaMDS(Metabolites_t, distance="bray") symbol=as.numeric(Metab_metadata$LandType) col.list <-

I do have some trouble with matrices. I want to build up a covariance matrix with a hierarchical structure). For instance, in dimension n=10, I have two subgroups (called REGION). NR=2; n=10 CORRELATION=matrix(c(0.4,-0.25, -0.25,0.3),NR,NR) REGION=sample(1:NR,size=n,replace=TRUE) R1=REGION%*%t(rep(1,n)) R2=rep(1,n)%*%t(REGION) SIGMA=matrix(NA,n,n) for(i in 1:NR){ for(j in

When I try to add the following annotation to a plot the entries are plotted one on top of the other. I'm trying to get something that looks like " eta = 0.2 " where the Greek letter is used on the plot. I realize that expression( eta == 0.2) is one solution, but ultimately I'd like to use this in a legend that uses a loop to fill the entries so I don't want

Dear friends of lme, After so many year with lme, I feel ashamed that I cannot get this to work. Maybe it's a syntax problem, but possibly a lack of understanding. We have growth curves of new dental bone that can well be modeled by a linear growth curve, for two different treatments and several subjects as random parameter. By definition, newbone is zero at t=0, so I tried to force the

glm(S ~ -1 + Mdif, family=quasipoisson(link=identity), start=strt, sdat) gives error: > Error in glm.fit(x = X, y = Y, weights = weights, start = start, etastart > = > etastart, : > cannot find valid starting values: please specify some strt is set to be the coefficient for a similar fit glm(S ~ -1 + I(Mdif + 1),... i.e. (Mdif + 1) is a vector similar to Mdif. The error

Dear all, I have a multiline plot with each line labeled with a different letter. But I'm not able to make the legend display the same kind of pattern '-a-', instead the letter is overwritten by the line. A simpler legend with only the letter is not very visible and the pt.bg does nothing with letters. Any idea? plot(1:10,10:1,lty=1,type='b', lwd=2,pch='a')

Cordial saludo. Estoy tratando de hacer un gráfico con qplot de la libreria ggplot2, en ella tengo dos vectores un "x" y un "y" de tal forma que el gráfico seria qlot(x,y...), ahora quiero añadir un nuevo gráfico con un vector "y1" pero en el mismo cuadro algo parecido a lo que hace la funcion "lines", es decir, que en en mismo cuadro me queden las dos

Hi R People: If I have a fitted values from a model, how do I plot the (1-alpha)100% confidence intervals along with the fitted values, please? Also, if the intervals are "shaded" gray, that would be nice too, please? I check confint, but that doesn't seem to do what I want. Thanks in advance, Sincerely, Erin -- Erin Hodgess Associate Professor Department of Computer and

I have two problems with svyhist. I hope this is easy, and it may be that "simple R syntax" will do what I need. 1) I'm able to get a single plot out, but I need to put two distributions in the same graphic to visually compare them. 2) I get uniform breaks at intervals of 10. I can plot the unweighted data using MyBreaks with hist as shown in the second code snippet