search for: ncku

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,

...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 list -- Read http:/...

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, 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, 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-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

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

I thought I was imaging it but... I have a friend (honestly) that likes to clear his inbox every month by copying older into a folder (e.g. June-2009). Every month he complains that all the mail in the new folder has disappeared. This has happened three times now, and he's sort-of-right. I now believe it's a bug rather than him. The setup is Dovecot 1.1.3 on FreeBSD 7.0 and the client is

Hi all, I have a Samba box (Mandrake 8.1, Samba 2.2.2). I have another hard disc that has Mandrake 9.0 + Samba 2.2.7 on it. I am looking to install this other hard disc into the Mandrake 8.1 box and migrate everything over. Do I need to do anything other than copy the contents of /etc/samba (smb.conf, secrets.tdb etc) over to the new hard disc to ensure that our Windows 2000 Pro workstations

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

Hello!! Does there any packages can simulation data from Markov Random Field? Thanks much Ming Hung Chen [[alternative HTML version deleted]]

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="$") ->

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,