similar to: Newbie question about plotting density objects

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:

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 am using an older version of R (1.6.2) to run a Monte Carlo simulation, generating 10,000 samples per 'run'. When I plot histograms I get the expected 'bins' on the x-axis and the frequency distribution on the y-axis. However when I ask R to plot the SAME data set with a density curve the x-axis remains the same but the y-axis can generate values of up to 1e8 etc. Can anyone

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 =

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

How do you draw two histograms in one figure? Two separate uses of the "hist" function always produce two separate figures... Thanks, Roman -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

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

Hello all, I wonder if there is a package including a program to display a 3dim histogram. Thanks. --- D.Trenkler --- **************************************************************************** ***** Dr. Dietrich Trenkler (dtrenkler at nts6.oec.uni-osnabrueck.de) Statistik / Empirische Wirtschaftsforschung Universitaet Osnabrueck Rolandstrasse 8

Hi R lovers! Are there any means to define and plot a function given the equation that specifies the function? For example I'd like to plot and work with the Gumbel Distribution density defined by Lambda(x)=exp(-exp(-x)) My question may appear very simple but I haven't got an idea yet about how to do that. I could plot something with x a vector/set of value but I don't know how to

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

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;

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 wish to produce a weighted (Gaussian?) kernel density estimate, in 2 dimensions, where the weights are permitted to be ***negative***. (I can ASSURE you that there are perfectly legitimate reasons why I want to do this. :-)) Clearly it is not really a density that I am trying to estimate. I had a go at using sm.density() from the sm package, but this requires that the weights be positive. Is