search for: mahbub

Hello list: How to extract the value of "Rho" from a gnls() object. I am using gnls() function similar to res <- gnls(y~SSmicmen(),correlation=corCompSymm(form~1|b),data=dat) Thanks in advance, Mahbub. -- Mahbub Latif School of Mathematical Sciences Queen Mary, University of London United Kingdom [[alternative HTML version deleted]]

...configuration failed for package 'RODBC' ** Removing '/usr/local/lib/R/site-library/RODBC' R > version platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor 9.0 year 2004 month 04 day 12 language R Thanks, Mahbub.

Dear friends There is a list of arrays comprising different no of rows and columns even sometimes NULL, such as [[2]] given below. How can we ignore [[2]] or others like this in the complete list. Any help in this regard is needed. Thanks [[1]] [,1] [,2] [1,] 3 1 [2,] 3 1 [3,] 3 1 [[2]] NULL [[3]] [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 3 1

...222) > x <- sample(x=1:100000, size=10000, replace=T) > max(x) [1] 99992 > summary(x) Min. 1st Qu. Median Mean 3rd Qu. Max. 6 24510 49950 50020 75040 99990 Which value is the correct one? I had this problem while I was working with microarry data. Thanks. Mahbub

Hi, I want to use (similar to) las options in lattice (bwplot) plot. Actually I want to have x-axis labels as vertical instead of default horizontal. Thanks in advance for your help. Mahbub

Hi there, I have four panels in a lattice bwplot. I want to have two different ylim for the panels, for example panels [1,1] and [1,2] with ylim=c(0,200) and panels [2,1] and [2,2] with ylim=c(0,100). Thanks for help in advance. Mahbub.

Hi, How can I use portrait/landscape option in lattice bwplot? Is there any option in trellis.device where I can define this? Thanks in advance, Mahbub.

Hi, I was trying to install XML package in a linux (dabian) machine and got the following error message. I am not sure whether there is any error in my linux installion. I appreciate suggestions to install XML properly in this machine. Thanks in advance. Mahbub. ######Errors#### $ R CMD INSTALL -l /usr/local/lib/R/library/ XML_0.93-4.tar.gz * Installing *source* package 'XML' ... creating cache ./config.cache checking for gcc... gcc checking whether the C compiler (gcc ) works... yes checking whether the C compiler (gcc ) is a cross-compiler....

...is that it does not show separate grid for each category. p <- qplot(category,value, data=dat.m, geom="bar", ylab="Percentage of eco change",stat="identity",fill=variable, position="dodge") p + coord_flip() Can you help me generate my vision plot? -- Mahbub Majumder Graduate Student Dept. of Statistics Iowa State University [[alternative HTML version deleted]]

Hi, I am looking for C/C++ codes for computing generalized inverse of a matrix. Can anyone help me in this regard? Thanks, Mahbub.

Dear R-helpers, i have a question on how to vectorize this problem: i have a dataframe: tester <- data.frame(groups=c("A","A","B","B","C","C"), one=c(1,1,2,2,3,3), two=c(6,6,7,7,8,8)) # i split it into a list tester.L <- split(tester, tester$groups) # And want to keep only the first item in each: lapply(tester.L, function(x) x

Hi, I use lme to fit models like R> res1 <- lme(y~A+B, data=mydata, random=~1|subject) R> res2 <- lme(y~B+A, data=mydata, random=~1|subject) (only difference between these two models are the sequence in which the indep variables are written in formula) where y is continuous and A, B, and subject are factors. To get ANOVA table I used R> anova(res1) R> anova(res2) and found

...irst panel and (2,-2.5) in the second panel. I need help on this.</FONT></P> <P><FONT size=4></FONT>&nbsp;</P> <P><FONT size=4>Thanks.</FONT></P> <P><FONT size=4></FONT>&nbsp;</P> <P><FONT size=4>Mahbub.</P> <P><BR></P></FONT> <P><FONT size=4>&nbsp;</P></FONT> <P>&nbsp;</P>

Dear User, what's the correct way to obtain a multiple return from a function? for example creating the simple function: somma <- function (a, b) { c <- a+b return (a, b, c) } when I call it, it runs but returns the following output: > somma(5, 7) $a [1] 5 $b [1] 7 $c [1] 12 Warning message: return multi-argomento sono deprecati in: return(a, b, c) i.e. multi-return is

Hello, I have written a function that demonstrates the CLT by generating samples following the exponential distribution, calculating the means, plotting the histogram, and drawing the limiting normal curve as an overlay. I have the title of each histogram state the sample size and rate (1/theta) for the exponential (the output is actually 4 histograms), but I can't get the greek letter theta

Hi, I have a dataframe (say myData) and want to get a list (say myList) that contains a matrix for each row of the dataframe myData. These matrices are calculated based on the corresponding row of myData. Using a for()-loop to do this is very slow. Thus, I tried to use apply(). However, afaik apply() does only return a list if the matrices have different dimensions, while my matrices have

Dear all, Is it possible to somshow plot text as italic AND bold. I tried font=c(2,3) in text(), but it doesn't work. It seems like the latter value is used. Thanks in advance! Sincerely, Tord ----------------------------------------------------------------------- Tord Sn?ll Avd. f v?xtekologi, Evolutionsbiologiskt centrum, Uppsala universitet Dept. of Plant Ecology, Evolutionary Biology

Hi, there. Does anybody know how to plot a smooth density plot for some data simulated from certain distribution? Thanks. Yulei $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Yulei He 1586 Murfin Ave. Apt 37 Ann Arbor, MI 48105-3135 yuleih at umich.edu 734-647-0305(H) 734-763-0421(O) 734-763-0427(O) 734-764-8263(fax) $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

Hi, I've tried to compile 3.0.22 on HPUX 10.20 and get the following error msg: No options are used when running 'configure'. checking configure summary ERROR: No locking available. Running Samba would be unsafe configure: error: summary failure. Aborting config ****** The server details: [qhus029@HP-UX]>uname -a HP-UX qhus029 B.10.20 A 9000/898 1684683361 two-user license

hi: apologies for taking up everyone's time. my problem is probably documented somewhere, but I again cannot find it. (which reminds me: I cannot find a search engine that allows me to search the archives of this very useful mailing list.) * it seems that lines(lowess()) fails to plot certain line segments. (and, what does it do at the x-min and x-max of a data set?) Rather than