similar to: lmomco package and confidence limits?

Hi, I am trying to automate the creation of tables for some simply analyses. There are lots and lots of tables, thus the creation of a user-defined function to make and output them to excel. My problem is that some of the analyses have convergence issues, which I want captured and included in the output so the folks looking at them know how to view those estimates. I am successfully able to do

1. I did not attempt to sort through your voluminous code. But I suspect you are trying to reinvent wheels. 2. I don't understand this: "I've failed to find a solution after much searching of various R related forums." A web search on "error handling in R" **immediately** brought up ?tryCatch, which I think is what you want. If not, you should probably explain why it

You can capture warnings by using withCallingHandlers. Here is an example, its help file has more information. dataList <- list( A = data.frame(y=c(TRUE,TRUE,TRUE,FALSE,FALSE), x=1:5), B = data.frame(y=c(TRUE,TRUE,FALSE,TRUE,FALSE), x=1:5), C = data.frame(y=c(FALSE,FALSE,TRUE,TRUE,TRUE), x=1:5)) withWarnings <- function(expr) { .warnings <- NULL # warning handler will

tryCatch() is good for catching errors but not so good for warnings, as it does not let you resume evaluating the expression that emitted the warning. withCallingHandlers(), with its companion invokeRestart(), lets you collect the warnings while letting the evaluation run to completion. Bill Dunlap TIBCO Software wdunlap tibco.com On Tue, Mar 6, 2018 at 2:45 PM, Bert Gunter <bgunter.4567 at

Hi guys, I'm trying to use lmomco package. first I did the manual calculation on what is the estimates scale and location parameter given L-CV=0.2, L1=1000 L-moments and k (shape parameter) =- 0.1. so what i get is: location: 821.0445 scale: 260.7590 shape: -0.1000 #I assign this as GEV vectors using vec2par GEVpara2<-vec2par(c( 821.0445 , 260.7590 ,-0.1),'gev') #then I

Dear List I?m using the ?lmomRFA? package to fit different distributions to my data sample. To calculate the error bounds I used: regsimq(?) and sitequantbounds(?) So my questions are: Are error bounds and confidence intervals the same thing? And: Does regsimq(? boundprob = c(0.05, 0.95)) calculate the 90 or the 95% confidence interval? If error bounds and confidence intervals are

Dear R forum >From the given data, I have estimated the parameters of Wakeby distribution using lmomco package as library(lmomco) (amounts <- read.csv("input_S.csv")$amount) # ___________________________________________________________ # Wakeby distribution - Parameter estimation N                      = length(amounts) lmr                    = lmom.ub(amounts)

Does anyone know how to do a fry plot in R? I have 600-800 points per image, and I really don't want to attempt that manually. Thank you! Cassie [[alternative HTML version deleted]]

I am trying to age standardize using the svystandardize package in R. I have successfully managed to hit my SUDAAN based targets for estimates by sex, but not the total. The total is only a little different, but I'd like some help knowing why it isn't exact. I've included the SUDAAN code that generates the targets and my R script (and output) that I have so far. I can't supply the

Hi, I am using parnor function of lmomco package. I believe it provides mean and std. dev for the set of data. But the std. dev provided does not match with the actual std. dev of the data which is 247.9193 (using sd function). Am I missing something here? > lmr <- lmom.ub(c(123,34,4,654,37,78)) > parnor(lmr) $type [1] "nor" $para [1] 155.0000 210.2130 >

Dear R forum, I am bit confused and please guide me - (1) Is "Pearson Type III Distribution" as given in lmomco package same as Three Parameter Pearson 5 Distribution? If not, how do I estimate the parameters of Three Parameter Pearson 5 Distribution? (2) Is there any other R forum dealing with only Statistical queries? Kindly guide Regards Katherine [[alternative HTML version

Hi R-users, I'm trying to repeat the same procedure to 1000 data set. I know this is very easy, but I got stuck finding the right and fastest way in running it. IID50=Riidf[1:50,1:1000] #where IID50 is a dataframe consist of 1000 time series(as column) and 50 time scales (row). #what I tried to do: estIID50=rep(NA,1000) for (i in 1:1000) estIID50[i]=pargev(lmom.ub(IID50[1:50,i])) #warning

Hello All, When I use bw.SJ (based on Sheather & Jones, 1991) in R to estimate the bandwidth for a highly skewed data, I get the following message: "sample is too sparse to find TD". I played around with the parameters such as no. of bins (nb), lower, upper (range over which to minimize) to no avail. My sample size is 50,000. Can anyone tell me what this means and of some

respectd sir,  i am working on regional frequency analysis of flood flow data and want to use packages lmomRFA, lmom, imomco,etc andby using these i want to find RMSE , errorbounds and absolute biase but could not make it clear to use, please help regards atta muhammad asif assistant professor atta_ycc@yahoo.com [[alternative HTML version deleted]]

Can anyone show me how to add a log Pearson type III plot using the evdistq() command to an extreme value plot using the lmom package? Attached sample code below... Thanks in advance, Dave library(lmom) # annual maximum daily streamflows Mackenzie River mackenzieRiver = c(26600, 30300, 34000, 32000, 29200, 28300, 28600, 26400, 28300, 28800, 29000, 22100, 32900, 31800, 21600, 32100, 27000,

Hi List! I?m using regsimq() from the ?lmomRFA?-package to calculate error bounds for diverse distributions. For example: regsimq(gumfit$qfunc, nrec = lmom.data$n, f = lcdfgum, boundprob = c(0.025, 0.975)) Several times I got this error massage: Fehler in quantile.default(ou, probs = boundprob, type = 6) : missing values and NaN's not allowed if 'na.rm' is FALSE So my question

Dear R users, I am new in R and I want to use the nls package to analyze some experimental data. The data is in the attached file "data". It is the response "Sav" measured at different "C0". Basically, the "C0" is a function of C1, K2, and r, and the "Sav" is a function of C0, C1, K2, and r. The math equations are shown in the attached

Dear colleagues, I'm conducting a meta-analysis of studies evaluating adherence of HIV-positive drug users into AIDS treatment, therefore I'm looking for some advice and syntax suggestion for running the meta-regression using proportions, not the usual OR/RR frequently used on RCT studies. Have already searched already several handbooks, R-manuals, mailing lists, professors, but... not

Hola a todos, Estoy utilizando la libreria e1071 para clasificar unos datos. Me gustaría poder guardar los modelos en el disco duro y no en memoria de R. He visto que hay una función: save : que guarda el modelo en memoria y load: que carga ese modelo Ejemplo: #saving the best model save(calibrate.rf.model1, file=''bestmodel.rda'') #loading the best model

Hello everyone I hope Im not bothering you all again. I have just begun to use R and so Im not yet familiarized with it.. I ve got an assignment which consists in calculating the area of a circle given a certain radius and center using the monte carlo method, which means that I have to plot a circle given its parameters. Limit the area inside it...with as many sample points as possible...and