similar to: boxplot by factor (Package base version 2.1.1) (PR#7976)

The issue is not with boxplot, but with split. boxplot.formula() calls boxplot(split(split(mf[[response]], mf[-response]), ...), but look at what split() returns when there are empty levels in the factor: > f <- factor(gl(3, 6), levels=1:5) > y <- rnorm(f) > split(y, f) $"1" [1] 0.4832124 1.1924811 0.3657797 1.7400198 0.5577356 0.9889520 $"2" [1] -1.1296642

There is still issues with the reversed y-log scale plot: # Test case A: works as expected plot(10:100,log="y",ylim=c(100,11)) grid() par("yaxp") # Test case B: grid does not have horizontal lines; par("yaxp") is different plot(1:100,log="y",ylim=c(100,10)) grid() par("yaxp") In the second test case, axTicks for the horizontal lines (in

Dear all, I would like to have the same color for the all boxplots from the same panel, but my code below shows the two colors alternating. Thanks! set.seed(42) D1 <- rnorm(200) D2 <- factor(sample(letters[1:2],200,TRUE)) D3 <- factor(sample(letters[3:5],200,TRUE)) DF <- data.frame(x=D1,a=D2,b=D3) print(bwplot(b~x|a,data=DF,col=c("black","black"),

I am experiencing some problems at working with lists at high level. In the following "coef" contains the original DWT coefficients organized in a list. Thorugh applying the following two commands: coef.abs <- lapply(unlist(coef,recursive=FALSE,use.names =TRUE),abs) coef.abs.sorted <- sort(unlist(coef.abs),decreasing=TRUE) I get vector "coef.abs.sorted" containing

On 10/30/2017 03:19 PM, Alexa Catalin wrote: > System: Cenots Linux 6.9 > Application: nut-2.7.5-0.20170613gitb1314c6? [with usb 0.1 from distro] > Device:??Mustek PowerMust 1060 LCD > Comunication log file: dump.txt > > We are looking at the possibility of successful communicating with > this device UPS?Mustek PowerMust 1060 LCD. > > PS: wolfy on the list gives me

System: Cenots Linux 6.9 Application: nut-2.7.5-0.20170613gitb1314c6 [with usb 0.1 from distro] Device: Mustek PowerMust 1060 LCD Comunication log file: dump.txt We are looking at the possibility of successful communicating with this device UPS Mustek PowerMust 1060 LCD. PS: wolfy on the list gives me assistance and i can install any new compiled nut version from sources. Thanks, Catalin.

I'm relatively new to R, and have looked through the documentation and FAQ, and have not been able to find out how I can accomplish something. If someone can point me in the right direction, I would greatly appreciate it. I am doing system performance tests, and end up with large volumes of data that consists of transaction timings, coupled with system timings. The analysis process is

Dear R-list members, I've had a hard time trying to solve a non-linear system (nls) of equations which structure for the equation i, i=1,...,4, is as follows: f_i(d_1,d_2,d_3,d_4)-k_i(l,m,s) = 0 (1) In the expression above, both f_i and k_i are known functions and l, m and s are known constants. I would like to estimate the vector d=(d_1,d_2,d_3,d_4) which is solution

=BE=B4=C6=F4=D5=DF=A3=BA =A1=B0=B5=DA=CB=C4=BD=EC=B1=B1=BE=A9=B3=AF=D1=F4=B9=FA=BC=CA=C9=CC=CE=F1= =BD=DA=A1=B1=BD=AB=D3=DA2003=C4=EA10=D4=C2=D4=DA=B1=B1=BE=A9=BE=D9=B0=EC=A1= =A3=BD=EC=CA=B1=A3=AC=B9=FA=BC=D2=C1=EC=B5=BC=C8=CB=A1=A2=CD=E2=B9=FA=D5=FE= =D2=AA=BA=CD=C0=B4=D7=D4=C3=C0=A1=A2=B5=C2=A1=A2=B7=A8=A1=A2=B0=C4=B4=F3=C0=

Also, just in case - are you in a virtualized environment? Is there (intentional or not) USB pass-through to the UPSes? My hint is, we had cases where a host or guest grabbed the device but it was not apparent from the part of system running NUT. Jim On Fri, Mar 24, 2023, 23:23 Jim Klimov <jimklimov+nut at gmail.com> wrote: > Sounds like some other program is holding the port. Have

Otras variantes con y sin paquetes adicionales... > sapply(split(datIn$Gain, as.factor(datIn$Diet)), mean) d1 d2 d3 280 278 312 > by(datIn$Gain, datIn$Diet, mean) datIn$Diet: d1 [1] 280 -------------------------------------------------------------- datIn$Diet: d2 [1] 278 -------------------------------------------------------------- datIn$Diet: d3 [1] 312 > > library(dplyr) >

Hi Berend, I see you are one of the contributors to the rbecnhmark package. I am sorry that I am bothering you again. I have tried to run your code (slightly tweaked) involving the benchmark function, and I am getting the following error message. What am I doing wrong? Error in benchmark(d1 <- s1(df), d2 <- s2(df), d3 <- s3(df), d4 <- s4(df), : could not find function

I am working with amino acid sequences changing each letter to numbers.I have a data from acf transformation called Zm as shown below. I would like to get Indices D1 to D10 and then create F1 to F10 as indicated below. Is there anyway I can do that in R without typing each of them one by one. For example to get D1 to D3, I have to type D1<-c(Zm[[1]])[1][[1]], D2<-c(Zm[[2]])[1][[1]], and

Sounds like some other program is holding the port. Have you stopped other NUT drivers for the device (e.g. via auto-resuscitating services) before starting this one? Does udev, ugen or similar facility have the configuration to hand off this device to NUT run-time user? (BTW, if you are now testing a custom build - was it configured to use same accounts as pre-packaged variant)? On Fri, Mar 24,

Hello, Another way would be n <- nrow(expand.grid(1:nrow(D1), 1:nrow(D2))) D5 <- data.frame(distance=integer(n),difference=integer(n)) D5[] <- do.call(rbind, lapply(seq_len(nrow(D1)), function(i) t(sapply(seq_len(nrow(D2)), function(j){ c(distance=sqrt(sum((D1[i,1:2]-D2[j,1:2])^2)),difference=(D1[i,3]-D2[j,3])^2) } )))) identical(D3, D5) In my first answer I forgot to say that

- Use MAC16_16 macros instead of (sum += a*b) and unroll a loop by 2. It increase performance when using optimized macros (ex: ARMv5E). A possible side effect of loop unroll is that i don't check for odd length here. - Add NEON version of FIR filter and autocorr -- Aur?lien Zanelli Parrot SA 174, quai de Jemmapes 75010 Paris France -------------- next part -------------- diff --git

Me gusta la respuesta uqe has dado, pero si por ejemplo, alguno de los datos tiene datos faltantes, entonces devuelve NA. He probado con: sapply(split(datos$uno, as.factor(datos$dos)), mean(na.rm=TRUE)) pero da fallo. ¿Cómo se podría hacer para que devolviera además la media obviando los NA y que contara el numero de NA por categoria? > Date: Wed, 28 Oct 2015 00:13:45 +0100 > From: cof

I?m making some progress I believe? I switched off and switched on the UPS further to the laptop having started the daemon and I?m not getting the same error message which is probably due to the server working in the background still locking the device I think the daemon starts but gets into an infinite loop and doesn?t finish and hand back control. Here are the updated traces I got by

Hello, The obvious way would be to preallocate the resulting data.frame, to expand an empty one on each iteration being a time expensive operation. n <- nrow(expand.grid(1:nrow(D1), 1:nrow(D2))) D4 <- data.frame(distance=integer(n),difference=integer(n)) k <- 0 for (i in 1:nrow(D1)){ for (j in 1:nrow(D2)) { k <- k + 1 D4[k, ] <-

Hi List l am interested in developing price model. I have found a research paper related to price model of corn in US market where it has taken demand & supply forces into consideration. Following are the equation: Supply equation: St= a0+a1Pt-1+a2Rt-1+a3St-1+a5D1+a6D2+a7D3+U1 -(1) Where D1,D2,D3=Quarterly Dummy Variables(Since quarterly data are considered) Here, Supply