similar to: Winsorized mean and variance

Hello R list, I have a little problem with splom. I'd like to wrap it in a function, for example: multi.scatterplot <- function(data,groups,cols,colors) { splom(~data[,cols], groups = as.symbol(groups), data = data, panel = panel.superpose, col=colors) } and then call it like in multi.scatterplot(iris,"Species",1:4,c("green","blue","red"))

Hi All, I want to take a matrix (or data frame) and winsorize each variable. So I can, for example, correlate the winsorized variables. The code below will winsorize a single vector, but when applied to several vectors, each ends up sorted independently in ascending order so that a given observation is no longer on the same row for each vector. So I need to winsorize the variable but

Hi everyone! #I need a loop or a function that creates a X2 variable that is X1 without the extreme values (or X1 winsorized) by industry and year. #My reproducible example: firm<-sort(rep(1:1000,10),decreasing=F) year<-rep(1998:2007,1000) industry<-rep(c(rep(1,10),rep(2,10),rep(3,10),rep(4,10),rep(5,10),rep(6,10),rep(7,10),rep(8,10),rep(9,10), rep(10,10)),1000) X1<-rnorm(10000)

Dear R-Experts I am sure this might look simple question for experts, at least is problem for me. I have a large data frame with over 1000 variables and each have different distribution( i.e. have different quantile). I want to create a new grouped data frame, where the new variables where the value falling in first (<25%), second (25% to <50%), third (50% to <75%) and fourth quantiles

Can anyone please tell me if there is a function to calculate confidence intervals for the results of the quantile function. Some of my data is normally distributed but some is also a squewed distribution or a capped normal distribution. Some of the data sets contain about 700 values whereas others are smaller with about 100-150 values, so I would like to see how the confidence intervals change

Dear Group I have few errors while installing package rattle from CRAN i do the installing from the local zip files... I am using R 2.4.0 do i have to upgrade to R2.4.1 ? ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ utils:::menuInstallLocal() package 'rattle' successfully unpacked and MD5 sums checked updating HTML package descriptions > help(rattle) No

Hello, I have a question that is not directly related to R ... but I try to do it in R ;-) : I would like to generate a matrix Q satisfying (for a given Z, X and W) the two following conditions: t(Q)%*%Q=Z (1) XQ=W (2) where: Q is m rows and r columns X is p rows and m columns D is p rows and r columns C is r rows and r columns with m>p,r e.g: m=6, p=2 r=3

Hello, I'd need to compute the calibration error of posterior class probabilities p(y|x) estimated by using rpart as classification tree. Namely, I train rpart on a dataset D and then use predict(... type="prob") to estimate p(y|x). I've found the possibility to do that in the ROCR package, but I cannot find a link to a paper/book which explains the details of the

Hi everybody, I noticed a strange behavior when using loops versus apply() on a data frame. The example below "explicitly" computes a distance matrix given a dataset. When the dataset is a matrix, everything works fine. But when the dataset is a data.frame, the dist.for function written using nested loops will take a lot longer than the dist.apply ######## USING FOR ####### dist.for

Hi folks, I'm having some trouble understanding the intricacies of panel functions. I wish to create three side-by-side graphs, each with different data-- so far, so good: I rbind() the data, add a column of subscripts as a conditioning variable, load up the lattice package, specify either a c(3,1) 'layout' or work through 'allow.multiple' and 'outer' and I'm

hi I'm puzzled as to the relation between the REML score computed by gam and the formula (4) on p.4 here: http://opus.bath.ac.uk/22707/1/Wood_JRSSB_2011_73_1_3.pdf I'm ok with this for poisson, or for quasipoisson when phi=1. However, when phi differs from 1, I'm stuck. #simulate some data library(mgcv) set.seed(1) x1<-runif(500) x2<-rnorm(500)

Dear All, Currently i am using CentOS 4.4 Linux After i am rebooted the PC and I have getting problem like /etc/rc.d/rc.sysinit: line 126: 442 Segmentation fault /bin/dmesg -n $LOGLEVEL Initialization hardware... /etc/rc.d/rc.sysinit: line 167: 450 Segmentation fault LC_ALL=C fgrep -xq "$1" /etc/hotplug/blacklist >/dev/null /etc/rc.d/rc.sysinit: line 167: 455 Segmentation

Hi, I was wondering whether anyone has managed to compile the LLVM libraries to run on iPhone? After compiling and installing the iPhone toolchain on MacBook running Leopard (10.5.2). I run configure In the llvm2.3 directory: ./configure –host=arm-apple-darwin –target=arm-apple-darwin –enable-optimized –enable-targets=arm I run make, including an override for TBLGEN (because I obviously

Dear R members, I would like to know if there is a stepwise function that gives LRT value instaed of AIC? thanks in advance, Luciana Pozzi CENPAT-CONICET Argentina [[alternative HTML version deleted]]

Hello, I am very new to R. I would like to plot astronomy data by right ascension, declination, and various "distance" values, such as redshift and comoving distance, in 3D. Is there any 3D polar plotting functions? I can't seam to locate any information on whether it exists or not.

llvm-project/llvm> git rev-parse --git-common-dir --git-common-dir llvm-project/llvm> On Fri, May 3, 2019 at 3:48 PM Craig Topper <craig.topper at gmail.com> wrote: > I might have seen this issue. What does "git rev-parse --git-common-dir" > return for you? > > ~Craig > > > On Fri, May 3, 2019 at 12:35 PM Cameron McInally via llvm-dev < > llvm-dev

useRs ? I?m working with the attached data that contains one year?s worth of sub-daily observations of flow (?Q?) and specific conductance (?SC?, a surrogate for concentration) at a point in a stream. The R code posted below shows the extent of data processing thus far. My goal is to create a covariance matrix that takes on the following form: Q1 Q2 ? Q365 SC1 SC2 ? SC365 Q1 Q2 ? Q365

Estimados tengo 1 lista de valores aleatorios misma que necesito reordenar de acuerdo a determinadas condiciones. Ej: sea la siguiente lista (la primera fila hace referencia a la posición del valor de lista 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 1 7 3 4 8 6 1 4 5 2 3 4 10 15 5 necesito que me tome los valores que son mayores 3 y menores a 8, pero unicamente considerando los valores hasta la

Hello, I'm new using R. I'm trying to develop a K-means Clustering with R for some data I have, however each time I use that instruction with the same data my cluster means, clustering vector and within cluster sum of square change and I don't understand why because I use the same parameters and the same data. Can anybody explain me why does it happen? Thank you Act. Calef

Dear all, I had a R segmentation fault, and then invoked debug mode and ran step by step. When I reached "terms(Y~X1*X2*...*X16)", I would then have "segmentation" fault. However, if I just ran this under regular "R interactive" mode, it would be fine though taking long time. My questions are: 1. Is there a known limit of terms for a formula? 2. Why does the