search for: wildscop

...are also available in SPSS) or 5. 'Best subset selection (as MINITAB's BREG - BEST k command)'. Any response / help / comment / suggestion / idea / web-link / replies will be greatly appreciated. Thanks in advance for your time. _______________________ Mohammad Ehsanul Karim <wildscop at yahoo.com> Institute of Statistical Research and Training University of Dhaka, Dhaka- 1000, Bangladesh

...: t Y 1 0.80 4 0.45 16 0.04 Whatever i could do, is in http://www.angelfire.com/ab5/get5/nonlinear.PDF Any response / help / comment / suggestion / idea / web-link / replies will be greatly appreciated. Thanks in advance for your time. _______________________ Mohammad Ehsanul Karim <wildscop at yahoo.com> Institute of Statistical Research and Training University of Dhaka, Dhaka- 1000, Bangladesh

...t;Run test" (which tests whether a sequence of numbers might be random or not)? If yes, i'd like a numerical illustration of this test. Any response / help / comment / suggestion will be greatly appreciated. Thanks in advance. ------------------------------- Mohammad Ehsanul Karim <wildscop at yahoo.com> Institute of Statistical Research and Training (ISRT) University of Dhaka, Dhaka -1000, Bangladesh ------------------------------- ===== ____________________________ WilDKiD <wildscop at yahoo.com>

...[1] -0.01363636 [1] 0.002207792 [1] 0.001346154 [1] 0.0003787879 -------------------------------------------- Any response / help / comment / suggestion / idea / web-link / replies will be greatly appreciated. Thanks in advance for your time. _______________________ Mohammad Ehsanul Karim <wildscop at yahoo.com> Institute of Statistical Research and Training University of Dhaka, Dhaka- 1000, Bangladesh

...draw the type-A and type-B operating characteristic curve for the plan. Can i do it on R? ----------------- --------------------------- Any comment / suggestion / idea / web-link / replies will be gladly accepted. Thanks for your time. _______________________ Mohammad Ehsanul Karim <wildscop at yahoo.com> Institute of Statistical Research and Training University of Dhaka, Dhaka- 1000, Bangladesh

...ny 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. _____________________________ Mohammad Ehsanul Karim E-mail: wildscop at yahoo.com Web: http://snipurl.com/ehsan ISRT, University of Dhaka, BD

...guage" 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 <wildscop at yahoo.com> Institute of Statistical Research and Training University of Dhaka, Dhaka- 1000, Bangladesh

Dear List, Say, we generate data like this- dat<-rnorm(1000,1,2) hist(dat) How do i make the histogram, say, red (col = 2) before X = dat = 0, and rest say, green (col = 3) beyond X = dat = 0 in R? The resulting histogram could be like this http://ehsan.karim.googlepages.com/histogram.JPG (edited) Thanks in advance. Ehsan http://ehsan.karim.googlepages.com/diaryofastatistician

...outer, 30 Newton-Raphson Variance of random effect= 0.493 Degrees of freedom for terms= 0.5 0.6 1.7 12.1 Likelihood ratio test=47.5 on 14.9 df, p=2.82e-05 n= 76 ###################################################### 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

Dear R-list, I am wondering how to perform a bootstrap in R for the weighted time dependent Cox model? (Andersen?Gill format, with multiple observations from each patients) to obtain the bootstrap standard error of the treatment effect. Below is an example dataset. Would 'censboot' be appropriate to use in this context? Any suggestions/references/direction to R-package will be highly

...How to get the initial value? N.B.: The data file is available at http://www.angelfire.com/ab5/get5/klein.txt Any response / help / comment / suggestion / idea / web-link / replies will be greatly appreciated. Thanks in advance for your time. _______________________ Mohammad Ehsanul Karim <wildscop at yahoo.com> Institute of Statistical Research and Training University of Dhaka, Dhaka- 1000, Bangladesh

Dear List, After including cluster() option the coxreg (from eha package) produces results slightly different than that of coxph (from survival) in the following time-dependent treatment effect calculation (example is used just to make the point). Will appreciate any explaination / comment. cheers, Ehsan ############################ require(survival) require(eha) data(heart) # create weights

...Variance of random effect? 2. What changes do I have to make 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 and Training University of Dhaka # *************************************** "sim.data"<- function(g,m){ set.seed(123) group <- rep(1:m,rep(g,m)) Frailty <- rep(rgamma(m,100,1),rep(g,m)) covariate <- rbinom(g*m,1,.05) stimes...

Dear List, I have a theoretical question related to epidemiological data analysis: If the treatment status (tx = 0,1) changes over time for the patients in a non-randomized cohort, is there a way to estimate the treatment effect? (i.e., after joining the study, some patients may have to wait for a period of time before receiving the treatment, i.e., the situation of patient with id == 2 for the

...ance=1.e-07,method="exhaustive", size.max=ncol(X), nbest=choose(9,5), f.crit=4, plot=T, time=0.02) Any response / help / comment / suggestion / idea / web-link / replies will be greatly appreciated. Thanks in advance for your time. _______________________ Mohammad Ehsanul Karim <wildscop at yahoo.com> Institute of Statistical Research and Training University of Dhaka, Dhaka- 1000, Bangladesh

Dear R users, I'm a new user of R : in fact , i hadn't used R yet. To say the truth , I'm used to with graphical interface of windows environment. so R console reminds me Linux , of which i am afraid of much more. Whatever , i like to learn about the 7 topics below of "R" or "S+4.5" [is there any other similar program like these ?]. Thats why i need

Dear all, So far i could do (in an informal way) to draw a Standardized Resisual plot in the following way- --------------------- >x <- c(104.1, 106.6, 105.5, 107.5, 109.6, 113.3, 115.5, 117.7, 119.9, 122.1, 124.3, 126.5, 128.2) >y <- c(53732, 52912, 57005, 61354, 67682, 71602, 71961, 75309, 82931, 93310, 102161, 103068, 108927) >

Dear [R] members, Can you provide me with the exact URL where can i find ForTran Compilers [ free / commercial / trial / shareware / evaluation copy ] for WINDOWS [ Preferable ] or DOS ; and some source codes of different kind of Statistical calculation making programs made / coded in ForTran ? And how can i use a program (coded in ForTran) in R console ? Can i do it in S+

Dear List, In R we can plot multiple graphs in same page using par(mfrow = c(*,*)). In each plot we can set title using main and sub commands. However, is there any way that we can place an universal title above the set of plots placed in the same page (not individual plot titles, all i need is a title of the whole graph page) as well as sib-titles? Do I need any package to do so? Thank you

Dear List, Suppose, I have some observed and expected frequencies, such as following. I need to draw a graph where plots of observed and expected frequencies are merged into one. ------------------------------------------------ m <- c(1,2,3,4,5,6,7,8,9,10,12,13,17) k <- c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 19) ExpWW <- c(0.309330628803245, 0.213645190887434, 0.147558189649435,