similar to: problem with plot() and R 2.1.0

I have done a fair bit of spectral analysis, and hadn't finished collecting my thoughts for a reply, so hadn't replied yet. What exactly do you mean by normalize? I have not used the functons periodogram or spectrum, however from the description for periodogram it appears that it returns the spectral density, which is already normalized by frequency, so you don't have to worry about

Hola. Tengo un problema con el formato de salida de un objeto "numeric", que los expresa en formato científico o exponencial. > range(total$ImpTotal) [1] 5.66 806907887.10 > Valores <- quantile(total$ImpTotal, c(0, 0.01, 0.02, 0.025, 0.03, 0.05, 0.1, 0.25, 0.5, 0.75, 0.9, 0.95, 0.96, 0.97, 0.975, 0.98, 0.99), na.rm=TRUE) > class(Valores) [1] "numeric" >

Hello all - This code will run, but it bogs down my computer when I run it for finer and finer time increments and more generations. I was wondering if there is a better way to write my loops so that this wouldn't happen. Thanks! -Tyler ################# # Bifurcation diagram # Using Braaksma system of equations # We have however used a Fourier analysis # to get a forcing function

Dear R helpers, Thanks a lot for your earlier guidance esp. Mr Davind Winsemius Sir. However, there seems to be mis-communication from my end corresponding to my requirement. As I had mentioned in my earlier mail, I am dealing with a very large database of borrowers and I had given a part of it in my earlier mail as given below. For a given rating say "A", I needed to have the bad-wise

Hi - I am looking for suggestions for ITSPs for faxing with asterisk 1.8. We are based in the US, so would need an ITSP with US DIDs. #1) We would like to use Fax For Asterisk with asterisk 1.8.4 in order to receive faxes via T.38. Sending faxes is not a requirement. Does anyone have a working asterisk 1.8.4 configuration and ITSP provider that they can recommend? We have been trying T.38

