Hello all,
In order to find the kernel smoothing density function of a 3 dimensional data-set, I wrote this code in R: library(ks)defect <-read.table(file.choose(),sep=",")hat <-kde(defect)pdf("SampleGraph.pdf")plot(hat,drawpoints=TRUE)dev.off() But I have problem in interpretation of the result..I don't see the estimated values for density function (I need to have the values of f to use it in another code). The other problem is that the pdf file is not created correctly..Although the plot is created in R window, when I want to open the pdf file, it says:There was an error opening this document. This file cannot be opened because it has no pages. [[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list [hidden email] https://stat.ethz.ch/mailman/listinfo/r-sig-geo |
Hi,
On Thu, Jul 24, 2014 at 1:22 PM, Ferra Xu <[hidden email]> wrote: > Hello all, > In order to find the kernel smoothing density function of a 3 dimensional data-set, I wrote this code in R: > library(ks)defect <-read.table(file.choose(),sep=",")hat <-kde(defect)pdf("SampleGraph.pdf")plot(hat,drawpoints=TRUE)dev.off() > But I have problem in interpretation of the result..I don't see the estimated values for density function (I need to have the values of f to use it in another code). The other problem is that the pdf file is not created correctly..Although the plot is created in R window, when I want to open the pdf file, it says:There was an error opening this document. This file cannot be opened because it has no pages. > [[alternative HTML version deleted]] See why you shouldn't post HTML to this list? Trying to parse that, I see that I can't run your example because no data are provided (dput() is a good function to know). When I run the example code given in ?kde it works as expected. Does your plot command create a plot on other devices, like the default for your OS? Does the example code in ?kde work when you run it? Does it create a plot? As for the estimated values, did you carefully read the Value section of the help file? A kernel density estimate is an object of class ‘kde’ which is a list with fields: x: data points - same as input eval.points: points at which the estimate is evaluated estimate: density estimate at ‘eval.points’ If this density estimate is not the estimate you are looking for, you'll need to be more specific. Sarah -- Sarah Goslee http://www.functionaldiversity.org _______________________________________________ R-sig-Geo mailing list [hidden email] https://stat.ethz.ch/mailman/listinfo/r-sig-geo |
Dear Sarah Thank you for your response. I have attached the data that I'm using. I noticed that no estimate and eval.points are generated...Would you please take a look at it? It creates a plot in R environment but the pdf file doesn't work. Appreciate your help, On Thursday, July 24, 2014 12:35 PM, Sarah Goslee <[hidden email]> wrote: Hi, On Thu, Jul 24, 2014 at 1:22 PM, Ferra Xu <[hidden email]> wrote: > Hello all, > In order to find the kernel smoothing density function of a 3 dimensional data-set, I wrote this code in R: > library(ks)defect <-read.table(file.choose(),sep=",")hat <-kde(defect)pdf("SampleGraph.pdf")plot(hat,drawpoints=TRUE)dev.off() > But I have problem in interpretation of the result..I don't see the estimated values for density function (I need to have the values of f to use it in another code). The other problem is that the pdf file is not created correctly..Although the plot is created in R window, when I want to open the pdf file, it says:There was an error opening this document. This file cannot be opened because it has no pages. > [[alternative HTML version deleted]] See why you shouldn't post HTML to this list? Trying to parse that, I see that I can't run your example because no data are provided (dput() is a good function to know). When I run the example code given in ?kde it works as expected. Does your plot command create a plot on other devices, like the default for your OS? Does the example code in ?kde work when you run it? Does it create a plot? As for the estimated values, did you carefully read the Value section of the help file? A kernel density estimate is an object of class ‘kde’ which is a list with fields: x: data points - same as input eval.points: points at which the estimate is evaluated estimate: density estimate at ‘eval.points’ If this density estimate is not the estimate you are looking for, you'll need to be more specific. Sarah -- Sarah Goslee http://www.functionaldiversity.org _______________________________________________ R-sig-Geo mailing list [hidden email] https://stat.ethz.ch/mailman/listinfo/r-sig-geo kernel1.csv (13K) Download Attachment |
On Thu, Jul 24, 2014 at 4:06 PM, Ferra Xu <[hidden email]> wrote:
> Dear Sarah > Thank you for your response. I have attached the data that I'm using. I > noticed that no estimate and eval.points are generated...Would you please > take a look at it? Of course the estimate is generated, just as the help file promised. hat$estimate contains the estimate: all the named components given in the Value section of the helpfile are part of hat, and you can see that with names(hat) Have you read the Introduction to R that comes with your R installation? Perhaps you should. > It creates a plot in R environment but the pdf file doesn't work. > > Appreciate your help, Well, ?plot.kde states that for 3D data it produces an interactive RGL plot, and unlike for 1D and 2D there are no other options. What do you expect to be printed to PDF when you create an interactive plot? Reading ?rgl suggests that you might want to look into ?rgl.snapshot Sarah > > On Thursday, July 24, 2014 12:35 PM, Sarah Goslee <[hidden email]> > wrote: > > > Hi, > > On Thu, Jul 24, 2014 at 1:22 PM, Ferra Xu <[hidden email]> wrote: >> Hello all, >> In order to find the kernel smoothing density function of a 3 dimensional >> data-set, I wrote this code in R: >> library(ks)defect <-read.table(file.choose(),sep=",")hat >> <-kde(defect)pdf("SampleGraph.pdf")plot(hat,drawpoints=TRUE)dev.off() >> But I have problem in interpretation of the result..I don't see the >> estimated values for density function (I need to have the values of f to use >> it in another code). The other problem is that the pdf file is not created >> correctly..Although the plot is created in R window, when I want to open the >> pdf file, it says:There was an error opening this document. This file cannot >> be opened because it has no pages. > >> [[alternative HTML version deleted]] > > > See why you shouldn't post HTML to this list? > > Trying to parse that, I see that I can't run your example because no > data are provided (dput() is a good function to know). > > When I run the example code given in ?kde it works as expected. Does > your plot command create a plot on other devices, like the default for > your OS? Does the example code in ?kde work when you run it? Does it > create a plot? > > As for the estimated values, did you carefully read the Value section > of the help file? > > A kernel density estimate is an object of class ‘kde’ which is a > list with fields: > > x: data points - same as input > > eval.points: points at which the estimate is evaluated > > estimate: density estimate at ‘eval.points’ > > If this density estimate is not the estimate you are looking for, > you'll need to be more specific. > > Sarah > -- Sarah Goslee http://www.functionaldiversity.org _______________________________________________ R-sig-Geo mailing list [hidden email] https://stat.ethz.ch/mailman/listinfo/r-sig-geo |
Free forum by Nabble | Edit this page |