search for: maiya

There must be a very basic thing I am not getting... I'm working with some entropy functions and the convention is to use log(0)=0. So I wrote a function: llog<-function(x){ if (x ==0) 0 else log(x) } which seems to work fine for individual numbers e.g. >llog(0/2) [1] 0 but if I try whole vectors or tables: p<-c(4,3,1,0) q<-c(2,2,2,2) llog(p/q) I get this: [1] 0.6931472

Hi, I have downloaded openssh(http://www.openssh.com/) on my Windows host. Just wanted to know if openssh allows executing a batch script on Windows host. Batch script contains mapping of shares which is required for my scripts. My command run from hpux: $ ssh <windows server> myscript The myscript requires shares to be mapped on login via

OK, this is really weird! here's an example code: t1<-c(1,2,3,4) t2<-c(4,2,4,2) plot(t1~t2, xlab="exp1", ylab="exp2") dev.copy2eps(file="test.eps") that all seems fine... until you look at the eps file created, where for some weird reason, if you scroll down to the end, the code reads: /Font1 findfont 12 s 0 setgray 214.02 18.72 (e) 0 ta -0.360 (xp1) tb

Hi! (a complete newby, but will not give up easily!) I was wondering if there is any way to decouple the axis and tick mark widths? As I understand they are both controlled by the lwd setting, and cannot be controlled independently? For example I might want to create major and minor ticks, which I now know how to do by superimposing two axes with different at settings, but what if I also wanted

Hi! This is going to be a real newbie question, but I can't figure it out. I'm trying to plot densities of various functions of chi-square. A simple chi-square plot I can do with dchisq(). But e.g. chi.sq/degrees of freedom I only know how to do using density(rchisq()/df). For example: plot(1, type="n", xlab="", ylab="", xlim=c(0,2), ylim=c(0,7)) for (i

Hi there! I'm not sure I can create a minimal example of my problem, so I'm linking to a minimal .RData file that has only two objects: obs and exp, each is a 6x9 matrix. http://dl.dropbox.com/u/10364753/test.RData link to dropbox file (I hope this is acceptable mailing list etiquette!) Here's what happens: > obs[1, 1] [1] 118 > exp[1, 1] [1] 118 > obs[1, 1]-exp[1, 1] [1]

I am fairly new to log-linear modelling, so as opposed to trying to fit modells, I am still trying to figure out how it actually works - hence I am looking at the interpretation of parameters. Now it seems most people skip this part and go directly to measuring model fit, so I am finding very few references to actual parameters, and am of course clear on the fact that their choice is irelevant for

I'm trying to import a table into R the file is about 700MB. Here's my first try: > DD<-read.table("01uklicsam-20070301.dat",header=TRUE) Error: cannot allocate vector of size 15.6 Mb In addition: Warning messages: 1: In scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : Reached total allocation of 1535Mb: see help(memory.size) 2: In scan(file, what,

Hi everyone! This is a ridiculously simple problem, I just can't seem to find the solution! All I need is something equivalent to sum(is.na(x)) but instead of counting missing values, to count empty cells (with a value of 0). A naive attempt with is.empty didn't work :) Thanks! Maja Oh, and if the proposed solution would be to make all the empty cells into missing cells, that is

i appologise for the trivialness of this post - but i've been searching the forum wothout luck - probably simply because it's late and my brain is starting to go.. i have a frequency table as a matrix: orig<-matrix(c(40,5,30,25), c(2,2)) orig [,1] [,2] [1,] 40 30 [2,] 5 25 i basically need a random sample say 10 from 100: [,1] [,2] [1,] 5 2 [2,] 0 3 i

Hi! I have a dataset with three columns -the first two refer to x and y coordinates, the last one are odds ratios. I'd like to plot the data with x and y coordinates and the odds ratio shown as a fourfold plot, which I prefer to do using the stars function. Unfortunately the stars option in symbols is not as cool as the stars function on its own, and now i can't figure out how to do it!

...no-latex? (Jim Lemon) > 37. Survreg function for loglogistic hazard estimation (chenjiakai) > 38. Re: RPostgreSQL segfault with LEFT JOIN (Joe Conway) > 39. Re: Journal Articles that Have Used R (Jonathan Baron) > 40. stars (as fourfold plots) in plot (symbols don't work) (maiya) > 41. extract rows having negative values (Manisha Brahmachary) > 42. Convert a lis to matrix (Manisha Brahmachary) > 43. Re: extract rows having negative values (Bill.Venables at csiro.au) > 44. Re: EBImage not loading (Gregoire Pau) > > > ----------------------------...