similar to: returning multiple values

Dear R users, Can anyone tell me why I might get the error message "the array index is greater than the upper bound for t" in WinBUGS ? t is a vector I have defined. I have checked repeatedly that the array index for t does match with the upper bound till which the loop is run. I am facing this problem while calling the R2WinBUGS package and both in just WinBUGS as

Hello everyone, I need some help regarding calling WinBUGS from R. I have a model for WinBUGS and an R code which calls WinBUGS. On running the code I am being shown the error message : Error in file(con, "wb") : cannot open the connection In addition: Warning messages: 1: In file.create(to[okay]) : cannot create file 'c:/Program

Hi I am trying to install R-igraph, but when I issue the command : install.packages("igraph", lib="~/.R/library") in the R session, I get the following error: Warning in install.packages("igraph", lib = "~/.R/library") : 'lib = "~/.R/library"' is not writable Error in install.packages("igraph", lib =

Hi everyone, I am fitting a simple linear regression model in R. My command is j=lm( Y ~ Sex + begsal + time + int) Call: lm(formula = Y ~ Sex + begsal + time + int) Coefficients: (Intercept) Sex begsal time int 191.916 -241.805 3.969 5.003 3.040 Now I wish to access the values of these coefficients for other purposes

Hi, I need some answers regarding to R. In K-means clustering ,for K=2, How can I find the members of each cluster and What are the coordinates for the cluster centers? Please send me a reply soon. Thank You Sent from Windows Mail [[alternative HTML version deleted]]

Hello, I am trying to calculate the mean value of each row in a data frame (d), I am having troubles and getting errors using the code I have written. Below is a brief example of the code, any thought or suggestions would be great. Thank you for your time, Doug # Example Code: d <- data.frame(st1=c(1,2,3,4), st2=c(2,5,6,7), st3=c(5,5,NA,7), st4=c(6,5,7,8)) avg <-

On 2016-05-25 19:13, Kelly Lesperance wrote: > Hdparm didn?t get far: > > [root at r1k1 ~] # hdparm -tT /dev/sda > > /dev/sda: > Timing cached reads: Alarm clock > [root at r1k1 ~] # Hi Kelly, Try running 'iostat -xdmc 1'. Look for a single drive that has substantially greater await than ~10msec. If all the drives except one are taking 6-8msec, but one is very

Hi, I'm stuck on creating a plot with x tick labels for every Nth tick mark - how is that done? I don't see a simple solution to this in help(plot) or help(par) and what I've tried is not working, eg, the following does not work, although it seems intuitive to me that it should work: x <- seq(-100,1000,25) y <- x * x % find all the x values that are multiples of 100 tmp <-

useR's, I am trying to find a quick way to change some values in a list that are subject to a condition to be NA. Consider the 3x1 matrix: delta <- matrix(c(2.5,2.5,1), nrow = 1) And consider the list named v that has 3 elements > v v[[1]] [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] 4.25 3.25 2.25 1.25 0.25 0.75 1.75 2.75 3.75 4.25

Re: Heteroscedasticity and mgcv. (Collin Lynch) The GAMLSS package can model heterogeneity in the scale parameter (e.g. standard deviastion) [and also heterogeity in skewness and kurtosis parameters].of the response variable distribution. For parametric models a generalized likelihood ratio test can be used to test whether the heterogeity is needed. Alternatively a generalized Akaike

Hello, I have a x,y,z file.Z is not corresponding to a simple analytical function of x and y (see below). I am trying to find the minimum location of this surface and the z value corresponding to this location by a spline interpolation or from a polynomial fit. I tried with the akima package but as the location of the point I am looking for (the minimum) is outside of the convex hull it's

Estimados Este es el test de friedman que se logra asi library(PMCMR) y <- matrix(c( 3.88, 5.64, 5.76, 4.25, 5.91, 4.33, 30.58, 30.14, 16.92, 23.19, 26.74, 10.91, 25.24, 33.52, 25.45, 18.85, 20.45, 26.67, 4.44, 7.94, 4.04, 4.4, 4.23, 4.36, 29.41, 30.72, 32.92, 28.23, 23.35, 12, 38.87, 33.12, 39.15, 28.06, 38.23, 26.65),nrow=6, ncol=6, dimnames=list(1:6,LETTERS[1:6])) print(y)

Dear R helpers I have following dataset rate_number = matrix(c(5, 15, 60, 15, 5, 0, 20, 60, 20,0, 10, 20, 40, 20, 10), nrow = 5, ncol = 3) range_mat = matrix(c(6.25, 6.75, 7.25, 8.75, 9.25, 9.75, 8.5, 9, 9.5, 10.5, 11, 11.5, 4.25, 4.75, 5.25, 5.75, 6.25, 6.75), nrow = 6, ncol = 3) > rate_number        [,1]   [,2]   [,3] [1,]    5     0      10 [2,]   15   20      20 [3,]   60   60      40

Hi >i have this data > X [1] 5.79 1579.52 2323.70 68.85 426.07 110.29 108.29 1067.60 17.05 22.66 [11] 21.02 175.88 139.07 144.12 20.46 43.40 194.90 47.30 7.74 0.40 [21] 82.85 9.88 89.29 215.10 1.75 0.79 15.93 3.91 0.27 0.69 [31] 100.58 27.80 13.95 53.24 0.96 4.15 0.19 0.78 8.01 31.75 [41] 7.35 6.50

Hi all, I've got object *cumsumtab, *where the first row (1,2,3,4,5,6,8,13,14,17,40) are names (actually, categories), and the second row has numeric values. > cumsumtab 1 2 3 4 5 6 8 13 14 17 40 2 4 6 7 9 11 13 14 15 16 17 I'm trying to find a way of finding to which category ( in *names(cumsumtab) * ) belongs a given value (let's say Q1 <- 4.25). So far, I've

Dear all, I would like to convert the first column of a dataframe to a date (original format: year (4 digits) and month (last 2 digits)) >str(dat_FF) 'data.frame': 1022 obs. of 4 variables: $ date : int 192607 192608 192609 192610 192611 192612 192701 192702 192703 192704 ... $ Rm.Rf: num 2.69 2.52 0 -3.06 2.42 2.66 0 4.29 0.51 0.57 ... $ SMB : num -2.49 -1.25 -1.38 -0.2 -0.34

Hi nice people, i am trying to made a double "for" cycle, i wish that for cycle on k activates t times the for cycle on s the first cycle is this: s<-rep(1,10) s [1] 1 1 1 1 1 1 1 1 1 1 > for (n in 2:10) { + s[n]=1+s[n]+rnorm(1,mean=0,sd=1) + } > s [1] 1 4.75 4.86 4.05 4.09 4.56 4.63 4.65 4.12 4.01 now i wish another cycle that activates t times the before cycle gives

Hi all, I''m in the procedure of migrating my old puppet server (puppet-2.6.14-1.el5) with mongrel to a new one (puppet-2.7.14-1.el6.noarch) with passenger. I''m facing a really strange behaviour with the catalog compilation time. I''ve 3 servers: 1.-) Centos 5.8 puppet-2.6.14-1.el5 + rubygem-mongrel-1.0.1-6.el5 Client: Compiled catalog for td055.pic.es in environment

dear members, I am using arfima() from forecast package to model a time series. The following is the code: > LYGH[[202]] [1] 45.40 3.25 6.50 2.15 > arfima(LYGH[[202]]) Error in .fdcov(x, fdf$d, h, nar = nar, nma = nma, hess = hess, fdf.work = fdf$w) : NA/NaN/Inf in foreign function call (arg 5) I tried viewing .fdcov() with the following code:

>>>>> akshay kulkarni >>>>> on Wed, 31 May 2023 20:55:33 +0000 writes: > dear members, > I am using arfima() from forecast package to model a time > series. The following is the code: >> LYGH[[202]] > [1] 45.40 3.25 6.50 2.15 >> arfima(LYGH[[202]]) > Error in .fdcov(x, fdf$d, h, nar = nar, nma = nma,