[R] bivariate weighted kernel density estimator

Is there code for bivariate kernel density estimation? For bivariate kernels there is kde2d in MASS kde2d.g in GRASS KernSur in GenKern (list probably incomplete) but none of them seems to accept a weight parameter (like density does since R 2.2.0) -- Erich Neuwirth, University of Vienna Faculty

Re: [R] second try; writing user-defined GLM link function

Ok, I think I have things figured out. Predict.glm still doesn't work (with type=response); however, the problem with the strange behavior of the betas and the results of predict.glm with type=link was simply due to data issues. I had converted my Program MARK input file (which needs only three

Re: [R] Monitor() in R2WinBUGS

Doran, Harold wrote: I have an array from an mcmc simulation and I'd like to try and use the monitor function in R2WinBUGS. The R2WinBUGS package is in the search path and when I try help(monitor) the help page appears for the function. But, when I try and use it the following happens:

[R] help! A quetion about the Elasticnet package in R

Can anybody help me to run the Elasticnet package of R to build some model, i am a freshman to R language , when i use the Elasticnet package to my data, it always reture a error, but i can't settle that problem. I consider if there is any constrant of the data to that package? Can anyone help

Re: [R] how to draw a circle

You may want to have a look at the package pgirmess and the function polycirc() Kind regards, Patrick __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

[R] (no subject)

Ilija Barukcic, Jever, Germany. GMT + 1h. [EMAIL PROTECTED] Hallo, you are receiving this message from me as someone who potentially has bought the first Edition of my book Causality. New statistical methods. ISBN 3 - 8 3 3 4 - 3 6 4 5 - X. Your email address was obtained either through my

[R] Comparing GLMMs and GLMs with quasi-binomial errors?

Dear All, I am analysing a dataset on levels of herbivory in seedlings in an experimental setup in a rainforest. I have seven classes/categories of seedling damage/herbivory that I want to analyse, modelling each separately. There are twenty maternal trees, with eight groups of seedlings

Re: [R] bivariate weighted kernel density estimator [Broadcast]

locfit() in the locfit package can do that. Andy _ From: [EMAIL PROTECTED] on behalf of Erich Neuwirth Sent: Sun 4/23/2006 3:13 AM To: r-help Subject: [R] bivariate weighted kernel density estimator [Broadcast] Is there code for bivariate kernel density estimation? For bivariate

Re: [R] Extract AIC, BIC

Which Garch function are you using? help.search(garch) identified garch functions in fOptions, fSeries, and tseries. If you try RSiteSearch(garch, functions), you might find others. I just tried 'library(tseries)' and tried an example in help(garch). I modified this

[R] Question about bicreg

Dear Adrian and Ian (and r-helpers), I encountered a curious result in developing an example using the bicreg function in the BMA package: I noticed that pairs of models with equal R^2 and equal numbers of predictors had nevertheless different BIC values. Looking at the bicreg function, the

Re: [R] help! A quetion about the Elasticnet package in R

邓敏 [EMAIL PROTECTED] writes: Can anybody help me to run the Elasticnet package of R to build some model, i am a freshman to R language , when i use the Elasticnet package to my data, it always reture a error, but i can't settle that problem. I consider if there is any constrant of the data

[R] RGL and scales of axes

Hi, Is there anyway i can set the aspect ratio options in an RGL graphic? Like in lattice. panel.cloud has an option called aspect=c(r1,r2) where r1 is the y/x ratio and r2 is z/x ratio. I'm not sure if this is the right place to ask this question. Thanks for

Re: [R] info : Manova - eigenvector analysis and canonical analysis

Have you considered the sem package (for structural equation modeling), which (to me at least) is a generalizion of MANOVA with canonical analysis. Alternatively, have you considered partial least squares (e.g., packages pls or plsgenomics)? I haven't used them, but they sound like

[R] Cohen-Daubechies-Fauraue (CDF) 9/7 wavelet-transform

Hello, I have to make use of the Cohen-Daubechies-Fauraue (CDF) 9/7 wavelet-transform, but have not found these biorthogonal wavlet in the released packages waveslim, wavethresh or wavelets. Is there any package doing this transform? Or can anybody send me the function doing this? I would be

[R] lme: null deviance, deviance due to the random effects, residual deviance

A maybe trivial and stupid question: In the case of a lm or glm fit, it is quite informative (to me) to have a look to the null deviance and the residual deviance of a model. This is generally provided in the print method or the summary, eg: Null Deviance: 658.8 Residual Deviance: 507.3

