search for: cjlu

Dear R-Friends, How can I calculate a double integral like \int_a^b \int_c^y g(x, y) dx dy where a, b, c are constants, g(x, y), e.g., g(x, y) = tan(x + y). I tried to nested integrate() and adapt(), but none of them working, seemingly due to the limits can not be specified constants. Best regards, C. Joseph Lu Department of Statistics National Cheng-Kung University Voice:

In rw1050, I found that > pweibull(3:10, 2) [1] 0.9998766 0.9999999 1.0000000 1.0000000 NaN NaN [7] NaN NaN Warning message: NaNs produced in: pweibull(q, shape, scale, lower.tail, log.p) more surprisingly, > pweibull(3:10, 2.1) [1] 0.9999566 1.0000000 1.0000000 -Inf NaN NaN [7] NaN NaN Warning message: NaNs produced in: pweibull(q,

...June day 15 language R However, xlab=expression(alpha) correctly shows up on x-label by using previous version: R.1.0.1. Hope this is not a bug for R.1.1.0! Thanks for your help! C. Joseph Lu Department of Statistics National Cheng-Kung University Tainan, Taiwan E-mail: cjlu at email.stat.ncku.edu.tw __________________________________________________________ PC home ?K?O?q?l?H?c?A????????: http://www.pchome.com.tw ?| ?? ?? ?@ ?A ?x ?W ?? ?j ?? ?J ?f ?? ?? -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing...

Dear R-friends, For the following data: > xy x y i 1 731 0.313 2 2 739 0.340 2 3 790 0.373 2 4 855 0.451 2 5 980 0.608 2 6 575 0.156 1 7 608 0.207 1 8 630 0.249 1 9 670 0.332 1 10 838 0.377 1 11 964 0.466 1 > coplot(y ~ x|i, data=xy) coplot gives 3 panels, rather than 2, namely one for i=1 and one for i=2. Futhermore, when I extand data fram xy to have i=3 as follows:

Dear R-Friends, Sorry to bother everyone with my earlier question: "Do we have similar R function to work like Minitab's INDICATOR?" I found the way to make things work: e.g., > x <- c(2,2,5,3,6,5,NA) > model.matrix(~ factor(x) - 1) factor(x)2 factor(x)3 factor(x)5 factor(x)6 1 1 0 0 0 2 1 0 0 0

Dear R and ESS users, I got the following error when trying to load a source file (smart.R in directory E:\R\) into XEmacs through the menu iESS and "Load source file": Error in file(file, "r") : unable to open connection In addition: Warning message: cannot open file `E:Rsmart.R' It seems the slash "/" or "\" being missing, but I don't know

Dear R-Friends, Do we have a similar R function to work like Minitab's INDICATOR? INDICATOR C2, C11-C15 C2 C11 C12 C13 C14 C15 2 1 0 0 0 0 2 1 0 0 0 0 5 ------> 0 0 0 1 0 3 0 1 0 0 0 6 0 0 0 0 1 5 0 0 0 1 0 * * * * *

Dear R-friends, After change (and saved) font style to "bold" in GUI preferences from menu-bar Edit, R 1.2.0, Win98, the "buflines" in the file Rconsole changed to "bufiles" and that causes error in invoking R next time. Is this a bug? Merry Christmas & Best regards, C. Joseph Lu Department of Statistics National Cheng-Kung University Tainan, Taiwan

Dear R-friends, When I selected "Save as postscript" via mouse left-button on the active "R Graphics" window, those colored lines I made in the plot changed in the saved postscript file, e.g., color blue becomes green, red becomes blue, ... I have not changed any setting/options. What's going on? My version of R is > version _ platform i386-pc-mingw32 arch

Dear R-friends, I have whole bunch of objects, say, A, B, ..., Z, all of them are "list" and haveing "Title" (a character string) as their first element. E.g., > A$Title [1] "Airplane Component" What I like to do is to find a way to printout all the Titles at one time. I tried the following: > paste(LETTERS, "Title", sep="$") ->

Dear R-friends, As the binary windows version of R being distributed, what is the correct way to install contributed packages without using winst.exe as before? The answer in "R for Windows FAQ" says: "You can install pre-compiled packages either from CRAN or from a local .zip file by using install.packages: see its help page." But the information available in

Please allow me to add just a little more about this: nothing wrong with pweibull(), namely, the two cases I reported: pweibull(3:10, 2) and pweibull(3:10, 2.1), in rw1041 and earlier version. I wonder this might just due to the change from rw1041 to rw1050, however, I can't find anything relevant (seems to me) in the News or Readme. Thanks Sundar for the suggestion of using 1 -

Dear R-users, The numerical integration example given in Gentleman and Ihaka (2000), "Lexical Scope and Statistical Computing," JCGS, 9, 491-508, is very interesting and helpful in understanding how lexical scope is about. However, I got some questions that I just can't figure out. First all, allow me to copy the two functions given by the authors: midpoint <- function(f, a,