search for: ymax

Displaying 20 results from an estimated 84 matches for "ymax".

Did you mean: max
2013 May 18
3
bar plot with non-zero starting level
Hi, I want to plot grouped bars to compare 95% confidence interval estimates from two models. Each bar represents a 95% confidence interval estimate of a coefficient from one of the two models. Each group represents confidence interval estimates of the same coefficient from the two models. I think such a bar plot will nicely present whether 95% confidence interval estimates of the same
2003 Jul 24
5
inverse prediction and Poisson regression
...poisson()) (that's why you see the "dose[-1]" term. The "first" dose in the dose vector is 0. This is really a nice fit. I can obtain a nice slope (B) and intercept (A): log(Y) = B log(x) + A I do have a biological value for dose = 0 from my "control". i/e Ymax = some number with a Poisson error again So, what I want is EC50x : Y/Ymax = 0.5 = exp(B log(EC50x) + A) / Ymax exp((log(0.5) + Log(Ymax)) - A)/B) = EC50x That's all fine, except I don't have a clue on how to calculate the confidence intervals of EC50x or even if I can model...
2001 Oct 23
0
problems with postscript device
...Imedian <- quantile(meds, c(0.025, 0.5, 0.975)) CImean <- quantile(means, c(0.025, 0.5, 0.975)) hmean <- hist(means, freq = FALSE, plot = FALSE) hmeds <- hist(meds, freq = FALSE, plot = FALSE) split.screen(c(1, 2)) xlims <- range(c(hmean$breaks, hmeds$breaks)) ymax <- max(c(hmean$density, hmeds$density)) screen(1) plot(hmean, freq = FALSE, col = "red", xlim = xlims, ylim = c(0, ymax), xlab = paste("Sampling dist. of mean of ", xname), main = "") bringToTop() if (norm) { n <- length(...
2013 Feb 03
1
ggplot2 plotting errorbars.
Hi, i'm using this lines of code: dodge <-position_dodge(width=0.9) ggplot(dfm,aes(x = X,y = value)) + geom_bar(aes(fill = variable), position=dodge, stat="identity") + geom_errorbar(aes(ymin=value-er, ymax=value+er),width=0.25, position=dodge,stat="identity") to plot this data frame X variable value er 1 A X4 58.74 9.44 2 B X4 52.41 10.01 3 C X4 95.52 4.88 4 A X1 75.51 8.54 5 B X1 0.73 23.20 6 C X1 96.66 1.18 7 A X5 76.70...
2002 Feb 11
0
profile
...PopStd <- 90468 beta <- Beta DIs <- Dls MnmN <- Nmin OptN <- Nopt1 # define model function Y.model <- function(gN, MnmN, OptN, DIs, beta, eta1, eta2, Popn, Dmax, AWC, SumEp, PotYield3, Nsupply) { Ymax<- 1-ifelse(Popn<=PopStd, eta1, eta2)*log(Popn/PopStd) Ymax <- Ymax*PotYield3*Popn/1000 Ymax <- Ymax*ifelse(Dmax<=DIs*AWC, 1, 1 - beta*(Dmax -DIs*AWC)/SumEp) Nstar <- (Nsupply- MnmN*Ymax) / (OptN*Ymax - MnmN*Ymax)...
2012 Jun 02
1
mtext bug
...y margin text when i run the loop. when i test the code outside of the loop the margin text appears - however when i run the loop the margin text does not appear ... what gives? here is the code for (col in 1:ncol(xdf)){ par(mar=c(3.5, 4.5, 2, 1), oma=c(2,0,0,0)) ymin = min(xdf[,col]) - 20 ymax = max(xdf[,col]) + 20 png(paste(drive, names(xdf)[col], '.png', sep="")) plot(coredata(yr09_10[,col]), type='l', ylim=c(ymin, ymax), col=1, lwd=2, las=1, xlab="", ylab="", xaxt='n', main=names(xdf)[col]) lines(coredata(yr10_11[,col]...
2006 Jun 25
0
Fighting with GDD & lwd
...ing and installing the GDD package under my FreeBSD 6.1 box after eliminating ifdef DEBUG ? ?PKG_CFLAGS+=-DJGD_DEBUG endif in Makevars.in Now, with the following code .............................................. grafico1<- function() { plot(1:length(tabella[,1]),tabella[,2],ylim=c(Ymin,Ymax),type="h", lwd=6,col="red",main=paste("TEMPERATURE MEDIE GIORNALIERE\nMESE DI", toupper(format(ierid,"%B")),sep=" "), col.main="dark blue", xlab="Giorni del mese",ylab="?C", col.lab="blue",axes=FALS...
2011 Jul 01
1
highlighting clusters in a heatmap
I would like to draw horizontal or vertical lines on a heatmap to highlight the clusters at some specified cut depth of the dendrogram. As a hacked example, the following code would work if I could set the coordinates of the top and bottom of the false color image correctly (ymin and ymax), but the correct values seem to depend on the output device and its size. I realize that heatmaps use a 2x2 layout which makes the coordinate system non-obvious, but the result seems very difficult to customize. I would appreciate any suggestions for manual or pre-made solutions. Example: set....
2010 Feb 10
6
prompts and running means
...he deep end with a something my company want me to write. my code is as follows: kenttemp=read.csv("mnowak.11.1.csv") rows=nrow(kenttemp)-5 kent=kenttemp[1:rows,] #have to remove the last 5 lines of the graph as they interfere with rest of data max(kent[,16],na.rm=TRUE)->ymax ymax=ymax+200 #This is to get vertical scale to fit more accurately Ann=kent$Ann kent=as.matrix(kent) barplot(Ann,main="Annual Monthly Rainfail Data", xlab="Year",ylab="Rainfall (mm)",ylim=c(0,ymax),col="blue",space=0,names.arg=c(kent[,2]), cex.names=0.8)...
2011 Jan 11
0
modified FAST Script from package SensoMineR for the R community - Reg
...ali.sup, function(k) seq(nbre.grpe * (k - 1) + 1, length = nbre.grpe))), 1]) xlim <- c(xmin, xmax) * 1.1 } else { xmin = xlim[1] xmax = xlim[2] } if (is.null(ylim)) { ymin <- ymax <- 0 if (is.na(test.invisible[1])) ymin <- min(ymin, coord.ind[, 2]) if (is.na(test.invisible[1])) ymax <- max(ymax, coord.ind[, 2]) if (!is.null(res.mfa$ind.sup) & is.na(test.invisible[2])) ymin &lt...
2010 May 01
1
ggplot2's geom_errorbar legend
Hello, I create a simple ggplot that only shows a straight line. I then add three datasets of CI using the geom_errorbar function. The problem is that I can't find any way to have the legend showing up ... I need to show what each color of the CIs corresponds to i.e. which method. Can anyone advice please? TIA, Best regards, Giovanni
2011 Oct 20
4
Randomized Points in space/ saving model results
A bit new to R and I'm working on something a bit more challenging than I am used to- so here's whats going on: Data inputs: 9 different shapefiles (.shp) of different point locations (lat, long) shapefile bounding box (lat/long corner points (14, 1) (15,1) (14, 2) (15,2)) 1 .csv of lat/long coordinates for points
2013 May 14
1
Tamaño plots y calidad en grafico ggplot
...quot;a","Nivel 2","Nivel 3","c","d") niveles=data.frame(inicio,fin,Niveles) Jurisdiction=rownames(mtcars) for (i in 1:length(Jurisdiction)) if(Jurisdiction[i]=="Fiat 128"){k=i} h <- ggplot(tabla, aes(x=factor(Jurisdiction), y=Mean, ymin=245, ymax=305)) h + scale_fill_manual(values = c("#A0551E", "#FA964B","#E1EBFF", "#B4CDFF", "#055014","#A0551E"))+ geom_hline(yintercept = seq(245, 305, 10), linetype=2, colour="#E1EBFF")+ geom_vline(xintercept = seq(1, 33, 2.5), linetype...
2001 Feb 28
1
Spider diagrams
Hi, I am trying to create what are called "spider diagrams" in the geochemical literature using R. A spider diagram is basically a plot of the atomic number versus the concentration on a log scale. Lines are drawn from each atomic number for each sample. Right now, my data frame looks like: SAMPLE SITE V3 LA.NASC LAATNUM CE.NASC CEATMNUM 1 1A:001
2020 Oct 26
0
How to shade area between lines in ggplot2
Hi Put fill outside aes p+geom_ribbon(aes(ymin = slope_1*x + intercept_1 - 1/w[2], ymax = slope_1*x + intercept_1 + 1/w[2]), fill = "blue", alpha=0.1) The "hole" is because you have two levels of data (red and blue). To get rid of this you should put new data in ribbon call. Something like newdat <- trainset newdat$z <- factor(0) p+geom_ribbon(data=new...
2020 Oct 23
2
How to shade area between lines in ggplot2
...ot;) ``` Why there is a hole in the middle of the ribbon? and the color is not grey... On Fri, Oct 23, 2020 at 2:35 PM PIKAL Petr <petr.pikal at precheza.cz> wrote: > > Hi > > What about something like > > p+geom_ribbon(aes(ymin = slope_1*x + intercept_1 - 1/w[2], > ymax = slope_1*x + intercept_1 + 1/w[2], fill = "grey70", alpha=0.1)) > > Cheers > Petr > > > -----Original Message----- > > From: Luigi Marongiu <marongiu.luigi at gmail.com> > > Sent: Friday, October 23, 2020 11:11 AM > > To: PIKAL Petr <petr....
2016 Apr 21
2
Calcular Error en modelo lineal
...Confint(RegModel.1, level=0.95) Estimate 2.5 % 97.5 % (Intercept) 0.6100000 -6.8296312 8.049631 x 0.9328571 0.4919135 1.373801 Es decir, el intervalo de confianza de la respuesta y en función de los valores x, sería de la forma (ymin,ymax), siendo: ymin = -6,82 + 0,49*x ymax = 8,05 + 1,37*x Es esto correcto??? Gracias Jesús [[alternative HTML version deleted]]
2009 Feb 01
2
Extracting Coefficients and Such from mle2 Output
The mle2 function (bbmle library) gives an example something like the following in its help page. How do I access the coefficients, standard errors, etc in the summary of "a"? > x <- 0:10 > y <- c(26, 17, 13, 12, 20, 5, 9, 8, 5, 4, 8) > LL <- function(ymax=15, xhalf=6) + -sum(stats::dpois(y, lambda=ymax/(1+x/xhalf), log=TRUE)) > a <- mle2(LL, fixed=list(xhalf=6)) > summary(a) Maximum likelihood estimation Call: mle2(minuslogl = LL, fixed = list(xhalf = 6)) Coefficients: Estimate Std. Error z value Pr(z) ymax 19.2881...
2009 Apr 23
1
ggplot2/aesthetic plotting advice
...rary(ggplot2) nx <- 3 ngrp <- 5 nper <- 4 x <- rep(1:nx,ngrp*nper) y <- runif(nx*ngrp*nper) g <- factor(rep(1:ngrp,each=nx*nper)) dat <- data.frame(x,y,g) se <- with(dat,tapply(y,list(g,x),sd)/table(g,x)) means <- with(dat,tapply(y,list(g,x),mean)) limits <- aes(ymax=means+se,ymin=means-se) gg <- factor(rep(1:ngrp,nx)) xx <- rep(1:nx,each=ngrp) p <- ggplot(dat,aes(y=c(means),x=xx, group=gg,colour=gg)) p + geom_pointrange(aes(ymin=c(means-se),ymax=c(means+se))) xdodge <- xx + rep(seq(-0.1,0.1,length=ngrp),nx) p + geom_line() + geo...
2012 May 19
3
converting csv to image file
Hello everyone, I want to get a 1km by lkm grid raster image using my csv data. If I call latitude=a, longitude=b and preciptation=c. a<-(1,2,3,4,5) b<-(6,7,8,9,10) c<-(10,20, 30,40, 50) Then I found an example in r help which goes like pts = read.table("file.csv",......) library(sp) library(rgdal) proj4string(pts)=CRS("+init=epsg:4326") # set it to lat-long pts =