[R] echo output to screen

2006-02-09 Thread Chun-Ying Lee
Hi all, I want to echo output to screen before next step, like the following when i type library(pkname), library(pkname) a few description about the package how should I do to get this? Thanks in advanced!! __

[R] echo output to screen

2006-02-09 Thread Chun-Ying Lee
On Thu, 09 Feb 2006 10:02:49 +0100, Uwe Ligges wrote Chun-Ying Lee wrote: Hi all, I want to echo output to screen before next step, like the following when i type library(pkname), library(pkname) a few description about the package how should

[R] how to save output all together

2005-12-04 Thread Chun-Ying Lee
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) }

[R] How to catch value

2005-12-03 Thread Chun-Ying Lee
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)

[R] confine the value

2005-10-23 Thread Chun-Ying Lee
Dear R users: I wonder if it is possible to confine the value which optim optimized to be all positive. I use optim with method Nelder-Mead. Please give me some comments. Thanks in advance!! __ R-help@stat.math.ethz.ch mailing list

[R] monte carlo simulation

2005-10-11 Thread Chun-Ying Lee
Dear R user: I wonder if it is possible to run monte carlo simulation with dse2 package(MonteCarloSimulations function) using ordinary differential equation. How do I define the model? Or if there are any functions which can run monte carlo simulation using ordinary differential equation.

[R] summary nls output

2005-09-25 Thread Chun-Ying Lee
Dear R user: I bulid a package, and in the package I use the function nls to solve some questions. If I have two sets of data, and I want to summary these two data's nls output, I write the command in the package source code like: { .. summary(fm1) summary(fm2) } then i compiler the

[R] Extract data from edit chart

2005-09-20 Thread Chun-Ying Lee
Dear R users: I wonder if it is possible to use data.frame and edit to show the form like: parameter upper lower A B C when I type the command edit(abc), assume the data.frame named abc. And then I want to extract the number form it, like abc[1,2] which mean the upper level of

[R] how to overlook the zero in the denominator

2005-07-27 Thread Chun-Ying Lee
Dear R users: I have two set of data, as follow: x-c(0,0,0.28,0.55,1.2,2,1.95,1.85, 1.6,0.86,0.78,0.6,0.21,0.18) y-c(0,0,0,0.53,1.34,1.79,2.07,1.88, 1.52,0.92,0.71,0.55,0.32,0.19) i-1:length(x) I want to sum each (x[i]-y[i])^2/x[i] together, like: Sum -sum((x[i]-y[i])^2/x[i]) Sum [1]

[R] Michaelis-menten equation

2005-07-19 Thread Chun-Ying Lee
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

Re: [R] Michaelis-menten equation

2005-07-19 Thread Chun-Ying Lee
. Thanks. regards, ---Chun-ying Lee Hmm, sorry, no. I'm talking through a hole in my head there. Vm*y/(Km+y) makes OK sense. Vm is what you get for large y - passing from 1st order to 0th order kinetics. However, looking at the data plot(PKindex) abline(lm(conc~time,data=PKindex