similar to: coupled statistical models

Dear R-listers, I have to explain "y", a continous variable, with two variables "t" and "p" which are "factor variables": t is [1:12] and p is [1:18] (a value for y for each combination t*p). OK, then I try to adjust models like: y=a_p + b*t (I mean as many intercepts as levels of p) or y=a_p + b_p*t (as many intercepts and coefficients for t as levels of

Dear R experts, I would greatly appreciate your help with R's hypertext online documentation, which I can pass on to our system manager. We are running version 1.7 on sparc-sun-solaris2.9 using netscape7 browser for displaying the help. The system is working, but there are a few frustrating details that I hope might be fixed. Here are three problems we are having: 1. When the R session

Sorry to bother you all with what is probably a simple question. I am attempting my first go at a simple AGI application using PHP (Getting Asterisk to SAY PHONETIC ABC). I have dabbled with PHP but I am by no means a professional standard developer. My script seems to execute ok, and I can see asterisk playing the sounds but my phone goes from ringing to busy, and I don't hear the

Dear R listers, When using smooth.spline to interpolate data, results are generally good. However, some cases produce totally unreasonable results. The data are values of pressure, temperature, and salinity from a probe that is lowered into the ocean, and the objective is to interpolate temperature and salinity to specified pressures. While smooth.spline provides excellent values at the

Sorry to re-post this but I'm sure it's something simple that someone has found before. To summarise: Dial plan answers the phone AGI script executes AGI debug in console show phonetics ABC - However no audio at all on the phone and this step is less than 1 second. Dial plan Busy Phone hangs up. Total time less than a second. This is such a simple AGI script do I need

The University of Miami's Rosenstiel School of Marine and Atmospheric Sciences (RSMAS) is seeking a research associate to work with hydrographic data. The objective of the project is to develop statistical models for estimating salinity and density from observations of temperature together with knowledge of location and time. Because the relationship between salinity and temperature can be

The University of Miami's Rosenstiel School of Marine and Atmospheric Sciences (RSMAS) is seeking a research associate to work with hydrographic data. The objective of the project is to develop statistical models for estimating salinity and density from observations of temperature together with knowledge of location and time. Because the relationship between salinity and temperature can be

Hello, Apologies if this question has already arised, hope you can help me to the find the solution to this or point the place to look at. I have a multidimensional piece-wise regression linear problem, i.e. to find not only the regression coefficients for each "interval" but also the beginning and ends of the intervals. To simplify it to the one dimensional case and two intervals,

...as promised. This describes the mdct used in my d.m.l patch. I think it is the same as the Lee fast-dct. I typed it in a kind of pseudo-TeX, 'cause the ascii art would kill me. Hope you can read TeX source; if not, ask someone who can to make a .ps/.gif/.whatever of the TeX output, and put it on a webpage or something. I'm to lazy to do it (and besides, I don't have access to TeX,

This posting has nothing to do with R (except maybe that I am using R very heavily in writing the paper to which the question pertains.) I simply wish to draw upon the impressive knowledge and wisdom of the R community. Since this question is way off topic, if anybody has the urge to reply, they should probably email me directly: rolf at math.unb.ca rather than via this list. My question

Sorry for starting another topic, this is actually a reply to Segher's post on Sun Oct 22 on the 'mdct question' topic. I wasn't subscribed properly and so I didn't get email confirmation and thus can't add to that thread. So Segher, if the equation is indeed what you say it is, then replacing mdct_backward with this version should work, but it doesn't. Am I applying

Hello, I was wondering if someone knows the formula used by the function lm to compute the t-values. I am trying to implement a linear regression myself. Assuming that I have K variables, and N observations, the formula I am using is: For the k-th variable, t-value= b_k/sigma_k With b_k is the coefficient for the k-th variable, and sigma_k =(t(x) x )^(-1) _kk is its standard deviation.

<head> <meta http-equiv="Content-Type" content="text/html; charset=windows-1252"> <meta name="GENERATOR" content="Microsoft FrontPage 4.0"> <meta name="ProgId" content="FrontPage.Editor.Document"> <title>Do you want to capitalize on the Biotech Revolution</title> </head> <body> <p

Dear all, I've only just joined this list, but I can't find any mention of this idea anywhere else, so I thought I'd just post here before getting too deep into programming and possibly reinventing the wheel. Here at Aber, we have around 30 unix and linux servers doing core services. Each one is maintaining its own logfiles and, for various reasons, we want to keep these on the

Is there a weighted median function out there similar to weighted.mean() but for medians? If not, I'll try implement or port it myself. The need for a weighted median came from the following optimization problem: x* = arg_x min (a|x| + sum_{k=1}^n |x - b_k|) where a : is a *positive* real scalar x : is a real scalar n : is an integer b_k: are negative and positive scalars

Hi R-helpers, I'm new to time series modelling, but my requirement seems to fall just outside the capabilities of the arima function in R. I'd like to fit an ARMA model where the variance of the disturbances is a function of some exogenous variable. So something like: Y_t = a_0 + a_1 * Y_(t-1) +...+ a_p * Y_(t-p) + b_1 * e_(t-1) +...+ b_q * e_(t-q) + e_t, where e_t ~ N(0, sigma^2_t),

Dear R experts, When using R 2.1 under Red Hat Enterprise Linux 4 with Mozilla 1.7.6 as browser, help.start() almost works correctly. When I follow one of the search result links everything is fine, but when I go back to follow another link, I find them all inactive. So I have to redo the search. Any suggestions? Thanks, Carlisle Thacker > version _ platform

Puppet 3.0.1 is a bugfix release for the 3.x series of Puppet. To see a list of the issues addressed by this release, check out the 3.0.1 version in our issue tracter at: https://projects.puppetlabs.com/versions/328. This release contains contributions from the following people: Andrew Parker, Jeff McCune, Josh Cooper, Matthaus Owens, Moses Mendoza, Nick Fagerlund, Niels Abspoel, Patrick

I am exploring the use of loess for oceanographic applications and would like to plot the locations (longitude and latitude) points where the models (salinity~temperature*longitude*latitude,parametric="temperature") are fitted. Chambers and Hastie(1993) explains the locations are nodes of a k-d tree. but I cannot find anything about accessing this information. It would be useful to

Dear readers, Is it possible to specify a model y=X %*% beta + Z %*% b ; b=(b_1,..,b_k) and b_i~N(0,v^2) for i=1,..,k that is, a model where the random slopes for different covariates are i.i.d., in lmer() and how? In lme() one needs a constant grouping factor (e.g.: all=rep(1,n)) and would then specify: lme(fixed= y~X, random= list(all=pdIdent(~Z-1)) ) , that?s how it's done in the