search for: aldi

Hi, If I use x<-1:10 y<-rnorm(10,0,1) ### pdf(file="c:\\aldi\\test.pdf") plot(x,y) segments(x,y,x+2,y+2) segments(x,y,x+0.5,y+0.5,col=3) ### dev.off() ### q() Is there a way that I can imbed in the graph plot for each point defined by x to x+0.5 and y to y+0.5 (and colored in green) a different hyperlink? For example point one (together with the g...

...ould like to have it broken based on the indicator. If the line of the plot reaches the last observation of indicator=1 then the line needs to stop; the next line will start at location 22.5 and continue up top 27.3; the next line goes from 35.1 up to 52.9. > x<-read.table(file='c:\\aldi\\testgraph.csv',sep=',',header=T) > x Location indicator otherinfo 1 1.2 1 2.2 2 2.5 1 2.5 3 3.7 1 2.3 4 20.1 2 4.3 5 22.5 2 5.2 6 25.0 2 3.4 7 27.3 2...

Hi, I am trying to install the genetics package on a server with Linux, Fedora. I installed it in a PC and worked fine. In a server since I am not used with R, I am not sure what do I need to change so genetics pkg can see some package dependencies: Any suggestion is appreciated, Aldi Note: genetics expects gregmisc and mvtnorm to be installed already. gregmisc creates gdata etc dependencies. So here is what I have done: Retrieved the base R and installed and compiled it with make under : /users/genetics/aldi/r/R/R-2.0.0 Under it directories: /mvtnorm /gregmisc /genetics are...

Hi all, I tried to follow an online tutorial to run openBUgs but the package BRugs has been removed from R repository. Could someone provide a link for where to download BRugs? Thanks a lot! -- View this message in context: http://www.nabble.com/Where-to-download-BRugs-tp17411410p17411410.html Sent from the R help mailing list archive at Nabble.com.

...the moment that I introduce in the ph1 a missing value in the ph1 as NA x ph1 cm1 0.2345 2 NA 1 1.2222 1 ............. the above transforms into [1] 1.4137000 NA 0.6150000 Question: is there a way I can protect this calculations from the NA values in the ph1 (some kind of: na.rm=T)? TIA, Aldi --

...ected from the vector c(0,1,2). Any good reason? See below the error. > xlrmN1 <- sample(c(0,1,2),400 ,prob=c(0.02 ,0.93 ,0.05 )) Error in sample(length(x), size, replace, prob) : cannot take a sample larger than the population when 'replace = FALSE' Execution halted TIA, Aldi --

...eralities it reports >library(rgenoud) Error in library(rgenoud) : there is no package called 'rgenoud' Execution halted So my question is why under cygwin in a batch mode, it does not find the installed package, which is already installed in my laptop's R? Thank you in advance, Aldi

...; followed by a request for username/password. When I enter a name/password pair I get:- '\\Kate is not accessible. The account is not authorised to login from this account' I can connect OK using smbclient and testparms give me an OK output. Any thoughts gratefully received. Kim Aldis Creative Director The Aldis Animation Co. Ltd. Kim@Aldis.co.uk http://www.aldis.co.uk

...ibution [1999/09/08 17:28:05, 0] lib/fault.c:fault_report(43) =============================================================== [1999/09/08 17:28:05, 0] lib/util.c:smb_panic(2383) PANIC: internal error Nothing immediately obvious in the other logs. ===================================== = Kim Aldis (Kim@Aldis.co.uk) = Creative Director = The Aldis Animation Co. Ltd = http://www.aldis.co.uk = 0171 287 7173 ======================================

...7), c(7, 8, 9)) how can I produce the mean, or the sum of the elements specified in the index iy from the vector Y1? expecting something like this for the sum: Y2 19 19 31 24 5 15 12 I thought lapply function may perform this, but does not work: Y2<-lapply(Y1[iy],sum) Any suggestion? TIA, Aldi --

