similar to: Density Plots

Hi, I found that the R package "KernSmooth" can deal with only 1D and 2D data. But now I have a collection of 4-dimensional data (x1,x2,x3,x4) and would like to estimate the "mode" of the underlying density. What can I do about it ? Thanks a lot. -- Ying-Chao Hung Assistant Professor Graduate Institute of Statistics National Central University Chung-Li, Taiwan TEL:

Sorry for the off-topic question, but I know there are some talented LaTeX users out there. Which bibliography style gives only the year in text citations (e.g "for further details, see Anderson (1992)" )? Thanks Jason -- Indigo Industrial Controls Ltd. http://www.indigoindustrial.co.nz 64-21-343-545 jasont at indigoindustrial.co.nz

Dear all, Our department uses the linux system and we are not allowed to submit job directly. We must make a batch to submit through "qmon". so, I make a foo.sh file, which only contains one line: nohup R --vanilla < foo.txt > foo.results foo is all my codes. It is a simulation of 200 times. I set the seed at the beginning. It is to estimate the success probability, which is

Dear List: I am unsure if my problem is with Sweave or LaTeX. Anyhow, I am using the MikTeX distribution and TexnicCenter. I can easily create Sweave files and all goes well until I try to incorporate graphics. I use the same code as found in the examples found in the users manual. In R, the graphics I want are created as Sweave is creating the .tex file. When I examine the .tex file

Paul E. Johnson wrote: > I handed out some results from glm() and the students ask "how many > observations were dropped due to missing values"? > > How would I know? > > In other stat programs, the results will typically include N and > the number dropped because of missings. Without going back to R and > fiddling about to find the total number of rows in

I think that there is a bug in the as.POSIXct function on Windows. Here is what I get on Win2000, Pentium III machine in R 1.8.1. > dd1 <- ISOdatetime(2003, 10, 26, 0, 59, 59) > dd2 <- ISOdatetime(2003, 10, 26, 1, 0, 0) > dd2 - dd1 Time difference of 1.000278 hours Now, the 26th of October was the day that change to the standard time occurred, so I suspect that this has

Dear List, I have an array of 6400 x 56 elements. I want to calculate the density function for each column and plot all 56 density functions in one plot. I have tried several procedures, but they all failed. What can I do? TIA Christian [[alternative HTML version deleted]]

Having trouble with netCDF since I upgraded to 1.2.1 (did directly from 1.1.1). Package was recompiled after the upgrade. Symptoms: open.netCDF returns cleanly, but read.netCDF causes R to segfault and dump core. What else do you need me to write to help? --please do not edit the information below-- Version: platform = i586-pc-linux-gnu arch = i586 os = linux-gnu system = i586, linux-gnu

OK. What if I have a time series which is collected every Monday, please? What is the proper way to use the start option within the ts command in order to indicate that this is Monday data, please? Thanks again! Sincerely, Erin

Within an user function, how are the warnings from an R function be trapped (such that some proper actions can be taken)? 'last.warning' is returned only at the top level. Pointers are appreciated. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

I've been developing a package and have been getting the following warning when running the check command: * checking S3 generic/method consistency ... WARNING plot: function(x, ...) plot.summaries: function(trees, sp) * checking for replacement functions with final arg not named 'value' ... OK * checking Rd files ... OK ...blah, blah, blah... * checking examples ... OK *

Hi, I am building r-patch from the sources (rsync-ed today). make check produced the following message: running tests of Internet and socket functions expect some differences make[3]: Entering directory `/usr/evahome/vograno/R/tests' running code in 'internet.R' ... OK comparing 'internet.Rout' to './internet.Rout.save' ...18c18 < Content type `text/plain;

Hi, I am trying to use piecewise linear regression to approximate a nonlinear function. Actually, I don't know how many linear functions I need, therefore, I want build an array of regression models to automate the approximation job. Could you please give me any clue? Attached is ongoing code: rawData = scan("c:/zyang/mass/data/A01/1.PRN", what=list(numeric(),numeric())); len =

I am new to R. Trying to find out how lda() {in MASS R1.8.0 Windows} was implemented in R. Does anyone know where to find out lda source code ? Thanks. Wei

Dear R experts, I installed R-1.7.1 on Linux (Red Hat 9.0) starting from R-1.7.1.tgz without a problem. Then I fired up R and tried things and found that help.start( ) does not work. I downloaded Sun Java j2re1.4.2, installed that and re-installed R-1.7.1 from scratch. I tried help.start( ) again and the browser (Mozilla 1.2.1) crashed. I read about copying libjavaplugin_oji.so into

I would like to produce a set of values, on the same chart, from an sql table that is structured like... species dbh ht expf DF 1.2 8.9 10.0 DF 2.4 17.3 12.4342 DF 3.1 20.9 56.76 PP 2.3 16.9 100.0 PP 12.8 97.3 40.3 PP 8.2 63.0 98.34 . . . SS blah, blah, blah... is it possible to, using a single command in the plot command to plot the different groups on the same plot or will I have to iterate

Background: OS: Linux Mandrake 9.1 release: R 1.9.0 editor: Xemacs 21.4 frontend: ESS 5.1.23 --------------------------------- Colleagues I have two time series (upwelling index and water temperature) of evenly spaced, daily data over 18 months, but the upwelling index series has a gap of about 2 months right in the middle of it. I want to do the acf, pacf, ccf, and a cross-spectral analysis

I sent this to the list last week, and haven't seen it pop up. Either I deleted it when it did appear, or possibly it was destroyed as spam...? If it did appear and I somehow missed it, appologies. In a nutshell, can the function FUN supplied to by() deduce what level of factor by() was on when FUN was called? I've been digging through the functions, and can't see where the

Dear R experts, density estimation on a 2 dimensional bounded domain --------------------------------------------------------------------- I am currently trying to estimate the probability density (PD) of cancers within the breast using the sm library with the routine sm.density Of course a practical PD must be limited by the curve of the breast outline. I don't have a clue after perusing

I understand that the tsplot function is no longer available in the latest verion. What should I do if I want to plot say the first five rows of a matrix?? -gyan -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the