similar to: NADA

Hello Rui, Thanks for your helpful suggestions. Just for illustration, let's use the well known Duncan dataset of prestige vs education + income that is contained in the "car" package. Suppose I wish to use boot function to bootstrap a linear regression of prestige ~ education + income and use the following script: duncan.function <- function(data, indices) {data =

Dear All, I have two simple data samples (no groups or factors, etc.) and would just like to compute the two-sample Wilcoxon Rank Sum test using the wilcox_test function contained in the coin package, which is reportedly better than the regular wilcox.test function because it performs some adjustment for ties. Would anyone know how to craft a script to perform this task? Much appreciated. Janh

Greetings! We are trying to obtain confidence and prediction intervals for a predicted Y value from bootstrapped linear regression using the boot function. Does anyone know how to code it? Greatly appreciated. Janh [[alternative HTML version deleted]]

Hello, Much clearer now, thanks. It's a matter of changing the function boot calls to return the predicted values at the point of interess, education = 50, income = 75. I have changed the way the function uses the indices a bit, the result is the same, it's just the way I usually do it. pred.duncan.function <- function(data, indices) { mod <- lm(prestige ~ education +

R-help is not a free coding service. We expect users to make the effort to learn R and *may* provide help when they get stuck. Pay a local R programmer if you do not wish to make such an effort. Cheers, Bert On Oct 14, 2017 7:58 AM, "Janh Anni" <annijanh at gmail.com> wrote: Greetings! We are trying to obtain confidence and prediction intervals for a predicted Y value from

Hello All, Would anyone be able to help me understand how R computes a quantile-quantile plot for comparing two data samples with unequal sample sizes? Normally, the procedure should be to rearrange the larger data sample into n equally-spaced parts using interpolation, where n is the sample size of the smaller sample, and then plot the matching data pairs. I tried using different plotting

The sample data sets that come with the NADA package are limited to one or two variables and a censored measurement indicator column. I try to mimic examples using my data but keep missing the target. My water chemistry data is available in two formats: long (as seen in a database table) and wide (as seen in a spreadsheet). The two structures are: str(chem) 'data.frame': 65349 obs. of

Hi, I try to boxplot following data on the subset of (V1,V3,V5,V7) and (V2,V4,V6,V8) V1 V2 V3 V4 V5 V6 V7 V8 2 4 6 7 12 33 43 53 how can I use boxplot function to plot it? thanks, William

I need a nudge in the right direction to get started using NADA. I bought Helsel's second addition and am currently reading it; NADA is installed in R. My data has been restructured with a couple of awk scripts. The data frame structure now has a flag if the quantity is censored (ceneq1 column) as well as a lower and upper limit for censored data. For present purposes, interval censoring

hi there, trying to follow the steps of Yan & Rossini 2003, I have two problems : first when I 'make CrossCompileBuild', I get : ****************************************************************************** make[1]: Entering directory `/home/ros/RWORK/CROSSSOMP/WinR/R-1.8.0/src/gnuwin32' make -f Makefile.docfiles make[2]: Entering directory

I need to make nonlinear regression with the posterior script, but how is the problem? I have error in library (nls), package 'nls' has been merged into 'stats'. I need help? What other forms I have to make nonlinear regression? and how I find to calculate statistics y residuals, scatterplot. thanks SCRIPT ros<-read.table("Dataset.csv",header=T,sep=",")

Hello all, I'm quite a newbie to wine, and have some basic things running under it like WT Library, & Seamonkey. I need to install the local Revenue software because they sent me all tax documentation by email now :-(. My PC: Athlon 2.6Ghz with 1 gig, 40G hd, running Fedora Core 5 ROS ((Irish) Revenue Online Service) System Requirements

I have water chemistry data with censored values (i.e., those less than reporting levels) in a data frame with a narrow (i.e., database table) format. The structure is: $ site : Factor w/ 64 levels "D-1","D-2","D-3",..: 1 1 1 1 1 1 1 1 ... $ sampdate: Date, format: "2007-12-12" "2007-12-12" ... $ preeq0 : logi TRUE TRUE TRUE TRUE TRUE

I want to get access to the code of an S4 method in order to possibly modify a function to accomplish my particular needs: in my case the function in is cenfit() from the package NADA Now, given my reproducible example: my.ex<-structure(list(TEC = c(0.21, 0.077, 0.06, 0.033, 0.014, 0.0072), LR = c(0L, 0L, 1L, 0L, 1L, 0L)), .Names = c("TEC", "LR"), class =

Check this first. https://www.itechtics.com/3-ways-showhide-security-tab-windows-10/ Greetz, Louis > -----Oorspronkelijk bericht----- > Van: samba [mailto:samba-bounces at lists.samba.org] Namens > Rowland Penny via samba > Verzonden: vrijdag 23 februari 2018 11:03 > Aan: samba at lists.samba.org > Onderwerp: Re: [Samba] R: Missing 'security' tab > > On

Rich, The documentation for cenboxplot states that the second argument must be logical and not integer. the function cenboxplot substitutes synthetic values for censored values using ros, hence the error message from the ros method. I also do not understand how you expect group = 'SO4' to work. It is not clear that the function will replicate 'SO4' to make a single group.

Hello, I am trying to conduct a Hausman test to choose between FE estimators and RE estimators. When I try to run: library(plm) fixed <- plm(ROS ~ DiffClosenessC +ZZiele + AggSK + nRedundantStrecken + Degree + KantenGew + BetweennessC + SitzKappazitaet, data=Panel,index=c("id","time"),model="within") summary(fixed) fixef(fixed) random <-plm(ROS ~

The cenmle function is used to fit two sets of censored data and test if they are significantly different. I can print out the results of the analysis on the screen but can't seem to figure out how to access these results in R and assign them to new variables, e.g., assign the slope calculated with cenmle to the variable m. Any suggestions? Tom -- View this message in context:

One set of data has censored (less-than detection limits) water chemistry concentrations for 80-100% of all observations. My initial trial-and-error attempts to apply the cenboxplot() method suggests that it has an upper limit to the percentage of censored observations. I do not see this limit in Dennis Helsel's second edition. Has anyone experience plotting censored data and can provide

Hi, I am using the cenmle function from the NADA package and some of my data are negative values. As a result the cenmle function will not work and NaN's are produced. I try to change the distribution to Gaussian, but it still will not run. Could somebody please help me with this? Thanks -- Shane [[alternative HTML version deleted]]