search for: correponding

Hi! I need to run ols regressions with Huber-White sandwich estimators and the correponding standard errors, without an intercept. What I'm trying to do is create an ols object and then use the robcov() function, on the order of: f <- ols(depvar ~ ind1 + ind2, x=TRUE) robcov(f) However, when I go f <- ols(depvar ~ ind1 + ind2 -1, x=TRUE) I get the following error: Error i...

Dear friends, In S-PLUS, we can use the following argument, but not in R. mode12 <- gam(score1 ~ lo(latitude) + lo(longitude)) I searched the help in S-PLUS, it says lo() Allows the user to specify a Loess fit in a GAM formula, but i didn't find the correponding function in R. Anybody knows how to do the similar task in R? Thanks very much. -- With Kind Regards, oooO::::::::: (..)::::::::: :\.(:::Oooo:: ::\_)::(..):: :::::::)./::: ::::::(_/:::: ::::::::::::: [***********************************************************************] Zhi Jie,Zhang ,PHD T...

...[8,] 3.9111111 105.35119 [9,] 4.4555556 105.50052 [10,] 5.0000000 111.35761 And I'd like to display the number in the first column, where its corresponding number is the minimum of the second column. In other words, I'd like to be able to extract out 1.1888889 from column 1, as its correponding number in column 2 (95.00240) is the smallest in the column. Cheers, Kevin ------------------------------------------------------------------------------ Ko-Kang Kevin Wang Postgraduate PGDipSci Student Department of Statistics University of Auckland New Zealand Homepage: http://www.stat.auckl...

Dear List, Does anyone know how to simulate data from a GLM object correponding to values of the independent (x) variable that do not occur in the original dataset? I have tried using simulate(), but it generates a new value of the dependent variable corresponding to each of the original x-values, which is not what I need. Ideally I whould like to simulate new values...

...(out1, function(m)summary(m)$coefficients[,4][2]) p1 <- do.call(rbind, p1) My ultimate objective is to apply this function to all the dataframes created (i.e. df1, df2, df3, df4, df5) and create five corresponding p-value vectors (p1, p2, p3, p4, p5). Then output would be a matrix of clvar and correponding p values clvar var1 var2 var3 var4 var5 1 2 3 4 Please help me ! Thanks NIL [[alternative HTML version deleted]]

