search for: eff

Hello list, I've seen the recent discussions documenting problems with lm. I have encountered the following problem. I use WinXP Pro with service pack 1, and R 1.9.0, on a XEON 2GHz, with 1GB of RAM. > eff.fro std.dev mean NSTRDSP 7.403749e-01 1.215686e-01 CPFGEP 9.056763e+00 1.815686e+00 WSWOLF 4.703588e+05 1.112832e+05 NPILGRIM 1.017640e+06 2.134335e+05 WSNMILE 1.367312e+07 1.892021e+06 WSHIDESL 1.830811e+07 1.892021e+06 > reg <- lm(log(mean) ~ log(std.dev), data=e...

Hello! I am trying to put two math expressions in the title of a plot. As you can see below, I can place correctly one expression at a time, but not both. Ideally I would like to have them separated by a comma. Any suggestions? > k <- 1 > n.eff <- c(20, 30) > ### this works > plot(0,0, main = substitute(n == k, list(k = k))) > ### this works > plot(0,0, main = substitute(N[eff] == neff, list(neff = n.eff[k]))) > ### this doesn't work > plot(0,0, main = substitute(n == k * N[eff] == neff, list(k = k, neff = n.eff[k...

...for(i in 1:n) { mat <- c(mat,((1-i):(n-i))) } mat <- abs(mat) mat <- corr(mat,d) mat <- matrix(mat, ncol=n) mat } corr <- function (x,d) { rho <- gamma(x+d)*gamma(1-d) rho <- rho/(gamma(x-d+1)*gamma(d)) rho } tr <- function(M) { sum(diag(M)) } eff <- function(d, n=100, tol=1e-7) { x <- 1:n X <- cbind(1,x) V <- corr.matrix(d, n=n) V1 <- qr.solve(V, tol=tol) GLS <- (t(X)%*%V1%*%X) GLS <- solve(GLS) %*%t(X)%*%X OLS <- solve(t(X)%*%X)%*%t(X)%*%V%*%X eff <- tr(GLS)/tr(OLS) eff } eff.all <- function...

...ximum Latency' is _too much high_ tiobench ======== Sequential Reads ext3 File Blk Num Avg Maximum Lat% Lat% CPU Kernel Size Size Thr Rate (CPU%) Latency Latency &2s &10s Eff ---------------------------- ------ ----- --- ------------------------------------------------------------ 2.4.20-pre10 8192 4096 1 37.28 21.64% 0.312 135.50 0.00000 0.00000 172 2.4.20-pre8 8192 4096 1 37.34 21.38% 0.311 133.76...

As far as I know, btrfs comes from "btree file system", but why does btrfs pronounce "butter-eff-ess"? -- Wang Shaoyan -- To unsubscribe from this list: send the line "unsubscribe linux-btrfs" in the body of a message to majordomo@vger.kernel.org More majordomo info at http://vger.kernel.org/majordomo-info.html

...The only thing this method does, is change xy into a matrix and then pass on to the next method using callGeneric again. Arguments are passed. Now this all works smoothly, as long as you stay in the global environment : require(raster) a <- raster() a[] <- 1:ncell(a) origin <- c(-80,50) eff.dist <- 100000 unlist(xyValues(a,xy=origin,buffer=eff.dist)) [1] 14140 14141 14500 14501 Now let's make a very basic test function : test <- function(x,orig.point){ eff.distance <- 100000 p <- unlist(xyValues(x,xy=orig.point,buffer=eff.distance)) return(p) } This giv...

...R user, finding my own way through; learning code etc. at the same time as learning the statistics I need to complete my project. What I have is the code from Stata and am trying to reproduce the same analysis in R - my program of choice. . glm count md ms rf sg, family(poisson) exposure(effort) eform I am lost at the point of finding the equivalent code for 'exposure'. Having looked at a few forums and 'googled'. I thought 'offset', used as offset=(log(Eff)) or the equivalent +offset(log(Eff)) would produce the desired effect. Incidentally my co...

