similar to: Fitting a Harmonic Function to Time Series Data

Hi Folks, I think I have searched exhaustively, including, of course R-help (D. Bates, S. Graves, and others) and but I remain uncertain about testing fixed effects with lmer(..., family=binomial). I gather that mcmcsamp does not work with Do we rely exclusively on z values of model parameters, or could we use anova() with likelihood ratios, AIC and BIC, with (or without)

Hi folks, Warning: I don't know if the result I am getting makes sense, so this may be a statistics question. The fitted values from my binomial lmer mixed model seem to consistently overestimate the cell means, and I don't know why. I assume I am doing something stupid. Below I include code, and a binary image of the data is available at this link:

Hi all, I am trying to re-compile some "unmaintained" (it seems) packages, namely rSoNIA and dynamicnetwork from: http://csde.washington.edu/~skyebend/rsonia/rsoniaDemo/ These packages predates R 2.10.0 so they need to be recompile. After split the single big file in /man in each packages into a file for each function + some minor fix, I successfully manage to recompile and load

I am looking for a mailing list for general statistical questions that are not R related. Do you have any suggestions for lists that are busy and helpful and/or lists that you use and recommend? Thanks in advance, Ralf

Hi everyone, I am trying to construct a glm and am running into a couple of questions. The data set I am using consists of 6 categories for the response and 6 independent predictors representing nutrient concentrations at sample point locations. Ultimately I'd like to use the probabilities for each response category in a simulation model such that these probabilities are used to define a

Hello. I'm running a generalized linear model and am interested in using the random effects that are output for further analysis. My random effect is interacting with two different fixed effects (which which are factors with two levels each). When I retrieve the random effects I get something like this: (Intercept) nutrient (Intercept) light (Intercept) Aa-0 0.59679192

Hello, I will try to explain the problem, sorry if it will be a little long... I'm using R to analyze results of cyclic mechanical testing, like this: - apply quasi-sinusoidal load - measure quasi-sinusoidal vertical and horizontal deformations (quasi-sinusoidal load means that load "should be" sinusoidal, but testing machine puts in some noise...) I enclose a sample of data at

Hello, I try to synthesize harmonic functions with a subset of frequencies determined by fast fourier transform. I wrote two different functions, a looped version and a matrix multiplication version. As an example the looped version takes 5 sec and the matrix-version takes 3-4 sec (R 1.4, Athlon 1.2 GHz, 256 MB, Win ME), but the latter needs huge amount of memory due to the matrices. So, none of

Dear R Users, I have a CO2 time series. I want to fit this series seasonal cycle and trend with fourth harmonic function, and then compute residuals. I am doing something like: file<-read.csv("co2data.csv") names(file) attach(file) fit<-lm(co2~1+time+I(time^2)+sin(2*pi*time)+cos(2*pi*time)+sin(4*pi*time)+cos(4*pi*time)+

I am looking for a package to perform harmonic analysis with the goal of estimating the period of the dominant high frequency component in some mono-channel signals. I guess there are presumably a number of CRAN packages allowing for such analysis. However, my search with keywords was not successfull. It brought up a lot of Fourier miscellanea but nothing specifically geared for my needs. I would

Dear all: Does anyone has harmonic regresssion analysis package written in R (to be used in Windows platform) ? Thanks _______________________________________ YM - Â÷½u°T®§ ´Nºâ§A¨S¦³¤Wºô¡A§AªºªB¤Í¤´¥i¥H¯d¤U°T®§µ¹§A¡A·í§A¤Wºô®É´N¯à¥ß§Y¬Ý¨ì¡A¥ô¦ó»¡¸Ü³£ÉN¨«¥¢¡C http://messenger.yahoo.com.hk [[alternative HTML version deleted]]

Hi, Sorry for the very short explanation about the problem of intersection. I have a wave function monitored from the heart beat in a particular interval of times.?Apart?fom that, there is?a line with positive slope (e.g: y = x+2) which lies across the wave and intersect on a number of points. My problem is i have no exact equation for?such a?complex harmonic wave produced by the heart

Hello list; I am working with a dataset containing animal locations over time for a large number of individuals. I would like to compute a center of activity by finding the minimum harmonic mean on a grid overlaying the points of an animal (a "standard" to express center of activity in animal ecology). I have searched the archives at Jonathan Barons server but have not been able

I' need to smooth a serie univariate with harmonic regression, but i do not know if it exists somepeople know in R this type of regression, not in S-plus. Thank's _________________________________________________________

Dear R-Experts. I collected different datas about Nitrogen content (mg/ml) in Dung. The dung was eighter fresh (day=0) or had different ages (15,29,47) to observe nutrient changes over time. Now I like to draw a boxplot. boxplot(nmgml~day) abline((nmgml~day) The Problem is, that the boxplot considers the day values as groups and not as time series (neighter when the days are numeric or

I don't know if there is a way to use the scale relation free argument in ggplot2 like in lattice. I have a feeling that there is not, but I would like to make a plea for this feature. It would be nice to be able to plot Total Inorganic Nitrogen Total Phosphorus and the ratio of the two- the numbers on the axis are not related, but the previous two are surely related to the last (this ratio

Hi, Is there a package to perform a sine function fitting to XY data? Thx, Ashz -- View this message in context: http://r.789695.n4.nabble.com/Sine-function-fitting-tp3000156p3000156.html Sent from the R help mailing list archive at Nabble.com.

Hey, all. I'm working on a data set with a broken stick linear regression where I know one of the two slopes. It is a negative linear function until the line intersects with the x-axis, at which point it becomes 0. It is not a nonlinear asymptotic function, and, indeed, using negative exponential or logistic types of fits as an approximation has tended to lead to an under or

Hi, I am interested in ggplot2 and I found this lattice code very interesting (http://addictedtor.free.fr/graphiques/graphcode.php?graph=48). Code: library(lattice) lattice.options(default.theme = canonical.theme(color = FALSE)) tmp <- expand.grid(geology = c("Sand","Clay","Silt","Rock"), species = c("ArisDiff",

Hi, I have 20*60 data matrix (with some NAs) and I wish to perfom a Pearson correlation coefficient matrix as well as simple linear regression equation and coefficient of determination (R2) for every possible combination. Any tip/idea/library/script how do to so. Thanks, As hz -- View this message in context: