[R] Single Title for the Multiple plot page

2007-05-03 Thread Mohammad Ehsanul Karim
is a title of the whole graph page) as well as sib-titles? Do I need any package to do so? Thank you for your time. Mohammad Ehsanul Karim (R - 2.3.1 on windows) Institute of Statistical Research and Training University of Dhaka __ R-help

[R] exclude the unfit data from the iteration

2007-04-24 Thread Mohammad Ehsanul Karim
rogram to continue untill fixed number of times (say, 100) it iterates with good data. Thank you for your time. Mohammad Ehsanul Karim (R - 2.3.1 on windows) Institute of Statistical Research and Training University of Dhaka __ R-help@stat.math.ethz

Re: [R] Extracing "Interval of Time" in seconds in R

2007-04-22 Thread Mohammad Ehsanul Karim
ot; Start time:", begin.times , "\n", "Finish time:", end.times, "\n", "Run time:", run.time, 'secs.\n') Thanks. Mohammad Ehsanul Karim __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

[R] Extracing "Interval of Time" in seconds in R

2007-04-22 Thread Mohammad Ehsanul Karim
POSIXt" objects > as.numeric(format(ending, "%H:%M:%S")) [1] NA How should I proceed? Thanks for your valuable time. Thanks in advance. Mohammad Ehsanul Karim Using Windows Xp, R 2.3.1 Institute of Statistical Research and Training, University of Dhaka _

[R] Approaches of Frailty estimation: coxme vs coxph(...frailty(id, dist='gauss'))

2007-04-20 Thread Mohammad Ehsanul Karim
estimate parameters of a positive stable frailty model) by Wassell et al (1999)? Thanks for your valuable time. Thanks in advance. Mohammad Ehsanul Karim Institute of Statistical Research and Training, University of Dhaka __ R-help@stat.math.ethz.ch m

[R] Hiding "Warning messages" in coxme output

2007-04-20 Thread Mohammad Ehsanul Karim
Dear list, I have been trying to use coxme in R 2.3.1. When I use coxme in the following data sim.fr1, i get "Warning messages: using 'as.environment(NULL)' is deprecated" Why does it occur? How can I hide such warning message, especially when coxme is under a loop? Mo

Re: [R] Extracting approximate Wald test (Chisq) from coxph(..frailty)

2007-04-17 Thread Mohammad Ehsanul Karim
) Mohammad Ehsanul Karim (R - 2.3.1 on windows) wildscop at yahoo dot com Institute of Statistical Research and Training University of Dhaka On Tue, 17 Apr 2007, Mohammad Ehsanul Karim wrote: You _can_ use tmp <- capture.output( print( ) ) and then furthe

[R] Extracting approximate Wald test (Chisq) from coxph(..frailty)

2007-04-17 Thread Mohammad Ehsanul Karim
Thank you for your time. Thanks in advance. Mohammad Ehsanul Karim wildscop at yahoo dot com Institute of Statistical Research and Training University of Dhaka __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-hel

