similar to: How to specify "newdata" in a Cox-Modell with a time dependent interaction term?

Hi all, I am looking at the tutorial/appendix from John Fox on ?Cox Proportional-Hazards Regression for Survival Data? available here: http://cran.r-project.org/doc/contrib/Fox-Companion/appendix-cox-regression.pdf I am particularly interested in modelling survival with time-dependent covariates (Section 4). The data look like this: > Rossi.2[1:50,] start stop arrest.time week arrest fin

Dear List, I?m trying to figure out *memory page sharing***using HVM on Xen 4.0.1, but I?m not really able to find some useful information or sources about this issue. I just found two links porviding these sources: http://knol.google.com/k/learning-grant-tables http://blog.chinaunix.net/space.php?uid=20286427&do=blog&id=109114 Both are not working.. I think they are made for an older

Perhaps I am making this too hard, but how does one regress y(t) on a constant, x(t-1) and y(t-1)? I've tried the manuals and until I get Dalgaard's book (just ordered on Amazon), I am stuck! Thanks to all in advance for your patience and consideration.

Hi there, I have the following problem Four data frames exist: data1 data2 data3 data4 Now I want to write a loop and temporarily store the data1, data2, data3, data4 in a variable called data. I tried the following... for (i in 1:4) { Data <- paste("data",i,sep="") ... .. } but it doesn't function. I think the problem is the definition of the mode of the pasted

Hi there, Does anyone know a method to calculate the number of different patterns in a given data frame. The variables are of polytomous type and not binary (for the latter i found a package called "countpattern" which unfortunately only functions for binary variables). V1 V2 V3 0 3 1 1 2 0 1 2 0 So, in this case, i would like to get "2" as output. Thanks

Simple usage questions that I ought to be able to figure on my own, but can't. 1.I'm able to produce a cross tabulation table showing counts with either table or xtabs. But I want column percentages for interpretation, and it seems stupid to sit there with a calculator figuring marginals and column percentages. How to make R do it after this: > x <- c(1,3,1,3,1,3,1,3,4,4)

I tried your code, but it's giving the following error.. Error in match.fun(FUN) : argument "FUN" is missing, with no default

Hi there, I just don't find the solution on the following problem. :( Suppose I have a dataframe with two predictor variables (x1,x2) and one depend binary variable (y). How is it possible to define a glm object (family="binomial") with a user defined logistic function like p(y) = exp(a + c1*x1 + c2*x2) where c1,c2 are the coefficents which I define. So I would like to do no

Hi User's, This might be a simple question but it is giving me a hard time as I am a new user. I installed R version 2.9.2 (2009-08-24) 1. I just copied a short script from Fox (2002) as a practice and wanted to save it as Rossi.R. The system saved it without complain but when I looked at using a windows explorer it is not *.R file but only Rossi. Why this is happening? 2. the script

Chris here, I'm protecting myself from receiving junk mail. Just this once, click the link below so I can receive your emails. You won't have to do this again. http://spamarrest.com/a2?ZGZmZGRjAmclp3yhL0OfnKA0pl5mLJ1vLF5ipzpj Spam Arrest - Take control of your inbox! http://spamarrest.com/affl?1331107 ------------------------------------------------------------ You are receiving

Hello, I have a multi-year dataset (see below) with date, a data value and a flag for the data value. I want to find the monthly median for each month in this dataset and then plot it. If anyone has suggestions they would be greatly apperciated. It should be noted that there are some dates with no values and they should be removed. Thanks Emily > print ( str(data$flow$daily) )

Hi, I'm new to R and using it for Cox survival analysis. Thanks to this great forum I learned how to compute the HR with its confidence interval. My question would be: Is there any way to get the p-value for a hazard ratio in addition to the confidence interval? Thanks, Thierry -- Thierry Panje Visiting Student Researcher Department of Psychology Stanford

Hi folks, Not sure what this sort of estimation is called. I have a 2-column time-series x(i,t) [with (i=1,2; t=1,...T)], and I want to do the following 'simultaneous' regressions: x(1,t) = (d - 1)(x(1, t-1) - mu(1)) x(2,t) = (d - 1)(x(2, t-1) - mu(2)) And I want to determine the coefficients d, mu(1), mu(2). Note that the d should be the same for both estimations, whereas the

Hi there, I have the following problem: I have a package called "polLCA" which has the following syntax: poLCA(formula, data) and needs the following formula definition: formula <- cbind(V1,V2,V3,...) So far so good. What I tried now was the following: #Get "data" with the "read.table" fuction data <- read.table("d:/ .....) #Select cols to use in the

Dear List, I use Debian Squeeze with Xen dom0 Kernel 2.6.32-5-xen-686, the Xen Hypervisor 4.1 and qemu 4.0. I have installed a HVM domU which contains a Debian Squeeze standard installation. Everytime I start the domU with the ''xm create ..'' command from dom0, the eth0 lan connection gets interrupted and rebuilts. When domU is booted, the dom0 system hangs every minute for

I believe the 911 is a serious issue if one does an asterisk installation in an office. How do you test 911? Won't they arrest you or something for dialing 911 for no reason and talking to one of their agents who could have taken a more important call? On the other hand what an emergency comes up (like someone got seriously injured) and on top of that asterisk crashed all of a sudden

Dear Germano, Thank you for your fast reply, In the above code, *MYData *is the actual data set. Do not we need to convert *MYData to *the dissimilarity matrix using *pam(as.dist(**MYData**), k = 10, diss = TRUE*)* code line?* *Regards.* On Thu, Aug 17, 2017 at 2:58 PM, Germano Rossi <germano.rossi at gmail.com> wrote: > try this > > MYdata <-

Sorry, I never use pam. In the help, you can see that pam require a dataframe OR a dissimilarity matrix. If diss=FALSE then "euclidean" was use.So, I interpret that a matrix of dissimilarity is generated automatically. Problems may be in your data. Indeed pam(ruspini, 4)$diss write a dissimilaty matrix while pam(MYdata,10)$diss wite NULL 2017-08-17 16:03 GMT+02:00 Sema Atasever

Folks- I have developed a package which I am planning on posting to CRAN. I include examples for each of the 40 or so functions in the package. However, since the examples are non-trivial and the package is using MCMC simulation methods, it takes about 20 minutes to run the complete set of examples on a 2.8 Ghz windows PC. The R process reaches a high water mark of about 100 MB of memory as

Dear R user, I want to estimate the parameter of ZETA or/and ZIPF distributions using R, given a series of integer values. Do you know a package (similar to MASS) or a function (similar to fitdistr) I can use to estimate the parameter of these distributions using MLE method? Otherwise do you know a function (which use MLE method to estimate distribution parameters) that allow me to specify a