search for: annuities

Hi, Imagine I have the following data frame: > a <- c(1,NA,3) > b <- c(2,NA,NA) > c <- data.frame(cbind(a,b)) > c a b 1 1 2 2 NA NA 3 3 NA I want to delete the second row. If I use na.omit, that would also affect the third row. I tried to use a loop and an ifelse clause with is.na to get R identify that row in which all records are missing, as opposed to the first

Hi: I have tried to find the appropriate clust algorithm for mixed type of data. The suggested way I see is: 1. use daisy to get the dissimilarity matrix 2. use PAM/hclust by providing the dissimilarity matrix, to get the clusters but by following this, when the data set grows bigger say 10,000 rows of data, the dissimilarity matrix will be O(n^2), and out of memory will occur. I am

Hi I'm trying to measure the effect of a policy intervention (Box and Tiao, 1975). This query has to do with the coding of the model rather than with the particulars of my dataset, so I'm not providing the actual dataset (or a simulated one) in this case, apart from some general description. The time series are of length n=34 (annual observations between 1977 and 2010). The policy

...QuotaPeriod 0 PageLimit 0 KLimit 0 </Printer> <Printer NOG01> Info NT Group HPLJ 8000N Location CP2-250 DeviceURI socket://192.168.8.51 State Idle Accepting Yes JobSheets none none QuotaPeriod 0 PageLimit 0 KLimit 0 </Printer> What have I done wrong? Paul Orwig Pacific Life Annuities Newport Beach, California ------------------------------------------------------------------------------ The information in this e-mail and any attachments are for the sole use of the intended recipient and may contain privileged and confidential information. If you are not the intended recipient,...

Hi all, After running kruskal.test I have got results (p<0,005) pointing to reject the hypothesis that the samples were draw from the same population. Howerver when I run the kruskalmc there are no significant differences in any of the multiple comparisons. Is that possible? Some clarification? Thanks, Humber <https://sites.google.com/site/humberandrade> [[alternative HTML version

Dear Luigi, Here's an option: boxplot(x, boxcol="white", whiskcol="white",medcol="white",staplecol="white") abline(h=c(-1,0,1)) grid(NA, 4, lwd = 2) boxplot(x,add=T) Regards, Jos? Iparraguirre -----Original Message----- From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On Behalf Of Luigi Sent: 18 September 2012 09:36

Hello R experts, Is there any possibility to perform exploratory structural equation modeling (ESEM) in R? Which package should I use? Thanks a lot for help, Krzysztof

Hello all, I tried to extract the p-value from adf.test in tseries; however, I got the error message such as > ht=adf.test(list.var$aa) > ht$p-value Error in ht$p - value : non-numeric argument to binary operator > ht Augmented Dickey-Fuller Test data: list.var$aa Dickey-Fuller = -2.3147, Lag order = 4, p-value = 0.4461 alternative hypothesis: stationary > ht$data [1]

Hi Couple of printing issues that I can't get around. Probably relevant conf sections at bottom. Minor one first: In the Add Printer Wizard (the place most users look first), the Samba print server is not listed. Manually typing the \\server\printer works fine, but this is "difficult" for users adding printers on the fly. Adding the printer from \\server\printers works fine too,

Is there any further development in R to identify the missing data Mechanism? Would you please help ? ======= Please find our Email Disclaimer here-->: http://www.ukzn.ac.za/disclaimer ======= [[alternative HTML version deleted]]

Hi everybody, I am quite new to data imputation, but I would like to use the R package ' Amelia II: A Program for Missing Data '. However, its unclear to me how the input for amelia should look like: I have a data frame consisting of numerous coulmns, which represent different experimental conditions, whereby each column has 3 replicates. I want amelia to perform an imputation across

In R I have 3 variables x,y and z that contain: | x y z [1,] 1 12 122 [2,] 1 13 113 [3,] 1 11 111 [4,] 3 13 123 [5,] 4 14 134 [6,] 5 15 155| I am creating a coplot with the following syntax: |given.depth<- co.intervals(z, number=3, overlap=.1) coplot(x~y|z,pch=19,col="red",cex=1.5,given.values=given.depth)| The resulting plot however has has bars in the

Hello, Is there a function in R by which one can run AR model with Bayesian Information Criteria (BIC)? To my knowledge, functions ar and ar.ols could select the order only by AIC. Thanks, Miao [[alternative HTML version deleted]]

Hi all, I want to test whether the MEAN of two different variables, (and different number of observations) are the same. I am trying to use the anova test but it doesn't seem to like that the number of observations are different: a=c(1:5) b=c(1:3) aov_test=aov(a~b) >>Error in model.frame.default(formula = a ~ b, drop.unused.levels = TRUE) : variable lengths differ (found for

Dear all, I have just sent an enquiry but probably I hadn’t expressed myself properly. Could anyone help me with the following? When I run the code on my data I get a boxplot with outliers identified by numbers 200 & 201. However, what I would like is to label these outliers with their corresponding “DATA$num” values of the data frame. In this example, the outliers should be labelled as:

Dear useRs, In a thesis, I found a mention of the use of pairwise deletion in linear regression and GLM (binomial family). The author said that he has used R to do the statistics, but I did not find the option allowing pairwise deletion in both lm and glm functions. Is there somewhere a package allowing that ? Thanks, Arnaud [[alternative HTML version deleted]]

Dear colleagues. I wonder if anybody knows about a procedure in R to calculate the Berger-Parker index and the Brillouin Diversity index. I searched the net but did not find anything good about it. Thanks a lot for any help Mayo, Elvira. [[alternative HTML version deleted]]

Dear useRs, With the following set of command i managed to create the copula density ================================================ >library(MASS) >library(evd) >X<-matrix(sample(1:3000), ncol=2) >U=cbind(rank(X[,1])/(nrow(X)+1),rank(X[,2])/(nrow(X)+1)) > mat1=kde2d(U[,1],U[,2],n=150) image(mat1$x,mat1$y,mat1$z,col= rev(heat.colors(1500)),xlab="",ylab="")

Hi I'm getting an error message with the roots() function in the vars package. Even the example in the help file comes up with an error: > data(Canada) > var.2c <- VAR(Canada, p = 2, type = "const") > roots(var.2c) Error in UseMethod("roots") : no applicable method for 'roots' applied to an object of class "varest" The error is odd, for

I have a data. frame to which you want to change the names to some of their columns. For example: > seba <- data.frame ('constant' = 3, 'numbers' = 1: 10, 'letters' = LETTERS [1:10], otros = 2:11) List their names: > names (Seba) [1] "constant" "numbers" "letters" I want to rename c the column called constant and b the column