In the course of exploring response prediction, I stumbled upon a small discrepancy between the CIs produced by predict.lm() and all.effects() require(mlmRev) require(effects) hsb.lm <- lm(mAch ~ minrty * sector, Hsb82) hsb.new <- data.frame( minrty = rep(c('No', 'Yes'), 2), sector = rep(c('Public', 'Catholic'), each = 2)) hsb.eff <- all.effects(hsb.lm) cbind( hsb.new, pr...

Hi All, I am trying to use optim() to minimize a function with a penalty function term. This is a simple model bioeconomic model of a fishery. The penalty function constrains the amount of effort (f) at 9. This works fine. The code is: ********** nfleets<-2 M<-1 M<-array(M,dim=c(nfleets)) N<-1000 cost<-c(30,30) cost<-array(cost,dim=c(nfleets)) Price<-2 Price<-array(Price,dim=c(nfleets)) q<-array(0.1,dim=c(nfleets)) f<-1 f<-array(f,dim=c(nfleets)) f1<-f...

I have tried to set parameters of an effects plot with moderate success, but I'm finding that although the plot(effect()) method uses trellis graphics (via the lattice package), not all aspects of the plot can be controlled by standard trellis graphics. I am able to set any parameters that are given in trellis.par.get() easily enough,...

Ok. It's a Firefox Add-on: https://www.eff.org/https-everywhere Questions: 1) But: Why can't i find it on the offical Firefox Add-ons site?: https://addons.mozilla.org/en-US/firefox/ 2) Did anyone audited the "HTTPS Everywhere" code? 3) Can someone trust this Add-on? Is it safe to install/use? 4) If it's so great why...

...0.4332, 0.652, 0.6524, 0.6527, 0.653, 0.6534, 0.6537, 0.6541, 0.6544, 0.6547, 0.6551, 0.6554) hc <- c(rep(1.05e-02, 2), rep(2.09e-02, 3), rep(4.19e-02, 15)) hd1B <- c(0.1109, 0.1381, 0.1485, 0.1637, 0.2446, 0.2497, 0) test.example <- PwrGSD( EfficacyBoundary=LanDemets(alpha=0.05, spending= ObrienFleming), FutilityBoundary=LanDemets(alpha=0.1,spending=ObrienFleming), RR.Futility = 0.82, sided="<",method="A",accru =7.73, accrat=9818.65, tlook =tlook, tcut0 =t0, h0=h0, tcut1=t0, rhaz=rhaz,...

...0.4332, 0.652, 0.6524, 0.6527, 0.653, 0.6534, 0.6537, 0.6541, 0.6544, 0.6547, 0.6551, 0.6554) hc <- c(rep(1.05e-02, 2), rep(2.09e-02, 3), rep(4.19e-02, 15)) hd1B <- c(0.1109, 0.1381, 0.1485, 0.1637, 0.2446, 0.2497, 0) test.example <- PwrGSD( EfficacyBoundary=LanDemets(alpha=0.05, spending= ObrienFleming), FutilityBoundary=LanDemets(alpha=0.1,spending=ObrienFleming), RR.Futility = 0.82, sided="<",method="A",accru =7.73, accrat=9818.65, tlook =tlook, tcut0 =t0, h0=h0, tcut1=t0, rhaz=rhaz,...

Hi, I have some values in a list format generated by the following: Path_Number <- 0010 ID.Path <- formatC(0001:Path_Number, width=4, flag=0) # Make vector of ID's. No_of_Effectors <- sample(1:550, length(ID.Path), replace=TRUE) # Define Number of Effectors each individual gets. Effectors <- split(sample(1:10000, sum(No_of_Effectors), replace=TRUE), rep(ID.Path, No_of_Effectors)) # Generate effectors and dish them out. Effectors And I've written a chunk which...

...essful Centos 'brand name' :-) The user-agent string is one of the items used in uniquely identifying/fingerprinting a user/machine, so the more generic it is the better. Including the details of the OS add to the "bits of identifying information" available to trackers. See the EFF testing site for more details: <https://panopticlick.eff.org/>

...sterisk*CLI> core show channel CBAnn/207-0000067f;1 -- General -- Name: CBAnn/207-0000067f;1 Type: CBAnn UniqueID: 1398809161.20186 LinkedID: 1398809161.20186 Caller ID: (N/A) Caller ID Name: (N/A) Connected Line ID: (N/A) Connected Line ID Name: (N/A) Eff. Connected Line ID: (N/A) Eff. Connected Line ID Name: (N/A) DNID Digits: (N/A) Language: en State: Up (6) NativeFormats: (nothing) WriteFormat: unknown ReadFormat: unknown WriteTranscode: No ReadTranscode: No Time to Hangup: 0 Elapsed Time: 0h1m3s Bri...

Dear everybody! Need help with graphics. I am runnig a simple lm and then using allEffects from 'effects' package: require(effects) model<-lm(Y~A+B, data=mydataframe) I am trying to build (for each predictor - A and then B) a plot of means on Y. I was successful doing it like this - in one swoop: ml.eff<-allEffects(ml1, se=F) plot(ml.eff,ylab="Title of Y")...

Greetings. I'm trying to replicate an analysis I did a few years ago, then in SPSS, using the SPSS GLM command: GLM n_diffpt WITH age_i inc_i join_i work_i educ_i give_i cs_i eff_i age_a inc_a join_a work_a educ_a give_a cs_a eff_a /METHOD = SSTYPE(3) /INTERCEPT = INCLUDE /PRINT = PARAMETER ETASQ RSSCP GEF /CRITERIA = ALPHA(.05) /DESIGN = age_i*age_a inc_i*inc_a join_i*join_a work_i*work_a educ_i*educ_a give_i*give_a cs_i*cs_a eff_i*eff_a...

...(0,4), rep(.96,3), rep(1.92,2), 2.88) z<-c(.65,1.78,2.4,2.54,3.04,2.22,2.97,3.56, .745,1.287,1.98,1.91,2.02,2.24,2.95, .47,.80,.90,1.32,1.57,2.51, .23,.51,1.38,2.09,2.74, .14,.88,1.87,2.84, .40,1.61,2.48, .78,1.93, .87) These data sample one triangular corner of a matrix. I want a plot of the efficiency zideal<- abs(x-y)/.3 eff<- (z/zideal)^2 But I can't figure out how to get a persp() plot of eff; this doesn't work: persp(x,y,eff,col="grey",theta=45,phi=30) I tried a trick using unique(x),unique(y) which didn't work. If I define x, y and z like this: x<-c...

..., 1] smbd/server.c:main(628) smbd version 2.0.5a started. Copyright Andrew Tridgell 1992-1998 [1999/08/26 12:33:30, 1] smbd/files.c:file_init(216) file_init: Information only: requested 10000 open files, 246 are available. [1999/08/26 12:33:30, 0] smbd/uid.c:become_gid(105) Couldn't set effective gid to 65534 currently set to (real=0,eff=0) [1999/08/26 12:33:30, 0] smbd/uid.c:become_gid(107) Looks like your OS doesn't like high gid values - try using a different account [1999/08/26 12:33:30, 0] smbd/service.c:make_connection(463) Can't become connected user! I can not con...