similar to: Robust variance estimation with rq (failure of the bootstrap?)

I am interested in comparing the fit of robust (i.e., S and MM) and non-robust (i.e., OLS) estimators when applied to a particular data set. The paper entitled "A comparison of robust versions of the AIC based on M, S and MM-estimators" (available at: http://ideas.repec.org/p/ner/leuven/urnhdl123456789-274771.html) presents formulas for robust Akaike information criteria (AIC) for the M,

Hi all, I am attempting to create a new variable based on values of other variables. The variable is called pharm. It basically takes the numeric code of 1 as yes and 0 to be No from the variable B20_C1 (a question on a survey). However, I would also like to have a level for non-respondents and these are captured in the variable nr.B20C, which is a 1 when there is a non-response on the whole group

May I inquire about how to obtain the coefficients from a loess fit of multiple covariates? I have tried terms and coef, but it does not provide me with the coefficients to identify the terms of the entire regression equation. Many thanks Paul -- Paul R. Hutson, Pharm.D. Associate Professor UW School of Pharmacy 777 Highland Avenue Madison WI 53705-2222

Hello R users, I am a fairly new R user, however, one of the problems I am having is the use of applying if-, and-, then-, else- statements in R against datasets & dataframe. Rtips mentions it, but without simple examples I am not quite sure what my syntax should look like. I am particularly interested in syntax that can be used with the transform function. for example if foo[,1] ==

I'm using the lme function in nmle to estimate the variance components of a fully nested two-level model: Y_ijk = mu + a_i + b_j(i) + e_k(j(i)) lme computes estimates of the variances for a, b, and e, call them v_a, v_b, and v_e, and I can use the intervals function to get confidence intervals. My understanding is that these intervals are probably not that robust plus I need intervals on the

Hi all: After calculating in R I want to show the answer and some explanations not graphic plots in another new device and then save it as txt.file. However, I couldn't find any package or command to do it yet. I know some commands about generating graphics on different types of display or printing device. ex : windows(), postscript(), pdf() and etc. But these all works for graphic plots,

Hello, Please always cc the list. As for the question, yes, it does. If you want to remove just the ones with exactly 73.1 use the pattern grep("^73\\.1$", etc) Explanation: Beginning of string: ^ End of string: $ Escape special characters: \\ (needed because the period is a special character.) Hope this helps, Rui Barradas On 5/22/2018 12:50 PM, Ahmed Serag wrote: > Thank

One of my colleagues has a data set from a two-level nested design from which we would like to estimate variance components. But we'd like some idea of what the inevitable outliers are doing, so we were looking for something in R that uses robust (eg Huber) treatment and returns robust estimates of variance. Nothing in my collection of R robust estimation packages (robust, robustbase and MASS

Hello I am fitting a Cox PH model using the function coxph(). Does anyone know how to obtain the estimate of the covariance matrix under the null hypothesis. The function coxph.detail() does not seem to be useful for this purpose. Thanks, KD. [[alternative HTML version deleted]]

Dear R-experts, How can I remove a certain feature or observation by a part of its name. To be clear, I have a matrix with 766 observations as a rows. The row names are like this 70/556 71.1/280 72.1/556 72.1/343 73.1/390 73.1/556 Now I would like to remove all the rows that contain the text 73.1 Any ideas or suggestion please ? Regards ********************** Ahmed Serag Analytical

In theory ParkAndAnnounce has a lot of usefulness, however, that we've had very little success with application... Our application is similiar to the local Walgreens pharmacy.. Dr. Calls in, selects the "Im a doctor with a prescription option"... call is parked, and announcement overhead is given.. "Doctor holding on ${EXTEN}" exten =>

Dear All, I am experiencing a difficulty with the "lsoda" function in ODESOLVE library; my differentiation function contains a covariate changes with time, and I don't know how to incorporate that without using loops and so on. Does anyone know how to do it? Your help will be very much appreciated!! In-Sun In-Sun Nam School of Pharmacy and Pharmaceutical Sciences University of

On an Asterisk 1.2.12.1 system using a TDM400P with two FXO and two FXS modules, servicing two POTS lines: When dialing a number, such as a bank, or pharmacy, where it is required to enter a long series of numbers via the phone's keypad, an unexpected hangup occurs. The hangup does *not* occur when entering the numbers, as one might initially expect, but after the called end begins to

Dear R Mailing Subscribers, I just have a question of how R handles interacting variables in model creation using glm for instance. if I write : >> model=glm(solution~descriptor1+descriptor2+descriptor3,family=binomial(link="logit")) I should end up with coefficients for a logistic model that I can introduce easily in the mathematical form of such a model (weighted sum of

Hi, I want to be able to rename a folder using R, similar to file.rename. I want to paste " - done" onto the folder name. The reason for this is we run a loop on a large number of folders on a server and it would be nice for people to be able to log in and instantly see if their data has been processed so they can remove it. I have searched the help list but have not found anything and

Hello, I am using Samba on Linux as a CLIENT which is accessing Windows Server 2008 and I am trying to setup roaming profiles on the Linux/Samba client. My Linux distribution is openSUSE 11.3 RC1/Factory and YaST does a very good job at setting up Kerberos/Samba to join the domain. It all seems to work, except the Linux profile isn't synchronized back to the Windows domain server. I assume

I am new to R and am interested in using the program to fit quantile regression models to data collected from a multi-stage probability sample of the US population. The quantile regression package, rq, can accommodate person weights. However, it is not clear to me that boot.rq is appropriate for use with multi-stage samples (i.e., is capable of sampling primary sampling units instead of survey

Hi, I'm plotting PDFs and have a problem. If I have more than 8 sources of data the colours are repeated. These plots are used to remove poor data from the sets so it would be helpful if I could expand the colour range. Is there any way to do this? The plots are coloured by defining a vector colsVec<-1:(length(mzXMLfiles)) And defined in the loop for (file in 1:length(mzXMLfiles)) as

Dear R users & experts, I'd like to create a model using lm (or glm) under some constraints of how coefficients for each component could look like (sort of a range of coefficients that should be allowed). So let's go for an example : model=lm(age ~ eyecolor + height, data=inputdata) So let's suppose that R pops out a model with positive estimates for the coefficients eyecolor

Hi, I'm trying to make a legend with a justified list of numbers, so I thought I would use phantom(0) to align the 3-digit numbers properly with the 4-digit ones: legend(x, y, xjust=1, yjust=1, lty=c(1,2,3,4,5), adj=c(0,0.5), legend=expression(phantom(0)*300*plain(K), phantom(0)*550*plain(K), phantom(0)*830*plain(K),