I could use some help understanding how nls parses the formula argument to a model.frame and estimates the model. I am trying to utilize the functionality of the nls formula argument to modify garchFit() to handle other variables in the mean equation besides just an arma(u,v) specification. My nonlinear model is y<-nls(t~a*sin(w*2*pi/365*id+p)+b*id+int,data=t1,

HTB performance improvement Hi all ! i''m looking at the performance of the HTB algorithm/implementation because i would like more packets/sec !! this is the scenario of the performance test: i''m using an embedded system with: SPEED CPU: 399,999 MHz RAM: 128 MByte FLASH: 16 MByte EEPROM: 8Kbyte PROCESSOR MPC8272 a lan to lan 10/100 and in particular we are sending

I've been looking thru manpages and online for help..... I *think* I've got winbind mostly setup (somehow) but it seems to be using the wrong Domain. Some info: I am running RH7.3 Kernel 2.4.18-3 Samba 2.2.7 In my smb.conf file I have: [global] security = Domain workgroup = LMC winbind separator = + winbind uid = 10000-20000 winbind gid = 10000-20000 winbind

All, I am looking for an optimization library that does well on something as chaotic as the Schwefel function: schwefel <- function(x) sum(-x * sin(sqrt(abs(x)))) With these guys, not much luck: > optim(c(1,1), schwefel)$value [1] -7.890603 > optim(c(1,1), schwefel, method="SANN", control=list(maxit=10000))$value [1] -28.02825 > optim(c(1,1), schwefel, lower=c(-500,-500),

All, I am looking for an optimization library that does well on something as chaotic as the Schwefel function: schwefel <- function(x) sum(-x * sin(sqrt(abs(x)))) With these guys, not much luck: > optim(c(1,1), schwefel)$value [1] -7.890603 > optim(c(1,1), schwefel, method="SANN", control=list(maxit=10000))$value [1] -28.02825 > optim(c(1,1), schwefel, lower=c(-500,-500),

Hola Jorge, muchas gracias por tu pronta respuesta, no me di cuenta que el formateo podría causar problemas, envío de nuevo el código sin formatos. La idea básica es para un set de números de columnas (desordenados) y un set de numeros de fila el loop lo que hace es ir a la fila y columna correspondiente de data, tomar el valor y luego hacer la media sobre esos.

Full_Name: Morten Welinder Version: snapshot OS: Submission from: (NULL) (65.213.85.129) Two things are wrong. 1. There is nan test outside IEEE_754. 2. The meat part of the function should really be something like... if (!lower_tail) x = -x; if (fabs (x) > 1) { double temp = atan (1 / x) / M_PI; return (x > 0) ? R_D_Clog (temp) : R_D_val (-temp); } else

Hola de nuevo. Resulta que al actualizar a la versión 0.9.2 de ggplot2 , ahora en mi eje continuo y , los valores se etiquetan en formato científico , cuando antes era normal. Es decir, mis valores máximos están sobre 300000, pero en el eje ahora pone 3e+05 y no he encontrado como cambiarlo. He visto que se puede cambiar utilizando el paquete scales y utilizando por ejemplo , siendo p un

Hi All, I am running Samba3.0 beta 3 installed from RPM's provided on the Samba.org FTP server. My platform is Red Hat 9. I am seeing a lot of errors such as: Aug 6 17:50:18 server3 winbindd[31138]: [2003/08/06 17:50:18, 0] nsswitch/winbindd_sid.c:winbindd_lookupname(103) Aug 6 17:50:18 server3 winbindd[31138]: could not find domain entry for domain AND smbd[27004]: [2003/08/06

Hola, ¿qué tal? Lo que te pasa no es tan raro: set.seed(1234) muestra <- abs(rnorm(100)) sd(muestra) #[1] 0.5811866 muestra.ceros <- c(muestra, rep(0, 100000)) sd(muestra.ceros) #[1] 0.03196273 En una muestra de números positivos, añadir un cero (sobre todo si está lejos de la media) sube la varianza. Si añado otro, posiblemente también. Pero cuando añado muchísimos ceros, la varianza

library(zoo) library(chron) t1 <- chron("1/1/2006", "00:00:00") t2 <- chron("1/31/2006", "23:45:00") deltat <- times("00:15:00") tt <- seq(t1, t2, by = times("00:15:00")) d <- sample(33:700, 2976, replace=TRUE) sin.zoo <- zoo(d,tt) #there are ninety six reading in a day d.max <- rollapply(sin.zoo, width=96, FUN=max)

I am getting ready to experiment with the Tellabs 2752 echo canceller. I have a 255D shelf (and power supply), but am struggling a little on connecting the echo canceller to a PRI. The shelf has 4 25-pair amphenol connectors. The two on the line side are marked "Receive In" and "Send Out". The 2 connectors on the drop side are marked "Send In" and "Receive

Dear all, I think I must be going crazy.... If I do > mean(rpois(1000000, 14)) I get: [1] 13.50420 and again: > mean(rpois(1000000, 14)) [1] 13.49896 > mean(rpois(1000000, 14)) [1] 13.50161 > mean(rpois(1000000, 15)) [1] 14.49250 > mean(rpois(1000000, 15)) [1] 14.49897 > mean(rpois(1000000, 14.5)) [1] 13.99689 > mean(rpois(1000000, 14.5)) [1] 13.9963 This was on a platform

I am currently testing a new Asterisk installation. The server has a T100P connected to a PRI, and about 50 Polycom IP600 phones connected via the local network. Every couple hours, Asterisk randomly stops responding to all calls, both incoming on the PRI and calls from the SIP phones. I'm not sure how or where to start debugging it. When Asterisk stops working I can still connect to

Hola, ¿qué tal? Tengo un df relativamente chiquito: > dim(loans) [1] 15172 22 > object.size(loans) 7488144 bytes Sin embargo, empresas <- split(loans, loans$number_id) length(empresas) #8748 object.size(empresas) # 50643907728 bytes Es decir, el objeto resultante de partir un df de 15000 filas en 8700 dfs ocupa 5GB. ¿Alguna idea de por qué ocurre esto? (Uso R 3.0.1 en un