similar to: R: writing data from one matrix into another with keeping NA's

Hi, I am using the package "demography" from Rob Hyndman for the Lee-Carter-Model. It is an amazing powerful tool but I am struggling with one issue: *I want to compute different quintiles for the cumulative survival probability derived from the Lee-Carter-Forecast (e.g. the 50%-quintile, 75%-quintile and 99%-quintile) for the next 10 years. * I am sure the package possess this

Dear All, I am economist and working on poverty / income inequality. I need descriptive statitics like the ratio of education expentitures between different income quintiles where each household has a different weight. After a bit of google search I found 'Hmisc' and 'quantreg' libraries for weighted quantiles. The problem is that these packages give me only weighted quintiles;

Hi All: Greetings, and best wishes from the festive times here at Kolkata, India -- the time of Durga Puja celebrations. I seek your advice as I try plotting lines for my data. The problem: I have created a dataframe that looks like this (name: myFrame): lowest second third fourth highest significance INAS 0.107 0.115 0.123 0.115 0.166 0.000 MMA 0.091 0.107 0.115

Inline. -- Bert On Tue, Apr 26, 2016 at 2:25 PM, Andr? Luis Neves <andrluis at ualberta.ca> wrote: > Ok. I`m trying to run a Poisson glmm with an observation-level random > intercept. But I`m getting the following error for the 'Baci' variable: > > 'Error: (maxstephalfit) PIRLS step-halvings failed to reduce deviance in > pwrssUpdate'. I guess this message

I want to use function cbind.na at library(qpcR) I install package qpcR and I can use functions such m1 <- pcrfit(reps, 1, 2, l5) > AICc(m1) [1] -102.5843 but when i try cbind.na(1, 1:7) i take message Error: could not find function "cbind.na" Thanks -- View this message in context: http://r.789695.n4.nabble.com/library-qpcR-cbind-na-tp4023339p4023339.html Sent from the

Ok. I`m trying to run a Poisson glmm with an observation-level random intercept. But I`m getting the following error for the 'Baci' variable: 'Error: (maxstephalfit) PIRLS step-halvings failed to reduce deviance in pwrssUpdate'. I guess this message is because the baci variable is not a an integer, and cannot be transformed into an integer as R has a threshold of 2x10^9 even in

Hi, I am trying in R to indicate in which quintile a value of a variable is for every month of my data frame in this case based on volatility. For each month I want to know for each stock if it is in the most volatile quintile of if it is in one of the others. So far I have come up with the following function (see below). Unfortunately, the function only works in some cases and often gives the

I am using R on a Windows XP professional platform. The following code is part of a bigger one CODE press=function(y,x){ library(qpcR) models.press=numeric(0) cat("\n") dep=y print(dep) indep=log(x) print(indep) yfit=dep-PRESS(lm(dep~indep))[[2]] cat("\n yfit\n") print(yfit) yfit.orig=yfit presid=y-yfit.orig press=sum(presid^2)

Dear all, Is there anyway too generate MA plot for 2 qPCR assays (an array of 2x 400). -- View this message in context: http://r.789695.n4.nabble.com/Bioconductor-MA-plot-for-qPCR-array-tp4182805p4182805.html Sent from the R help mailing list archive at Nabble.com.

Hi! Is There a way to manually create an affybatch object from qPCR array data? -- View this message in context: http://r.789695.n4.nabble.com/Creating-affybatch-objects-from-matrix-data-from-qPCR-array-tp3921559p3921559.html Sent from the R help mailing list archive at Nabble.com.

Hi, I have a mysql table with fields userid,track,frequency e.g u1,1,10 u1,2,100 u1,3,110 u1,4,200 u1,5,120 u1,6,130 . u2,1,23 . . where "frequency" is the number of times a music track is played by a "userid" I need to turn my 'frequency' table into a rating table (it's for a recommender system). So, for each user, I need to categorise the frequency of tracks

Hi, Many apologies for sending this twice. I accidentally hit the send button before I finished writing my mail. I am new to R and I hope someone can help me with my problem. I am trying to draw a side by side barplot. There is a main experiment and there are many sub experiments within the main experiment. I would like to draw a bar plot showing the number and type of sub_experiments done for

Hej, Calling newuoa (from the minqa package) makes R crash when the package rgl is loaded first. This however only on certain selected data. The data used for testing (saved to 'bugs.R'): xvals = c(1,2,4,5,7,8,9,10,11,12,14,15,16,18,19,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36) yvals =

> On Apr 10, 2016, at 9:38 AM, David Winsemius <dwinsemius at comcast.net> wrote: > >> >> On Apr 10, 2016, at 3:11 AM, Murray Efford <murray.efford at otago.ac.nz> wrote: >> >> Martin - >> Thanks, but although hatvalues() is useful for calculating PRESS, I can't find anything directly relevant to my question in the influence help pages. After

Hello, Hope that someone could help me plotting longitudinal data below: 7213 3333330001 0.8300 13.05.09 1 1 3333330001 0.8700 09.02.05 NULL 4797 3333330001 0.7700 21.03.07 NULL 2399 3333330001 0.7800 12.04.06 NULL 2400 3333330002 NULL 27.03.06 NULL 7230 3333330002 0.8200 14.05.09 0 2 3333330002 0.8400 09.02.05 NULL 4798 3333330002 0.8700 20.03.07 0 4799 3333330003 0.9000 20.03.07 13 2401

Hi, I'm having a problem getting the panel.average function to work as I expect it to in a lattice plot. I wish to draw lines between the averages of groups of y-values at specific x-values. I have created a dataset below which is similar to my real data. I also show an example of using panel.loess in place of panel.average; it performs in a manner similar to what I want panel.average to do

Hi I get the following error when I try and get the AICc for a gls regression using qpcR: > AICc(gls1) Loading required package: nlme Error in n/(n - p - 1) : 'n' is missing My gls is like this: > gls1 Generalized least squares fit by REML Model: thercarnmax ~ therherbmax Data: NULL Log-restricted-likelihood: 2.328125 Coefficients: (Intercept) therherbmax 1.6441405

I am currently trying to write a program that minimises the amount of work required for “auditable” qPCR data. At the moment I am using an Excel (.csv) spreadsheet as source data that has been transposed to the column format required for R to read. Unfortunately, this means I have* *to manually confirm the whole data set prior to doing any analysis, which is taking a considerable amount of time!

Can you explain why you need them as 'integer', A floating point representation can hold a value upto ~4.5e15 as an "integer" keeping the precision that you might need. Jim Holtman Data Munger Guru What is the problem that you are trying to solve? Tell me what you want to do, not how you want to do it. On Tue, Apr 26, 2016 at 1:11 PM, Andr? Luis Neves <andrluis at

> On Apr 10, 2016, at 3:11 AM, Murray Efford <murray.efford at otago.ac.nz> wrote: > > Martin - > Thanks, but although hatvalues() is useful for calculating PRESS, I can't find anything directly relevant to my question in the influence help pages. After some burrowing in the literature I'm doubting there is an answer out there (PRESS R^2 is always presented in a fairly