similar to: help with writing function

Hello, I am getting output from anova() and summary(aov()) that depends on the order of the factors in the fitted model object, and this has me baffled. I see this dependency with the data.frame below but not with an example (table 6.4) from Montgomery's DOE book. This is with R 1.3.0 on Debian GNU-Linux. Where have I gone wrong? > centerpts run sample CH50mg 1 day1 dev126 0.56 2

I am trying to write a function that will simulate observed counts X and Y as below, I want the function to be able give many replicates, can somebody advise? thanks Oarabile #alpha,n sdx and sdy are constant, and N is a vector of length n. unstructured<-function(n,N,alpha,sdx,sdy){ Vx<-rnorm(n,mean=0,sd=sdx) Vy<-rnorm(n,mean=0,sd=sdy) thetax<-exp(alpha+Vx) thetay<-exp(alpha+Vy)

Tried to look for package GLLMGibbs-(package for fitting mixed models by Gibbs samppling contributed by Myles' and Clayton) in CRAN but not found it, Does anybody know if it is still available? thanks Oarabile

I am trying to determine whether two samples are identical or not. I'm aware that somebody can use the Kolmogorov-Smirnoff test to compare empirical distributions, but since my samples have ties I'm not sure if I'm getting the right p-values for the comparison. Can the Kolmogorov-Smirnoff test be adjusted for the case when ties exists and are there any functions that already

Hello, I recently upgraded my glibc from 2.0.7 to 2.1.1 and now I cannot get R to start. I get R : Copyright 1999, The R Development Core Team Version 0.64.1 (May 8, 1999) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type "?license" or "?licence" for distribution details. R is a collaborative

Hi, I am using R to plot the graph and the problem I am facing with my graph is that I have lots of points concentrated in one area and It is creating a visualization challenge. Is there any commands in R I could use to solve this problem. Even if there is no command, do you know how I could tackle this problem...(I want to separate these points so I could see each of them...) Any help is much

-------------- next part -------------- An embedded message was scrubbed... From: Martin Maechler <maechler at stat.math.ethz.ch> Subject: Re: producing maps Date: Fri, 25 Feb 2005 00:18:59 +0100 Size: 2511 Url: https://stat.ethz.ch/pipermail/r-help/attachments/20050225/75bb1c86/producingmaps.mht

I'm trying to produce empirical bayes estimates based on the lognormal model in disease mapping Is there a way this can be done in R? thanks Oarabile

I want to do kernel smoothing on the standardised mortality rates at specified location(x-y coordinates) and be able to produce contour map. I've tried to use the package called fields in R (function smooth.2d) but the problem is I get negative results after smoothing and I do not understand why? Your contributions will be highly appreciated or if there is any package that can do this I

Hello, I have installed Venables & Ripley's VR_5.3pl025-1 in RHOME/library/VR and added the path to my RLIBS environmental variable. In R (0.62.2) typing "library()" lists these packages, but typing e.g. "library(nnet)" gets "Warning: Package 'nnet' contains no R code". It appears to me the code is there and I have checked the permissions and so

How do you prepare a *.Rd documentation file? What is that format? TIA, Robert Burrows rbb at nebiometrics.com -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To:

Hello, What can I read to learn more about interpreters in general and the R interpreter in particular? As for my level of computer competence, I can maintain my Linux system pretty easily and I can program a bit in C and Python. TIA, Robert Burrows rbb at nebiometrics.com -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Thank you all for your help :-) Much Much Aprreciated........ I tried with the "jitter" command and it worked :-) Cheers, Latha > -----Original Message----- > From: Latha Raja > Sent: Monday, February 07, 2005 3:21 PM > To: 'r-help at stat.math.ethz.ch' > Subject: Need your help with my R plot > Importance: High > > Hi, > > I am using R to

>From Tom: The subject is debugging a program written in the R language,under Windows. (Sorry, but I do not know either the Apple OS or *nix.) A computer program will usually not work on the first try, if only because of the risk of typos. Instead, it must be debugged. Roughly, here is the sequence: (1) One codes a program using the R language, and stores it on the hard drive, using

Dear list members! I've two questions concerning graphic export: a) I want to export my graphics as PostScript files. in this way I use the postscript() function. The tricky part is that they must have a pretended size (7 x 7 cm) and an absoulte font size (10pt). b) how can i (permanent) change the size of the graphic device window? Best regards mirca heli --

Hello - I am a fairly new user of R. I currently have R1.7 installed. I tried to install the KernSmooth package but a message tells me I need R1.8. Is there a way to update to R1.8 while in R1.7 in the same manner that I can use "update.packages", or do I have to install R1.8 in a separate directory and re-install the many packages I have already installed in R1.7? Thanks in advance for

The ices.log says INFO playlist-builtin/playlist_read Currently playing "<dateiname>.mp3" WARN playlist-builtin/playlist_read Corrupt or missing data in file "<dateiname>.mp3" then the thread is shut down <p><p>-------- Original Message -------- Subject: Re: Re-2: [icecast] streaming MP3s? (25-Sep-2003 13:31) From: remco@rc6.org To:

I'm pretty sure it was in the default ices.conf, but I guess I deleted it. How do I setup a 2nd stream? The first plays at 24k, I want same audio for a box on the local lan, so bandwidth is not an issue; so either don't decode/reencode or re-encode at 300+K, or whatever the highest is. Carl K http://www.personnelware.com/carl/resume.html --- >8 ---- List archives:

Dear R-users, I have some very simple data where 1's represent events and zeroes non-events, e.g. temp <- rbind(c(0,1,0,0,1,1,1,0), c(0,0,0,1,0,0,0,0)) For each row in the matrix, I would like to replace a singelton event by a 0. That is, any 1 surrounded by zeroes (010) should be replaced by a zero (000). Sequences of 1's should be left unchanged. So the modified matrix shoud look

Hello Remco, I found these decumentations before but they say nothing about setting the start time and end time of each mp3 file entry in the .m3u file Best Regards - Johann "Remco B. Brink" schrieb: > <quote who="Johann Soukup"> > > > is there any documentation on the .m3u file definition ? > > Google is your friend (tm): > >