Re: [R] help! A quetion about the Elasticnet package in R

Dear Leon: How good is your Manderin? I am pleased and humbled that Mr. Dengmin has made the effort to put his question in English. I don't know his background, but my Eastern-language enabled email provided a rendition of Richard Dengmin in the form or two Chinese characters,

Re: [R] help! A quetion about the Elasticnet package in R

Spencer Graves [EMAIL PROTECTED] writes: Dear Leon: How good is your Manderin? I am pleased and humbled that Mr. Dengmin has made the effort to put his question in English. I don't know his background, but my Eastern-language enabled email provided a rendition of Richard

[R] Kullback Liebler

Dear R Users, Does anyone have an implementation of KL apart from what is in reldist ? Something with flexibility to take in empirical data for two sets, or their CDF or their density, flagged separately for each ? Thanks __ R-help@stat.math.ethz.ch

Re: [R] help! A quetion about the Elasticnet package in R

Dear Leon: Please forgive me. As I said, I'm pleased and humbled that so many others learn English. I reacted as I did because I've heard too many people (fortunately NOT on this list) criticize or make fun of others whose use of language was obviously not native. Best

Re: [R] help! A quetion about the Elasticnet package in R

Dear Spencer, No worries:-) -- Leon Spencer Graves [EMAIL PROTECTED] writes: Dear Leon: Please forgive me. As I said, I'm pleased and humbled that so many others learn English. I reacted as I did because I've heard too many people (fortunately NOT on this list) criticize or

[R] distribution of the product of two correlated normal

Hi, Does anyone know what the distribution for the product of two correlated normal? Say I have X~N(a, \sigma1^2) and Y~N(b, \sigma2^2), and the \rou(X,Y) is not equal to 0, I want to know the pdf or cdf of XY. Thanks a lot in advance. yu [[alternative HTML version deleted]]

[R] Reorganizing rows and columns

I'm sure this is a simple task, but how to do it has escaped me. I have imported data from two separate files (each file contains the results from an information retrieval algorithm) organized into a list. They are organized by File,Query, and Rank (in that order): [[1]] Doc Query Rank 5

Re: [R] Reorganizing rows and columns

Use merge: # test data both - list(structure(list(Doc = c(5, 9, 7, 5, 7, 9), Query = c(1, 1, 1, 2, 2, 2), Rank = c(1, 2, 3, 1, 2, 3)), .Names = c(Doc, Query, Rank), class = data.frame, row.names = c(1, 2, 3, 4, 5, 6)), structure(list(Doc = c(4, 5, 9, 8, 5, 7), Query = c(1, 1, 1, 2, 2, 2), Rank =

[R] question about cumulating random effects in lmer

Hello: I am studying the effect of schools on student achievement growth over time. School effect is random. The effects of schools in prior years are assumed to be persistent till the current year. Thus, the total school effect in the second year is like J=J2+J1. Does anyone know how to model

Re: [R] Question about bicreg

Dear Adrian and Ian (and r-helpers), To answer my own question: The difference is simply due to rounding; that is, there are actually small differences in the R^2 values that account for the differences in the BICs. Sorry for the false alarm, John -Original Message- From: John Fox

[R] bivariate weighted kernel density estimator

Erich Neuwirth writes: Is there code for bivariate kernel density estimation? ... but none of them seems to accept a weight parameter In package 'spatstat' the function density.ppp performs weighted kernel smoothing, including bivariate kernel density estimation. At the moment it only

[R] garch warning

Dear r users, Few days ago I posted the same topic but unable to receive any suggestion. So I am asking this same question. I was trying to fit a garch(1,1) model to my dataset. But while executing I got a warning message NaNs produced in: sqrt(pred$e). And got the estimated sd's

Re: [R] Reorganizing rows and columns

Suppose your data frames are A and B: AB - merge(A, B, c(Doc, Query), all=TRUE) AB[order(AB$Query, AB$Doc),] gets the answer you are asking for. (Not sure why you want to sort it to use a correlation test, as those are indifferent to ordering.) On Sun, 23 Apr 2006, [EMAIL PROTECTED] wrote:

Re: [R] Comparing GLMMs and GLMs with quasi-binomial errors?

AIC is only valid for maximum likelihood fitting, so not for PQL and also not for quasi-binomial. In particular, who said AIC() was valid when applied to a glmmPQL fit? (Certainly not the book it supports!) You say TreeID is a `blocking factor', but you have a treatment (PLATEAU) that is