similar to: use abline() for regression model in the plot

Hi, i try to plot my graph into different device using x11(), but when i do this comes up: > x11(print(plot(A5e$ECAB,A5e$EXPEND,type='p',main='Per capita expenditure > against economic + ability index without outliners',xlab='ECAB',ylab='EXPEND',xlim=c(0,150),ylim=c(150,400), + col='red',col.axis='blue'))) NULL Warning message:

Hi, I was wondering if there is a ready made function or parameter for indicating the sample size in boxplots? Here's what I came up with so far: library(ISwR) data(energy) attach(energy) boxplot(expend~stature) sample.size <- tapply(expend, stature, length) sample.size <- paste("N=", sample.size, sep="") mtext(sample.size, at=1:length(unique(stature)), line=2,

> corruption happens only in this cases: > > - volume with shard enabled > AND > - rebalance operation > I believe so > So, what If I have to replace a failed brick/disks ? Will this trigger > a rebalance and then corruption? > > rebalance, is only needed when you have to expend a volume, ie by > adding more bricks ? That's correct, replacing a brick

I am using barplot, and barplot2 in the gregmisc bundle, in the following way: barplot2(sort(xtabs(expend / 1000 ~ theme)), col = c(mdg7, mdg8, mdg3, mdg1), horiz = T, las = 1, xlab = "$ '000", plot.grid = T) The problem is that the values of 'theme', which is a factor, are in some cases rather long, so that I would like to wrap/split them at a space once they

Hello, I am relatively new to R and I am working through the code that is provided in the book Numerical Ecology with R and I have run across an error message that I can't seem to figure out. I am using the vegan, ade4, gclus and cluster packages. The code is as follows: # Ordered community table # Species are ordered by their weighted averages on site scores or <- vegemite(spe,

Hi all, >From a list of strings, I desire to filter out the followings: 1. Digits at the beginning of the strings 2. Character "SPE" following the digits (if it exists) 3. Any characters followed by hyphen The following produces the desired result, but would like to know whether this can be done more efficiently. Any suggestions would be much appreciated. dat <- c("2148

Just to clarify as i'm planning to put gluster in production (after fixing some issue, but for this I need community help): corruption happens only in this cases: - volume with shard enabled AND - rebalance operation In any other cases, corruption should not happen (or at least is not known to happen) So, what If I have to replace a failed brick/disks ? Will this trigger a rebalance and

Dear all, I have run the confint() function according to below and I get the following error message: > confint(stepAIC.glm.spe.var.konn2.abund, level=0.95) Waiting for profiling to be done... Error: missing value where logical needed In addition: Warning message: NaNs produced in: sqrt((fm$deviance - OriginalDeviance)/DispersionParameter) or > confint(stepAIC.glm.spe.var.konn2.abund,

Instead of > m <- c(400000000, 50000000000) > paste("A", m, "B", sep="") [1] "A4e+08B" "A5e+10B" I want "A400000000" and "A50000000000"

On Mon, Jun 3, 2019 at 6:54 PM Marc Schwartz wrote: > I am on macOS primarily, albeit, I have run both Windows and Linux routinely > in years past. With all due respect, then you have no business in this thread. > That being said, these days, I do run Windows 10 under a Parallels VM on > macOS, as I have a single commercial application that I need to run for > clients now and

I think that the time has come for me to back off a bit from my involvement with Shorewall. I just don''t have enough cycles (or energy) to keep up the pace of the last several months. As a consequence, I''m going to do the following: 1. I''m going to stop personally supporting the entry level tools (samples and quick start guide). These tools are a source of constant

Hi, Sorry for the confusion. I would like to estimate a model wherein the marginals of z with respect to w1 and w2 are smooth functions of x and y. I have data on z, x, y, w1 and w2. so E[dz/dw1] = f(x,y) and E[dz/dw2] = g(x,y) and I would like to estimate f(x,y) and g(x,y) I suppose I could try to fit something more general using projection pursuit, but the nature of the problem suggests

Is anybody getting acceptable performance with a Rails application on Dreamhost? By this, I mean response times no greater than 3 seconds, and no large numbers of FCGI processes that you have to kill manually. If so, how? I''ve ended up running a constant ping script (once every ten minutes), and still have to kill some number of excess dispatch.fcgi processes every day. Along with

Hi Chandler, I am now working on a bisecting tool to find mis-optimization on LLVM. I found OptBisect a very useful option and hope to make it work on the new pass manager. I have several questions about it. 1. Any plans to apply codegen stage with new pass manager? IIUC, new pass manager only works for opt stage. However the OptBisect option tries to accumulate pass counts through opt stage

Dear R-gurus, Just supose I have a dara.frame that looks like myDF<-read.table(stdin(),head=T,sep=",") codID,namesp,k1,k2,k3,k4 1,spA,2,5,6,3 2,spB,4,5,4,6 3,spC,2,1,5,6 4,spC,5,4,3,2 5,spD,1,2,3,4 6,spE,2,4,3,1 I need to update the columns k1-k4 with the namesp, but considering the math between Kx and codID. My desired output must looks like: codID,namesp,k1,k2,k3,k4

Dear list, I want to calculate SPE (squared prediction error) in x-space, can someone help? Here are my codes: fit.pls<- plsr(Y~X,data=DAT,ncomp=3,scale=T,method='oscorespls',validation="CV",x= T) actual<-fit.pls$model$X pred<-fit.pls$scores %*% t(fit.pls$loadings) SPE.x<-rowSums((actual-pred)^2) Am I missing something here? Thanks in advance. Stella Sim

Dear all, Is there any function for calculating confidence limits for coefficients in an lm() object? I know of the confint() function in the MASS library working very well on my binomial GLMs and I have tried it (using glm () , family=gaussian) but it gives NAs according to below. Does the confint() function not accept gaussian GLMs? Could there be convergence problems in the GLM? Note the

Hello, I am using vegan to do an NMDS plot and I would like to suppress the labels for the loading vectors. Is this possible? Alternatively, how can I avoid overlap? Many thanks for the help. Example code: #perform NMDS using metaMDS() function spe.nmds<-metaMDS(data, distance='bray',k=2 , engine = "isoMDS", autotransform=F, trymax=1000) #calculate the loading (i.e.,

Hello, I've got an ISDN PRI circuit terminating in a Cisco AS5350, which in turn is talking to an Asterisk server via SIP for call origination and termination. Seems simple enough, and it works for the most part, but: 1) Caller ID name data comes in on the PRI, but doesn't appear to get handed off to the Asterisk server via SIP, at least not in any format that Asterisk

Hi, I am a new R user and am trying to construct a palaeoenvironmental transfer function (weighted averaging method) using the package rioja. I've managed to insert the two matrices (the species abundance and the environmental data) and have assigned them to the y and x values respectively. When I try and enter the 'WA' function though, I get an 'Error in FUN' message (see