search for: ziberna

Hi All, I have conducted the following survival analysis which appears to be OK (thanks BRipley for solving my earlier problem). > surv.mod1 <- survreg( Surv(timep1, relall6)~randgrpc, data=Dataset, dist="weibull", scale = 1) > summary(surv.mod1) Call: survreg(formula = Surv(timep1, relall6) ~ randgrpc, data = Dataset, dist = "weibull", scale = 1)

...#of course, the "pattern" in this case is wrong I know I can do text[grep(pattern="MY", x=text)][grep(pattern="TARGET", x=text[grep(pattern="MY",x=text)])] However I hope there exists a more elegant way. Thanks in advance for any suggestions! Best, Ales Ziberna

...ages) from which I had taken the code or only to the references they use. What about regarding code that was sent to the list, usually as a response to one of my problems. I assume that in this case it is best to consult the author? Any comments and opinions are very welcomed! Best regards, Ales Ziberna

...function that should pass argument "m" to the next function, however it does not! Please have a look at the function below that shows a problem and tell me what I am missing. As you can see, the "blocks" argument is passed corectly, while "m" is not. Best, Ales Ziberna opt.par.new<-function( #function for optimizig partition in blockmodeling M, #matrix clu, #initial partition maxiter=50, #maximum number of iterations trace.iter=FALSE, #save a result of each iteration or only the best (minimal error) switch.names=is.null(BLOCKS), #should partitions...

...get "ok" a #I do not want "a" in global enviorment, so here I should get #Error: Object "a" not found ff<-function()f() #here I also want to get "ok" - it should not matter if the parent fuction has any parents Thank you in advance for suggestions! Ales Ziberna

...(Mb) Ncells 165659 4.5 350000 9.4 350000 9.4 Vcells 61135 0.5 786432 6.0 283043 2.2 [1] 1.14 0.00 1.40 NA NA This does (usually) not happen in my case, since most of functions I consider are more complex. Thanks in advance for any suggestions! Best regards, Ales Ziberna

Does anybody know how to produce a East-European character "??" - "c with a v-shaped hat " in R (in text or plot)? I do know how to get "??,??" - "s,z, with a v-shaped hat", but not this one. Thanks in advance for any suggestions, Ales Ziberna P.S.: I am using Windows XP and R version 1.9.1 (details below) $platform [1] "i386-pc-mingw32" $arch [1] "i386" $os [1] "mingw32" $system [1] "i386, mingw32" $status [1] "" $major [1] "1" $minor [1] "9.1&...

...and left of plotting region), but non of the functions can plot there (I tried 'legend', 'text', 'points', 'lines', 'mtext', 'axes',...). I would be very grateful if someone could help me put something in that corrner, with any fuction. Thanks, Ales Ziberna

...trix to a elements of a list. #For example, if I have mat<-matrix(1:100,ncol=5, nrow=20) #I can do: list<-apply(mat,1,list) list #however this is not quite what I want. To get what I want, I have to do: list<-lapply(list,function(x)x[[1]]) list Is there a faster way? Best regards, Ales Ziberna

...classes are (I hope I have not missed anyone): 1,1,1,2 (this can be read as {1,2,3},{4}) 1,1,2,1 1,2,1,1 2,1,1,1 1,1,2,2 1,2,1,2 1,2,2,1 The partitions 1,1,2,2 and 2,2,1,1 are the same and are therefore not two unique partitions. Thank you in advance for any suggestions. Best, Ales Ziberna

...ich I would like to build on WinXP. I have all suggested tools/programs for bulding R packages on Windows (except latex). What is the best way of using these subrutines? Does sombody mybe know any translation tools for converting Free-form to fixed-form Fortran? Thanks for any suggestions, Ales Ziberna

Hello, it seems that the main results of survival analysis with package survival are shown only as side effects of the print method. If I compute e.g. a Kaplan-Meier estimate by > km.survdur<-survfit(s.survdur) then I can simply print the results by > km.survdur Call: survfit(formula = s.survdur) n events median 0.95LCL 0.95UCL 100.0 58.0 46.8 41.0 79.3 Is

...d! C:\Ales\Statistika\Blocno modeliranje\dr\blockmodeling>R CMD check * checking for working latex ...Error: environment variable TMPDIR not set (or s et to unusable value) and no default available. at D:\PROGRA~1\R\rw2011\share\perl/R/Utils.pm line 72 Thank you in advance for any help! Ales Ziberna P.S.: I am runing R 2.1.1 on Win XP, SP2. I installed rtools, mingw, perl as suggested in the manuals.Here are some details about R: platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 1.1 year 2005 month 06 day 20 language R

Hi, I am a newbie with R, so I hope my question isn't too stupid. I am trying to generate dose-response curves using the "drfit" package. I have formatted my CSV files to the correct format, and have no trouble running drfit to get a summary of my data. The problem is that when I try to use "drplot" to graph my data I get an error. The message is: Error in

...le to have this thing plotted with the same look as in R 1.9.1? I apologize if the instruction are in NEWS, since I am not sure how to interpret exactly this section (especially : "just add parentheses -- or braces in the case of plotmath". Thank you in advance for any replies! Ales Ziberna I am running binary version of R 2.0.1 (and R 1.9.1) on Windows XP Here are the R 2.0.1 version details: $platform [1] "i386-pc-mingw32" $arch [1] "i386" $os [1] "mingw32" $system [1] "i386, mingw32" $status [1] "" $major...

...res the arguments for the second one. What is the best way to put all resoults of the first one into the second one? I tried attach, however the object in the "main" enviorment have a priority over the ones in list. An example is at the end. Thanks in advance for any suggestions! Ales Ziberna For example - I would like to use just "a" instead of "list$a" in fuction "second" first<-function(a,b){ if(length(a)!=1) a <- a[1] if(length(b)!=1) b <- b[1] list(a=a,b=b) } second<-function(list,c){ list$a + list$b + c } a<-c(2,3)...

...= 2. Therefore, we can say we have units 1, 2 and 3. I believe all possible samples are : {1,2},{2,3} and {1,3}. I would like to emphasize that I am not looking for the number of different samples, but a list of all possible samples. Thank you in advance for any suggestions. Best, Ales Ziberna

Several times my packages have fallen foul of the latex check at the uploading-to-CRAN stage, creating extra work for Kurt Hornik. I've never bothered trying to get latex working on my own packages (happy with vanilla help) but am told that "latex is the closest thing to a syntax checker for Rd files" and that I should set it up. Now, there is already a latex on my (Windows XP)

See Also http://finzi.psych.upenn.edu/R/library/caTools/html/combs.html Jarek Tuszynski -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch] Sent: Thursday, December 08, 2005 11:19 AM To: Ales Ziberna Cc: R-help Subject: Re: [R] Finding all possible partitions of N units into k classe On 08-Dec-05 Ales Ziberna wrote: > Dear useRs! > > I would like to generate a list of all possible (unique) > partitions of N units into k classes. For example, all possible > partitions of 4 units...

...ements only in this two areas: 1. If I am running two instances of R. 2. If I'm running several other programs on the computer beside R, the programs and R would run faster, since they would not "compete" for processor time (so much) Thanks in advance for any useful suggestions, Ales Ziberna P.S.: Useful posts on the list follow: It depends on the usage pattern. If you run multiple CPU-bound processes in parallel without too much coordination (parallel make is a good example, simulations another), then you get close to double up from a dual. For a single R process, you can get some...