similar to: Upi

I would like to simulate some SARIMA models, e.g. a SARIMA (1,0,1)(1,0,1)[4] process. I installed the package 'forecast', where the function simulate.Arima should do what I am trying to do. I am not able to understand how it works Could somebody help me with an example? thank you Stefano Sofia AVVISO IMPORTANTE: Questo messaggio di posta elettronica pu? contenere informazioni

Dear list users, if within a function I first define a matrix mymat <- matrix(NA, nrow=m, ncol=n) and somewhere afterwards by mistake mymat <- c(1,2,3) this second command deletes the first one. How can I make sure that within the function mymat will always remain the matrix defined at the beginning, without possibility of changing it throughout the code? Secondly, is it possible to define

Dear R-list users, this question is not strictly related to R, but hopefully somebody will be able to answer. In a schematic way, which is the algorithm to sample from a multivariate normal distribution using the singular value decomposition? thank you for your help Stefano AVVISO IMPORTANTE: Questo messaggio di posta elettronica pu? contenere informazioni confidenziali, pertanto ? destinato

Dear R users, two quick questions about MCMCpoissonChange of the MCMCpack: 1. my data is stored into a vector (e.g. my_data), and as input for MCMCpoissonChange I need either a formula or a data frame. If I simply transform my vector to a data frame with data.frame(my_data), it does not work. If I use the formula my_data~1, then I get the following error: 'You have to prior for lambda (c0 and

Dear R users, I am trying to use the dlm package, and in particular the dlmBSample function. For some reason that I am not able to understand, this function does not work properly and the plot of the result does not make sense, while dlmFilter works perfectly. I think that my_mod is correct, because the output of my_dlmFilter$mod is fine. Where is my mistake or my misunderstanding? This is the

Thank you Rui for your code. I basically understood all your suggestions. I am using an old version of R (version 3.6.3, installed in a server I am not allowed to control), and the new pipe operator does not work. I tried to run your code without the "|>" operator, but I get an error when I use apply. Could you please expand your code without the pipe operator? Thank you again

I would especially look at the CHANUNAVAIL dial status Since it sounds like you are probably qualifying your IAX trunk, that status will be the quickest way to overflow from IAX to TDM. On Sat, May 30, 2015, 11:35 PM Ashwin Surendran < Ashwin.Surendran at now-health.com> wrote: > Hi Matt, > > > > I was a bit concerned on the delay if there might be any when my iax link

Thank you to Ozgur, who gave me the hints for an easy solution. There are the commands that work for my case: mean(my_dataframe$Freq[my_dataframe$month == 5 & my_dataframe$day %in% c(1:7)]) mean(my_dataframe$Freq[my_dataframe$month == 5 & my_dataframe$day %in% c(8:14)]) ... Stefano Sofia PhD This is the original message: Dear list users, I have a data frame as below specified. From

Please keep the mailing list copied in replies. Yes, apparently the map_data function is in the ggplot2 package, not the maps package. -- Sent from my phone. Please excuse my brevity. On August 8, 2017 11:18:35 PM PDT, "??????????? ?????? ?????????" <Maksim.Poltoradnev at uralkali.com> wrote: >Dear Jeff, I figured out what the problem was: the map-data function >does not

Dear R users, I have difficulty to create a list within a list. Example: with > A <- vector(mode="list", 4) I create a list of 4 elements: > A [[1]] NULL [[2]] NULL [[3]] NULL [[4]] NULL In each element of this list I can store, for example, a matrix: A[[1]] <- matrix ... I need each element of A to be a list (all the lists of the same length), and then store matrices

Dear list users, I have a data frame as below specified. From the 1st of May to the 30th of September of several years (e.g. from 2004 to 2011) I have a frequency of accidents. I need the mean of accidents divided by weeks (i.e. the mean of accidents from the 1st to the 7th of May of all the years, from the 8th to the 14th of May, ..., from the 29th to the 31st of May, from the 1st to the 7th of

On 21/03/2018 11:44 AM, Stefano Sofia wrote: > Dear list users, > let me ask you this trivial question. I worked on that for a long time, by now. > Suppose to have a data frame with NAs and to sum some columns with rowSums: > > df <- data.frame(A = runif(10), B = runif(10), C = rnorm(10)) > df[1, ] <- NA > rowSums(df[ , which(names(df) %in%

?s 15:51 de 13/05/2023, Stefano Sofia escreveu: > Dear list users, > > I have to aggregate wind direction data (wd) using a function that requires also a second input variable, wind speed (ws). > > This is the function that I need to use: > > > my_fun <- function(wd1, ws1){ > > u_component <- -ws1*sin(2*pi*wd1/360) > v_component <-

Dear R list users, the singluar value decomposition of a symmetric matrix M is UDV^(T), where U = V. La.svd(M) gives as output three elements: the diagonal of D and the two orthogonal matrices u and vt (which is already the transpose of v). I noticed that the transpose of vt is not exactly u. Why is that? thank you for your attention and your help Stefano AVVISO IMPORTANTE: Questo messaggio di

Should not the result be NULL if you have removed the NA with na.rm=TRUE ? B. > On Mar 21, 2018, at 11:44 AM, Stefano Sofia <stefano.sofia at regione.marche.it> wrote: > > Dear list users, > let me ask you this trivial question. I worked on that for a long time, by now. > Suppose to have a data frame with NAs and to sum some columns with rowSums: > > df <-

Dear list users, I have to aggregate wind direction data (wd) using a function that requires also a second input variable, wind speed (ws). This is the function that I need to use: my_fun <- function(wd1, ws1){ u_component <- -ws1*sin(2*pi*wd1/360) v_component <- -ws1*cos(2*pi*wd1/360) mean_u <- mean(u_component, na.rm=T) mean_v <- mean(v_component, na.rm=T) mean_wd

Dear Jim and dear Enrico, thank you for your replies. Unfortunately your hints didn't solve my problem, and I am getting mad. Can I show you my whole process? I will be as quick as possible. I start from a data frame called Snow of the form year month day hh mm hs 2007 11 19 0 0 0.00 2007 11 19 0 30 0.00 2007 11 19 1 0 0.00 2007 11 19 1 30 0.00 2007 11 19 2 0 0.00 2007 11 19 2 30 0.00 2007 11

Dear R-users, suppose I have three dataframes like these df1: mydate min_temp 31032013 12 01042013 8 02042013 -999 df2: mydate min_temp 31032013 10 01042013 11 02042013 14 df3: mydate min_temp 31032013 4 01042013 3 02042013 5 where -999 means that the temperature data is not available (at the moment I cannot change it to NA because I am not the db administrator); suppose also that oggi is

Hi Matt, I was a bit concerned on the delay if there might be any when my iax link is down? It would be two dial steps right when my iax link is down. But I?m more than happy to try. Many Thanks, Ashwin. From: asterisk-users-bounces at lists.digium.com [mailto:asterisk-users-bounces at lists.digium.com] On Behalf Of Matt Riddell (lists) Sent: 30 May 2015 16:55 To: Asterisk Users Mailing List -

Hi Stefano, As the help page says: "The default for the format methods is "%Y-%m-%d %H:%M:%S" if any element has a time component which is not midnight, and "%Y-%m-%d" otherwise. This is because when the result is printed, it uses the default format. If you want a specified output representation: format(strptime(init_day, format="%Y-%m-%d-%H-%M"),"%Y-%M-%d