search for: vonmises

Hi, I'm dealing with wind data and I'd like to model their distribution in order to simulate data to fill-in missing values. Wind direction are typically following a vonmises distribution and wind speeds follow a weibull distribution. I'd like to build a joint distribution of directions and speeds as a VonMises-Weibull bivariate distribution. First is this a stupid question? I'm a newbie in statistics and R :) Is it possible to do it in R? Is there a way to...

...do this weekend to understand every bit but your code will prove very useful. Cheers, Aziz -----Original Message----- From: Dimitrios Rizopoulos [mailto:Dimitris.Rizopoulos at med.kuleuven.be] Sent: May 12, 2006 4:35 PM To: Chaouch, Aziz Subject: RE: [R] Maximum likelihood estimate of bivariate vonmises-weibulldistribution look at the following code: library(copula) par(mfrow = c(2, 2)) x <- mvdc(normalCopula(sin(0.5 * pi /2)), c("norm", "norm"), list(list(mean = 0, sd = 1), list(mean = 0, sd = 1))) contour(x, dmvdc, xlim = c(-2.7, 2.7), ylim = c(-2.7, 2.7)) x <- mvdc(...

Dear r-helpers, I want to show that time is flowing CCW in the following: require(circular) len <- 8 labl <- as.character(c(0, 1, 1, 1, 0, 0, 1, 0)) r <- circular(2*pi* (rep(c(1, 3, 6), each = 200)/len + rnorm(600, 0, 0.025))) r.dens <- density(r, bw = 25, adjust = 4, kernel = 'vonmises') plot(r, shrink = 2.5, axes = FALSE, ticks = FALSE, pch = 1, col = 'lightblue', stack = TRUE, bins = 12 * len) axis.circular(at = circular(seq(0, (len - 1) * 2 * pi/len, length.out = len)), label = labl) lines(r.dens, col = 2) I had imagined a directed arc with a smaller radius th...

...finds an almost flat plateau. How do I then test against the H0 - by definition, I'll never get a Kappa = 0 exactly; so I can't compare against that. One idea I had: Define not only a parameter Kappa, but also one of an indicator function, which acts as switch between a uniform and a vonMises distribution. Call that parameter d. I could then for example let d switch state with a 50% probability and then make usual acceptance tests. Is this approach realistic ? is it sound and solid or nonsense / suboptimal? Is there a common solution to the before mentioned problem ? [I suppose there...

...Function tobit Tobit Model trplot.qrrvglm Trajectory plot for QRR-VGLMs uqo Fitting Unconstrained Quadratic Ordination vgam Fitting Vector Generalized Additive Models vglm Fitting Vector Generalized Linear Models vonmises von Mises Distribution Family Function vsmooth.spline Vector cubic smoothing spline wald Wald Distribution Family Function weibull Weibull Distribution Family Function weightsvglm Prior and Working Weights of a VGLM fit yeo.john...