...ing if there is a way to figure out why in SAS random beta coefficients are 0 vs. in R the beta-s are non zero. The variables of the data are nidl, time, and sub (for subject). Time and nidl are continuous variables. I am applying random coefficients model. Any input is greatly appreciated, Thanks, Aldi 1. mixed model in SAS: ====================== ods output SolutionR = out1.randomnidltest2; proc mixed data = a1 ; class sub ; model nidl = time / solution ; random int time / sub = sub solution; run; ods output close; 2. mixed model in R: ==================== a1<-read.table(file=&quot...

...) cucumbers5<-rnorm(300,mean=0,sd=5) cucumbers6<-rnorm(400,mean=0,sd=6) par(mfrow=c(1,2)) hist(x, main="Left screen OK") boxplot(wheat1,wheat2,tomatos3,tomatos4,cucumbers5,cucumbers6) title ("Right screen: boxplot with plants") Thank you in advance for any suggestions, Aldi --

> Date: Mon, 20 Feb 2006 16:43:55 -0600 > From: Aldi Kraja <aldi at wustl.edu> > > Hi, > Using package gclus in R, I have created some graphs that show the > trends within subgroups of data and correlations among 9 variables (v1-v9). > Being interested for more details on these data I have produced also the > var-covar matrice...

...ries in R? I am thinking if there is any option in R to create a unique id which has its own unique .Rdata, although in the same directory? SAS for example to each batch job it assigns a different ID and a separate temp space, and does not mix it with another job running in parallel. Thanks, Aldi --

...f U2 as (e) which needs an NA in the second position; (d) needs a name bbb8888 as a second column and as a second row and correspondingly in its rows and columns use value 1.0; NO change for (c); adding 0.0 in the second position for (b); and NO change for (a). (see After) Thank you in advance, Aldi Before: ======= > UUU[2] $U2 $U1$scores (e) [1] -1.946707 -57.970488 $U2$cov (d) bbb1378 bbb1379 bbb1378 1.10362564 -0.01222695 bbb1379 -0.01222695 26.88805020 $U2$n (c) [1] 1802 $U1$maf (b) [1] 0.0002774695 0.0077691454 $U2$sey (a) [1] 13.3867 After: ==...

...er blocks, the jobs for the block 19 fail. Error: ReadItem: unknown type 98, perhaps written by later version of R Execution halted even in the one I just mentioned if I execute by hand goes well. Do you know what could be the cause of bsub submission to fail? Any remedy? Thank you in advance, Aldi --

...ile it says: 'qnorm' is based on Wichura's algorithm AS 241 which provides precise results up to about 16 digits. What happen after the 16th digits? If I am running R in a server 64-bit, can that improve the chances that beyond 16th digits to still have precision? Thanks, Aldi --

...ipchart in R (maybe from Dalgaard), but although searching for it I couldn't find it. boxplot(x$earned ~ x$movie) stripchart(x$earned ~ x$movie, vertical=T,data=x, method="jitter", pch=19) Any suggestions are greatly appreciated, otherwise I have to go to points function. Thanks, Aldi --

...12 0.25992 0.06260 4.152 3.35e-05 *** --- Signif. codes: 0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1 Residual standard error: 1.02 on 4244 degrees of freedom Multiple R-squared: 0.01353, Adjusted R-squared: 0.01074 F-statistic: 4.851 on 12 and 4244 DF, p-value: 5.466e-08 Thank you in advance, Aldi P.S. Sorry that I cannot distribute these data for a test. --

...ar-covar matrices for each subgroup, that differ for a treatment on one group (treatment A) vs (non-Treatment A). Is there a software that can compare if two var-covar matrices are statistically the same? Below are a pair of two matrices, from several others. Thank you in advance for any input. Aldi First group var-covar matrix (the data were under treatment a) v1 v2 v3 v4 v5 v6 v7 v8 v9 v1 730.87 3.406 -283.41 -74.68 107.57 -1355.13 -112.46 14.000 5.77...