search for: parametrised

Dear all, I do not seem to grasp how contrasts are set for ordered factors. Perhaps someone can elighten me? When I work with ordered factors, I would often like to be able to reduce the used polynomial to a simpler one (where possible). Thus, I would like to explicetly code the polynomial but ideally, the intial model (thus, the full polynomial) would be identical to one with an ordered factor.

...R-ish in that the results depend on the parametrisation of the model which by default in R is formed using the so-called "treatment" contrasts. To wander from R into statistics (sorry Bert) the problem arises because the "usual" parametrisation of the model is the "over-parametrised" form: Y_ij = mu + beta_i + E_ij (i = 1, ..., I, j = 1, ..., J_i) where Y_ij is the j-th observation corresponding to the i-th "treatment" or group. (Things get a bit more complicated in "multi-way" models; let's not go there.) The parameter "mu" is the...

This is really a statistical issue. What do you think the Intercept term represents? See ?contrasts. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Thu, Feb 15, 2018 at 5:27 PM, Marc Girondot via R-help < r-help at

Dear R-er, I try to get the standard error of fitted parameters for factors with a glm, even the reference one: a <- runif(100) b <- sample(x=c("0", "1", "2"), size=100, replace = TRUE) df <- data.frame(A=a, B=b, stringsAsFactors = FALSE) g <- glm(a ~ b, data=df) summary(g)$coefficients # I don't get SE for the reference factor, here 0:

Hi list, Reading the thread "can a class require an other class?" it''s been mentioned that perhaps one way forward for the Puppet language is to phase out the include keyword in favour of parametrised classes. I''m thinking of my Puppet and the several levels of include chaining I use and I''m wondering how on earth that''d be possible. Maybe I''m designing my modules very differently to other people? I''ll throw you an edited example from my environment...

Dear listers, I have written a function to facilitate the drawing of altitude profiles with x (distance), y (altitude) and a z parameter (altitude magnification). profplot<-function(x,y,z=10,...){ op <- par()$mai par(mai=c(0.95625,0.76875,0.76875,0.95625)) plot(x,y*z, type="l",asp=1,las=1,xlab="",ylab="",yaxt="n",...)

In Venables \& Ripley 3rd edition (p. 231) the proportional odds model is described as: logit(p<=k) = zeta_k + eta but polr apparently thinks there is a minus in front of eta, as is apprent below. Is this a bug og a feature I have overlooked? Here is the naked code for reproduction, below the results. ------------------------------------------------------------------------ --- version

...e* same puppet environment,* so I cannot use %{environment} I have a module that sets all the *basic* functionality for the OS, resolution, authentication, security, packages, etc I have a module for each application hosted. At the moment all the ''data'' is in Puppet, mostly in parametrised classes in site.pp. What I want to get is a hiera set-up that allows me to use this structure: :hierarchy: - global # source of application names (classes? modules?) and environments list - basic # data for the basic class - prod/%{application}/%{hostname} # hostname files...

> I am planning in doing in IR, but with target specific-passes (such as X86ExpandAtomicPass) > that just share some of the code This would more normally be done via target hooks in LLVM, though the principle is sound. > But it must be target-dependent as for example on Power a > seq_cst store has a fence before it, while on ARM it has a fence > both before and after it (per

Dear Sir/Madam: this is shuangge Ma, graduate student in UW-Madison statistics department. I have a computational question. I have a function f(x,y). I want to find the y(x) that maximize f(x,y) under the constraint y(x) is a non-decreasing step function. Is there any R package or algorithm I can use for this purpose? thanks a lot for your time and help, Sincerely, Shuangge Ma

I am learning glm function, but how do you fit data using exponential distribution with glm? In the help file, under "Family Objects for Models", no ready made option seems available for the distribution as well as for other distributions satisfying GLM requirements not listed there.

...very interested to know if this can be done too. I'm thinking of something which would work like the 'define' mechanism in C (and some other languages), where #define nis.na !is.na would cause simple substitution of "!is.na" for "nis.na" on input (and you can have parametrised defines too). With thanks, Ted. -------------------------------------------------------------------- E-Mail: (Ted Harding) <Ted.Harding at nessie.mcc.ac.uk> Fax-to-email: +44 (0)870 167 1972 Date: 20-Oct-03 Time: 20:38:01 ------------------------------ XFMail -----------------------------...

...ith special correlation structures, or through random effects in lme/lmer ? should I try to define a "new" correlation structure corMultiAR1 ? If so, where can I find help on how to write such a piece of code ( nlme:::corAR1 is not clear to me) ? Or is there a way to define a general parametrised covariance matrix in gls ? Olivier > obs6 <- matrix( c(1,2,3,4,5,6, 1,1,2,2,3,3), byrow=F, nc=2) > dimnames(obs6) <- list(NULL, c("y","time")) > obs6 <- data.frame(obs6) > obs6 y time 1 1 1 2 2 1 3 3 2 4 4 2 5 5 3 6 6 3 > gls (y~1,...

Hi All, I'm trying to use one of the (2D) numerical integration functions, which is not where the problem is. The function definition is as follows: adaptIntegrate(f, lowerLimit, upperLimit, ...) The problem is that I want to integrate a 3D function which has been parametrised such that it is a 2D function: eg. I want to integrate f(x_start+gradient_x*t1, y_start+gradient_y*t2) for t1 in [0,1] and t2 in [0,1]. To explain the example in another way I only have access to f, x_start, y_start and the gradients. Note that I want to be able to change x_start, y_start such tha...

Hi, Is it possible to redirect passing multi dimensional parameters? I need to send params[:current][:selview] to another action mantaining the same data structure if I try as follows: redirect_to :controller=>"main",:action=>"inquiry",{"current"=>{ "selview"=>params[:current][:selview]}} I receive "syntax error, unexpected

Dear all, apologies for interrupting your important work. The new R package MittagLeffleR is now available on CRAN. It computes the two types of Mittag-Leffler distributions, i.e. provides probability density, distribution function, quantile function and random variate generation for the Mittag-Leffler distributions, and the Mittag-Leffler function itself. It is based on the Laplace-inversion

Dear list, I've generated a list of 3D coordinates representing ellipsoids in arbitrary orientations. I'm now trying to obtain a 2D projection of the scene, that is to draw the silhouette of each object on a plane (x,y). The only way I could think of is to compute the convex hull of the (x,y) coordinates of each object and use this as the outline of the object. This is clearly

A non-text attachment was scrubbed... Name: not available Type: text Size: 1485 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980904/6c03ac2b/attachment.pl

Hello! I would like to use constrain to sum coeficients of a factor to 0 instead of classical corner contraint i.e. I would like to fit a model like lm(y ~ 1 + effectA + effectB) and say get parameters intercept effectA_1 effectA_2 effectB_1 effectB_2 effectB_3 where effectA_1 represents deviation of level A_1 from intercept and sum(effectA_1, effectA_2) = 0 and the same for factor B. Is

Hello, I'm a grad student in the Intelligent Transportation Systems lab at Portland State Univ. in Portland, OR, USA. I'm trying to learn the basics of R to run a negative binomial in the near future, and so I ran a test regression on roadway crash data obtained from "Statistical and Econometric Methods for Transportation Data Analysis" by Washington et al (p. 250). I ran the