similar to: Help on function please

Dear R users: I encountered difficulties in michaelis-menten equation. I found that when I use right model definiens, I got wrong Km vlaue, and I got right Km value when i use wrong model definiens. The value of Vd and Vmax are correct in these two models. #-----right model definiens-------- PKindex<-data.frame(time=c(0,1,2,4,6,8,10,12,16,20,24),

Dear R users: I have a problem about catch the value from function. I have following two functions (part): sbolus1 <- function() { ....... for( i in 1:Subject) { kel<-par1 Vd<-par2 PKindex<-sbolus1.out(PKtime,kel,Vd,defun,par1,par2,Dose,i) } savefile(PKindex) } sbolus1.out<-function(PKtime,kel,Vd,defun,par1,par2,Dose,i) { time<-PKtime$time

Dear All,   #I have the following code   Dose<-1000 Tinf <-0.5 INTERVAL <-8 TIME8 <-matrix(c((0*INTERVAL):(1*INTERVAL))) TIME7 <-matrix(c((0*INTERVAL):(2*INTERVAL))) TIME6 <-matrix(c((0*INTERVAL):(3*INTERVAL))) TIME5 <-matrix(c((0*INTERVAL):(4*INTERVAL))) TIME4 <-matrix(c((0*INTERVAL):(5*INTERVAL))) TIME3 <-matrix(c((0*INTERVAL):(6*INTERVAL))) TIME2

Dear All   I have the following data (somewhat simplyfied):   TINF <-1 a <-c(500,750,1000,1250,1500,1750,2000) b <-c(8,12,18,24,36,48,60,72,96)   following function:   infcprodessa <-function (D, tin, tau, ts)   (D * (1 - exp(-0.048 * tin))/(tin * (0.048*79) * (1 - exp(-0.048 * tau)))) * exp(-0.048 * (ts - tin)) z <-sapply(1:1, function(n) infcprodessa(1000,TINF,12,12-TINF))   is

Dear All,   Currently I am running the following code:   library(stats4) library(odesolve) library(rgenoud) Input<-data.frame(SUB=c(1),time=c(0.5,3,10,15),lev=c(2.05,12.08,9.02,8)) XD<-500 IT<-3 diffeqfun<-function(time, y, parms) {   if(time<=IT)      dCpdt <- (XD/IT)/parms["Vol"] -

I have a simple graph (oode below) which looks fine on the screen but when I save it in png format the title (actually the last mtext line) is cut off. I am pretty sure that I am doing something very stupid but other than playing around with the png height and width commands which don't seem to help I have no idea of what to look for. Any suggestions or cures would be appreciated. Thanks

Dear list, I'm migrating a project from Matlab to R, and I'm facing a relatively complicated problem for nls. My objective function is below: >> objFun <- function(yEx,xEx,tEx,gamma,theta,kappa){ yTh <- pdfDY(xEx,tEx,gamma,theta,kappa) sum(log(yEx/yTh)^2) } The equation is yTh=P(xEx,tEx) + noise. I collect my data in: >> data <-

https://bugs.freedesktop.org/show_bug.cgi?id=56692 Priority: medium Bug ID: 56692 Assignee: nouveau at lists.freedesktop.org Summary: Xorg freezes with vblank enabled on nouveau nvc8 QA Contact: xorg-team at lists.x.org Severity: normal Classification: Unclassified OS: Linux (All) Reporter: kel.p.doran at

Hi group, I've been thinking of calculating euclidean distance between each column of a data frames that each consists of standardized numerical columns. However, I don't know if there's a way of summarizing the overall distance by some kind of metrics. If anyone know a proper way of doing so and/or a package I would greatly appreciate your suggestions. Thanks very much! Kel --

Hi there, I am trying to use linear regression to solve the following equation - y <- c(0.2525, 0.3448, 0.2358, 0.3696, 0.2708, 0.1667, 0.2941, 0.2333, 0.1500, 0.3077, 0.3462, 0.1667, 0.2500, 0.3214, 0.1364) x2 <- c(0.368, 0.537, 0.379, 0.472, 0.401, 0.361, 0.644, 0.444, 0.440, 0.676, 0.679, 0.622, 0.450, 0.379, 0.620) x1 <- 1-x2 # equation lmFit <- lm(y ~ x1 + x2) lmFit Call:

Hi group, I did a quick search in the archive but still couldn't find a function that performs nested logistic regression. Please point me in the right direction. Thanks a million! Regards, Kel

CPU 99.3%: 2170 root 25 0 11656 3228 10m R 99.3 0.5 0:07.42 smbd We're having the same problems as many others with print queues, client lock-ups, etc. I think maybe this CPU problem is related as we have two servers in two locations with the same versions, and both had the problem immediately after upgrading to 3.0.10. We use the YUM facility to upgrade Samba. To react to these

Hi Group, I wonder if R has the function CensorReg which takes care of left truncation. I know it's in S. Thank you so much! Kel __________________________________________________ Do You Yahoo!? Get personalized email addresses from Yahoo! Mail - only $35 a year! http://personal.mail.yahoo.com/ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help

Hi group, I did some searches about this very simple question. Hope someone can help me out. If I have the following: a <- "x <- 2^2" a [1] "x <- 2^2" How do I evaluate the expression that gets me an answer of 4? I tried the following: > eval(a) [1] "x <- 2^2" > get(a) Error in get(a) : object 'x <- 2^2' not found (apparently

Dear R users: I have a problem about catch the value from function. I have following two functions (part): fbolus1 <- function() {......... par<-data.frame(Parameter=c("kel","Vd"),Initial=c(0)) check(par) .....} check<-function(par) { if (par[ ,2] <= 0){ cat("\nEnter again (y/n) ?\n\n") ans<-readline() if (ans == "n" ){

I am working on minimization of sum of squared errors for a problem that has 2 box-constrained parameters. I got the solution for this problem using "L-BFGS-B" method in optim function using an R code as res<-optim(par=c(parInit), fn=myFunction, method = c("L-BFGS-B"), lower = parMin, upper = parMax,

My institute has been heavily dependent on SAS for the past while, and SAS is starting to charge us a very deep amount for license renewal. Since we are a non-profit organization that is definitely not sustainable. The team is brainstorming possibility of switching to R, at least gradually. I am talking about the entire institute with considerable number of analysts using SAS their entire

https://bugs.freedesktop.org/show_bug.cgi?id=54437 Bug #: 54437 Summary: linux-nouveau2.6 (3.6.0-rc4) : GTX580 : Xorg freezes when using accel Classification: Unclassified Product: xorg Version: git Platform: x86-64 (AMD64) OS/Version: Linux (All) Status: NEW Severity: critical

Hello, I have a problem with loader. I recently upgraded from 6_rel to 7_rel. Now when I install world there is a problem booting. Here is what I do: cd /usr/src make buildworld make buildkernel KERNCONF=BLACK make installkernel KERNCONF=BLACK At this point I can reboot and all is good. After boot I install the new world: cd /usr/src mergemaster -p reboot into single user mode cd /usr/src make

I'm running ubuntu. and I'm trying to get qtracker running in wine. www.qtracker.com is the website. I'm a total noob when it comes to wine on linux. I'm really a noob to linux. But I can listen to someone talk fast so dont be afraid to help. It should be an easy install for someone who's already used wine. I'm just wondering I've gotten as far as $wine qtracker.exe