[R] Programming Problem (for loop, random # control, 3 dimentional graph)

2007-04-11 Thread Mohammad Ehsanul Karim
n the next run, where should i put set.seed()? 3. How to plot such plots in R having dimentions G, K and vectorized M? In matlab, mesh(K,G,M) does the trick. scatterplot3d(K,G,M) gives me error message. Thank you for your time. Thanks in advance. Mohammad Ehsanul Karim wildscop at yahoo dot com Ins

[R] Simulation of the Frailty of the Cox PH model

2007-04-08 Thread Mohammad Ehsanul Karim
in the function to run it under correlated frailty model? (may be in kinship package) 3. Is there any program to run frailty by Inverse gaussian or stable family in R? Thank you for your time. Thanks in advance. Mohammad Ehsanul Karim wildscop at yahoo dot com Institute of Statistical Research

[R] frailty in coxph

2006-09-23 Thread Mohammad Ehsanul Karim
s. Mohammad Ehsanul Karim Institute of Statistical Research and Training University of Dhaka [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posti

[R] $theta of frailty in coxph

2006-09-22 Thread Mohammad Ehsanul Karim
of freedom for terms= 1 1 3 0 Likelihood ratio test=17.6 on 5 df, p=0.00342 n= 76 > frailty.object$history[[1]]$theta [1] 5e-09 Thanks in advance. Mohammad Ehsanul Karim Statistical Research and Training, University of Dhaka [[alternative HT

[R] chisq.test using amalgamation automatically (possible ?!?)

2005-06-26 Thread Mohammad Ehsanul Karim
at we don't have to check each elements if they are less than 5 or not - the whole process being automatic, may be by means of programming)? Any hint, help, support, references will be highly appreciated. Thank you for your time. ------ Mohammad Ehsanul Karim

[R] r programming help III

2005-06-24 Thread Mohammad Ehsanul Karim
,p2N1,p3N1,p4N1,p5N1,p6N1, p11^N) # Use of transition probabilities to get expected probabilities d<-p11-p01 P<-p01/(1-d) pN<-(1-P)*pN0+P*pN1 # Expected Probabilities for number of wet days in a week is pN pN # gives the following output: [1] 0.05164856 0.05126159 0.10424380

[R] r programming help II

2005-06-23 Thread Mohammad Ehsanul Karim
6 7 1 6 9 13 8 17 15 7 Now, i know this is a huge problem and time consuming for us, but is there any way we can solve this automatically only by specifying K.JUN data variable (or its particles K.JUN1984,...,K.JUN2002) by means of R programming? Thank you for your time. Any hin

[R] r programming help

2005-06-22 Thread Mohammad Ehsanul Karim
WET[27]*DRY[3]+WET[28]*DRY[2]+WET[29]*DRY[1] ) Thank you for your time. PS: if the above formula can not be seen correctly, one can see details at http://www.angelfire.com/ab5/get5/alt.ren.txt -- Mohammad Ehsanul Karim Web: http://snipurl.com/ehsan I

[R] Calculating moments of a distribution

2005-06-11 Thread Mohammad Ehsanul Karim
Dear All, Is there any existing package or direct function in R / S-plus that calculates moments (raw or central); the usual measures of skewness and kurtosis, and/or pearsonian k criterion? Thank you for your time. -- Mohammad Ehsanul Karim Web: http

[R] Fitting Theoretical Distributions to Daily Rainfall Data

2005-06-08 Thread Mohammad Ehsanul Karim
estion, direction, references, help, replies will be highly appreciated. Thank you for your time. ------ Mohammad Ehsanul Karim Web: http://snipurl.com/ehsan ISRT, University of Dhaka, BD -- __

[R] "Graphics (for goodness of fit)" Question

2005-03-20 Thread Mohammad Ehsanul Karim
Any suggestion, direction, references, help, replies will be highly appreciated. Thank you for your time. Mohammad Ehsanul Karim Web: http://snipurl.com/ehsan Institute of Statistical Reseach and Training University of Dhaka, Dhaka - 1000, Ba

[R] Calculating lengths of runs of 0 or 1 sequences in meteorological data

2005-03-11 Thread Mohammad Ehsanul Karim
should i proceed? Is there any existing program/function/package to solve such problem (seems like the algorithm should be similar to statistical run test)? Any suggestion, direction, references, help, replies will be highly appreciated. Thank you for your time. _________ Moha

Re: [R] Non-Linear Regression (Cobb-Douglas and C.E.S)

2004-04-18 Thread Mohammad Ehsanul Karim
A * (L^(BETA1)) * (K^(BETA2)) data: klein.data ALPHA BETA1 BETA2 0.003120991 0.414100040 1.513546235 residual sum-of-squares: 3128.245 - Thanks in advance for your time and effort - and sorry for my late reply. ___ Mohammad Ehsanul

Re: [R] Non-Linear Regression (Cobb-Douglas and C.E.S)

2004-04-18 Thread Mohammad Ehsanul Karim
s a problem (doing a grid search over infinite parameter space is indeed time consuming). Thanks in advance for your time and effort - and sorry for my late reply. ___________ Mohammad Ehsanul Karim <[EMAIL PROTECTED]> Institute of Statistical Research and Training University of

[R] S Programming

2003-11-17 Thread Mohammad Ehsanul Karim
age" where control stuctures / loops / vectorization and necessery sequences of S programming are presented in an organized form. Any comment / suggestion / idea / web-link / replies will be gladly accepted. Thanks for your time. ___ Mohammad Ehsanul Karim &l