search for: cmax

Hello, I have a problem with ggplot and colors I used this function to draw somes lines and I want them to be all black (just to test) however, I dont get any black line in the figure. Do you have any experience with ggplot and manual colors?? ggplot(cmax, aes(cmax[,3], cmax[,6],colour=interaction(cmax[,1],cmax[,2]))) + geom_line() + geom_point() + scale_fill_manual(values=c("black", "black" , "black", "black", "black", "black", "black", "black", "black", &qu...

...# compute weight-based vd (L) cl.pt <- wtkg * cl # compute CL (L/hr) k <- cl.pt/vd.pt # compute k (hr^-1) cave <- round((twhi - twlo)/log(twhi/twlo),2) doserate <- (cl.pt * cave)/F # maintenance dose (Dm) total_dose_tau <- round(doserate * tau,0) Cmax <- total_dose_tau/vd.pt Cmin <- Cmax * exp(-k * tau) AR <- 1/(1 - exp(-k * tau)) Cmaxss <- Cmax * AR Cminss <- Cmin * AR Cfluctss <- (Cmaxss - Cminss) Ct <- (Cmaxss * exp(-k * t)) dose_loading_vd <- round(Cmaxss * vd.pt,0) dose_loading_dm <- round...

...,2,2,1,1,2,2,1,1,2,2,1,1,2,2,1,1,2,2,1,1,2,2,1,1) e<-c(1739,1633,1481,1837,1780,2073,1374,1629,1555,1385,1756,1522,1566,1643, ???? 1939,1615,1475,1759,1388,1483,1127,1682,1542,1247,1235,1605,1598,1718 ) KK<-data.frame(subj=as.factor(a), drug=as.factor(b), per=as.factor(c),? seqe=as.factor(d), Cmax=e) M<- lm(Cmax ~ seqe+ subj:seqe + per + drug , data=KK) anova(M) drop1(M, test="F") The result of Type I SS: Analysis of Variance Table Response: Cmax ????????? Df Sum Sq Mean Sq F value Pr(>F) seqe?????? 1??? 585???? 585? 0.0160 0.9014 per??????? 1? 63175?? 63175? 1.7293 0.2131...

...nction in the evd package appears to contain an error. I am using evd version 2.2-4 in R 2.14.0 with Windows XP. > data(portpirie) > mrlplot(portpirie) ## No Error > tlim <- c(3.6, 4.2) > tcplot(portpirie, tlim) ## Line from documentation Error in fpot(data, u[1], model = model, cmax = cmax, r = r, ulow = ulow, : `x' must be a non-empty numeric vector > tcplot(portpirie, tlim=c(3.6, 4.2)) ## explicitly specifying the threshold limits Error in fpot(data, u[1], model = model, cmax = cmax, r = r, ulow = ulow, : `x' must be a non-empty numeric vector tcplot(port...

...the + #    the value of cc that has expected power + # tol1   the allowance between the expceted power and target power + #--------------------------- + pc<-(0:((1-d)*nc))/nc + ncl<-length(pc) + ne.old<-ne + ne.old1<-ne.old + ### sweeping forward + for(i in 1:ncl){ + cmin<-cc[1] + cmax<-cc[2] + ### fixed cci<-cmax bug + cci <-1 + lhood<-dbinom((i:ncl)-1,nc,pc[i]) + ne[i:ncl]<-(1+(cci-1)*(lhood/lhood[1])) * ne.old1[i:ncl] + Ep0 <-Epower(nc, d, ne, pc, alpha) + while(abs(Ep0[i]-power)>tol1){ + if(Ep0[i]<power) cmin<-cci + else cmax<-cci + cci<-(cmax...

...5.80, 5.80, 5.77, 5.72, 5.70, 5.66, 5.62, 5.62, 5.61, 5.59, 5.57, 5.51, 5.47, 5.41, 5.39, 5.40, 5.40, 5.42, 5.45, 5.47, 5.50, 5.50, 5.49, 5.43, 5.39, 5.33, 5.26) day <- c(1:100) df <- data.frame(day,walevel) library(evd) clusters(df$walevel, u = 5.80, r = 1, ulow = 5.35, cmax = T, plot = T) clusters(df$walevel, u = 5.80, r = 50, ulow = 5.35, cmax = T, plot = T) What have I done wrong? Tonja