...as.dist(a), method = "ward", members = NULL)), dist.method = NULL), .Names = c("merge", "height", "order", "labels", "method", "call", "dist.method"), class = "hclust") I just need the every item with its correponding cluster in a more or less organizated way. Of course, there is not problem in using different funtcions or librarys (till now I have not found anything sweeting to my needs). Advices or orientations are welcome and appreciated! Thanks in advance. AJ [[alternative HTML version delete...

...ge is that the two configurations are exactly simmetric, but only one of the two server fail to call the internal phone of the simmetric asterisk server. What is more ever strange is that if I replace the authenticate mode with md5 in [sintel-user] and [mickymouse-peer] sections below, (and the correponding dial command), all works fine! Here is an extract from the server A (mickymouse) iax.conf: ... [sintel-peer] type=peer host=192.168.100.1 auth=rsa qualify=yes [sintel-user] type=user auth=rsa inkeys=sintel-voip context=sintel-int qualify=yes ... extensions.conf: [voip_outcoming] exten => _2...

Hello people, I''m using squeegy fleximage for my site right now and having some difficulty to achieve some things. Updates photos are saved in created at date generated fom like : '' /images/2008/05/14/22.png '' . However this route is not accessible through the browser and I don''t know how to add it in my routes.rb knowing that the dates will change everyday.

Thanks for your advice. I am not sure that I understood your comment "If you need something, there should be a __builtin that corresponds to the intrinsic." Is that a better way to define an intrinsic function in C? How do you do it? I am actually trying to add several intrinsic functions for my target machine so I am looking for a simple and workable way of doing it. Thanks again.

Hi I'm using the "image" function to produce a plot and I want to define the breaks using "cut" and the colors using "heat.colors". > image(interp(mat2[,2],mat2[,1],mat2[,3]),breaks=cut(mat2[,3],30),col=heat.colors(29)) Error in image.default(interp(mat2[, 2], mat2[, 1], mat2[, 3]), breaks = cut(mat2[, : must have one more break than colour The

Greetings, I am new to R and have been comparing CPH survival analysis hazard ratios between R and SAS PhReg. The binary covariates' HRs are the same, however the continuous variables, for example age, have quite different HRs although in the same direction. SAS PhReg produces HRs which are the change in risk for every one increment change in the independent variable. How do I

I am trying to model a clusterd ordinal response data (either 1, 2 or 3) called , the correponding physician of the patient is also in the data. Since it is ordinal, I used the ordinal logit model topbox[i]~discrete with probability P[j,1],p[j,2], p[j,3], j is the corresponding physician of the ith patient C[j] is the physician effect , a1 and a1+theta is the common cutpoints for all physi...

...at yahoo.com.cn> > Subject: an ordinal regression MCMC run high correlation > To: r-help at r-project.org > Cc: chen1984612 at yahoo.com.cn > Date: Tuesday, March 16, 2010, 12:19 PM > I am trying to model a clusterd > ordinal response data (either 1, 2 or 3) called , the > correponding physician of the patient is also in the data. > > Since it is ordinal, I used the ordinal logit model > > topbox[i]~discrete with probability P[j,1],p[j,2], p[j,3], > j is the corresponding physician of the ith patient > > C[j] is the physician effect , a1 and a1+theta is...

Hello! How can one select outgoing MSN when dialing out from ttyI-interfaces? I have successfully done this with CAPI e.g... exten => _XXXXX.,2,Dial,CAPI/60:bBYEXTENSION ...in extensions.conf. Currently correponding for my ISDN modem interface is... exten => _XXXXX.,2,Dial(Modem/g1:${EXTEN}) ...but this selects only MSN of outgoing group g1 for dialout MSN number. I also tried to create multiple groups wait different MSNs and same devices in modem.conf, but each device can be only in one group. Is there...

...edeed for you to undertand my problem. Straight to the point: i have this domain in my company running in Samba 3.0.2 My users are: hcoelho, jardim, gamito, yesenia, smatias, fqueiros, faugusto, vamaro, peixinho, aragao, dina, pinho. I have this shares with the users that can access them and the correponding Linux groups: [DAT]: hcoelho, jardim, fqueiros, gamito, faugusto => Linux group A [DID]: hcoelho, jardim, gamito, faugusto, peixinho, aragao, vamaro => Linux group B [DGM]: hcoelho, jardim, smatias => Linux group C [SAD]: hcoelho, jardi, yesenia => Linux group D [NTL]: Everybody =>...

Hello. Id like you to help me with something here, i m trying to create a log manager just for fun but I cant make it work since the beginning. I created my project, configured the db, initialized my model Log and then my controller Log, then im just trying to use scaffold :log for the controller but I cant make it work, the error says: undefined method `errors'' for

Hi all, How can I overcome the error "object of type 'closure' not subsettable" I ran the following script seq <- paste(seq(1914, 1916, by=1), "*.y", sep=".") # make sequence c <- 3 # total number of files d2 <- file # creates dummy file # Input sequence in loop for (i in 1:3){ list <- list.files("~/ukcp09/txt/x.djf", seq[[i]]) file

Hello, Does anybody know if the outscale option of randomForest yields the standarized version of the outlier measure for each case? or the results are only the raw values. Also I have notice that this measure presents very high variability. I mean if I repeat the experiment I am getting very different values for this measure and it is hard to flag the outliers. This does not happen with two other

Dear R-users Coming from a proc mixed (SAS) background I am trying to get into the use of (n)lme. In this connection, I have some (presumably stupid) questions which I am sure someone out there can answer: 1) With proc mixed it is easy to get a hold on the estimated variance parameters as they can be put out into a SAS data set. How do I do the same with lme-objects? For example, I can see the

Hello. I'm having trouble installing rgl. I have a theory as to the problem. First, the error message and session info. > install.packages("rgl") trying URL 'http://probability.ca/cran/src/contrib/rgl_0.81.tar.gz' Content type 'application/x-gzip' length 1636939 bytes (1.6 Mb) opened URL ================================================== downloaded 1.6 Mb *