...model in R, following the reccomendations of Soetart & Herman (A practical guide to ecological modeling). However, I've obtained a number of problems, that I'm unable to solve (even thoughI've tried during at least 2 hours). Can anyone help me? Many thanks > ccrit <- 0 > cmax <- 5 > cclass <- ccrit:cmax > nc <- length(cclass) > tmax <- 10 > times <- 1:(tmax-1) > nprey <- 3 > dive <- c(0.02,0.1,0.14) > pfood <- c(200,60,40) > f <- matrix(nrow=cmax,ncol=nc ,0) > bestpatch <- matrix(nrow=tmax-1,ncol=nc-1,0) >...

...on -8, which I'll be looking into soon. After some intense debugging, I found the problem. One block in the file triggered a very rare bug in the LPC coefficient quantizer caused by insufficient floating point precision. There is a snippet to compute the log(base 2) of a number: floor(log(cmax) / M_LN2) It turns out on at least some x86's compiling with gcc-2.91 (at least) the result of floor(log(8.0) / M_LN2) is 2.0 not 3.0. Anyway, I made a fix and also added code to detect and handle similar problems (related to FP accuracy). Josh _____________________________________________...

> From: Jean-Louis Abitbol > > Dear All, > > I would like to merge a data frame such as: > > > basetab > subject dose cmax > 1 1031 50 21.8 > .... > > with the result of a by processing such as: > > > tmax<-by(pkga,subject,f.tmax) > > tmax > subject: 1031 > [1] 6 > -------------------------------------------------- > > f.tmax being a function written by H. Nguyen (...

...et = 1 giving row names to the H data. 2 - Import H from the data file and set the row names using row.names(H) <- x where x is the data/time data. The problem is when I try to use the cluster function on either I get the following error: > Clusters <- clusters(H, 3.0, r=144, plot=TRUE, cmax=TRUE, row.names=TRUE) Error in xy.coords(x, y, xlabel, ylabel, log) : (list) object cannot be coerced to type 'double' In addition: Warning message: In min(diff(xdata)) : no non-missing arguments to min; returning Inf Does anyone know how to set date/time row names in the cluster func...

Hi R users, I have a question about adding uncertainty bars to stacked bar plots. DF: year A B C Amin Amax Bmin Bmax Cmin Cmax 2009 40 45 15 30 61 23 56 14 17 2010 36 41 23 26 54 22 51 22 24 I use the code below: DF.refm = melt(subset(DF[,c(1:4)]),id.vars='year',variable_name='Legend') fig1 = ggplot(data=DF.refm, aes(x=year,y=value,fill=...

...gt; #include <linux/pm_runtime.h> #include "pci.h" @@ -2518,14 +2519,29 @@ static unsigned int pci_scan_child_bus_extend(struct pci_bus *bus, { unsigned int used_buses, normal_bridges = 0, hotplug_bridges = 0; unsigned int start = bus->busn_res.start; - unsigned int devfn, cmax, max = start; + unsigned int devfn, fn, cmax, max = start; struct pci_dev *dev; + int nr_devs; dev_dbg(&bus->dev, "scanning bus\n"); /* Go find them, Rover! */ - for (devfn = 0; devfn < 0x100; devfn += 8) - pci_scan_slot(bus, devfn); + for (devfn = 0; devfn < 0x100...

...gt; #include <linux/pm_runtime.h> #include "pci.h" @@ -2518,14 +2519,29 @@ static unsigned int pci_scan_child_bus_extend(struct pci_bus *bus, { unsigned int used_buses, normal_bridges = 0, hotplug_bridges = 0; unsigned int start = bus->busn_res.start; - unsigned int devfn, cmax, max = start; + unsigned int devfn, fn, cmax, max = start; struct pci_dev *dev; + int nr_devs; dev_dbg(&bus->dev, "scanning bus\n"); /* Go find them, Rover! */ - for (devfn = 0; devfn < 0x100; devfn += 8) - pci_scan_slot(bus, devfn); + for (devfn = 0; devfn < 256;...

...gt; #include <linux/pm_runtime.h> #include "pci.h" @@ -2518,14 +2519,29 @@ static unsigned int pci_scan_child_bus_extend(struct pci_bus *bus, { unsigned int used_buses, normal_bridges = 0, hotplug_bridges = 0; unsigned int start = bus->busn_res.start; - unsigned int devfn, cmax, max = start; + unsigned int devfn, fn, cmax, max = start; struct pci_dev *dev; + int nr_devs; dev_dbg(&bus->dev, "scanning bus\n"); /* Go find them, Rover! */ - for (devfn = 0; devfn < 0x100; devfn += 8) - pci_scan_slot(bus, devfn); + for (devfn = 0; devfn < 256;...

On 18/06/17 12:10, lily li wrote: > Hi R users, > > I have a question about adding uncertainty bars to stacked bar plots. > > DF: > year A B C Amin Amax Bmin Bmax Cmin Cmax > 2009 40 45 15 30 61 23 56 14 17 > 2010 36 41 23 26 54 22 51 22 24 > > I use the code below: > > DF.refm = melt(subset(DF[,c(1:4)]),id.vars='year',variable_name='Legend') > > fig1 = ggp...

...=3,byrow=TRUE) 1 30.0 41.1 1 45.0 51.5 1 60.0 54.0 1 75.0 56.2 1 90.0 59.1 1 120.0 59.3 time_vib <- vib[,2] diss_vib <- vib[,3] outfile <- "work" cat("WORK","\n",file=outfile,"\n",append=TRUE) n_vib <- length(time_vib) p_vib <- 3 cmax <- 30 WLSiter <- 500 beta0_vib <- list(Wmax_star= 30, tau_star=15, gamma_star=.8) t0 <- 20 tol <- 10^-8 meanfunc_vib <- function(time_vib,Wmax_star,tau_star,gamma_star){ diff<-time_vib-t0 f_vib<-Wmax_star*(1-exp(-log(2)*((diff/tau_star)^gamma_star))) f_vib # meanf...

...developer this softwere is very much useful thanks for the everyones availability first. I need the programm to run an IVIVC for an IR BCS class 2 product. I have dissolution upto 45 mins and plasma profile upto 72 hours. For test and reference products of pilot study which is actually failed for cmax parameter. To predict the upcoming formulations of the same product with varying release rates I need atlease rank order IVIVC program please help me. If possible please share example how to compile the data in .csv file. Thanks in advance. Regards, Manimaran [[alternative HTML version deleted...

...- else - q = (FLAC__int32)(error - 0.5); -#else q = lround(error); -#endif + #ifdef FLAC__OVERFLOW_DETECT if(q > qmax+1) /* we expect q==qmax+1 occasionally due to rounding */ fprintf(stderr,"FLAC__lpc_quantize_coefficients: quantizer overflow: q>qmax %d>%d shift=%d cmax=%f precision=%u lpc[%u]=%f\n",q,qmax,*shift,cmax,precision+1,i,lp_coeff[i]); @@ -234,14 +241,7 @@ int FLAC__lpc_quantize_coefficients(const FLAC__real lp_coeff[], unsigned order, #endif for(i = 0; i < order; i++) { error += lp_coeff[i] / (1 << nshift); -#if 1 /* unfortunately...

Basic x86 support [1] for running Linux as secondary Jailhouse [2] guest is currently pending in the tip tree. This builds on top and enhances the PCI support for x86 and also ARM guests (ARM[64] does not require platform patches and works already). Key elements of this series are: - detection of Jailhouse via device tree hypervisor node - function-level PCI scan if Jailhouse is detected -

Basic x86 support [1] for running Linux as secondary Jailhouse [2] guest is currently pending in the tip tree. This builds on top and enhances the PCI support for x86 and also ARM guests (ARM[64] does not require platform patches and works already). Key elements of this series are: - detection of Jailhouse via device tree hypervisor node - function-level PCI scan if Jailhouse is detected -