search for: diam

...on()? Can you give me some pointers on how I would need to go about this? Thanks and regards, Lorenz -- Lorenz Gygax, Dr. sc. nat. Artificial Intelligence Lab, Department of Information Technology University of Zurich-Irchel, +41-1-635 67 17, gygax at ifi.unizh.ch one.female <- function (X, diam= 1, prop= 1, lw= 1) { ## X a vector giving the position of the symbol c (x, y) ## diam: diameter in the entities on the Y axis ## prop: proportion of x to y axis (so that symbols will be round ## lw: adjustment of line width for large symbols x.ring <- cos (seq (0, 2*pi, len= 500)) *...

...on()? Can you give me some pointers on how I would need to go about this? Thanks and regards, Lorenz -- Lorenz Gygax, Dr. sc. nat. Artificial Intelligence Lab, Department of Information Technology University of Zurich-Irchel, +41-1-635 67 17, gygax at ifi.unizh.ch one.female <- function (X, diam= 1, prop= 1, lw= 1) { ## X a vector giving the position of the symbol c (x, y) ## diam: diameter in the entities on the Y axis ## prop: proportion of x to y axis (so that symbols will be round ## lw: adjustment of line width for large symbols x.ring <- cos (seq (0, 2*pi, len= 500)) *...

...ate the mean; then select "m" rows from group 2 and calculate the mean as well. So far I've been using a for loop for doing it but when it comes to a large data set is rather inefficient. Any hint to vectorize this would be appreciated. toy = data.frame(group = c(rep(1,10),rep(2,8)), diam = c(rnorm(10),rnorm(8))) nsel = c(6,4) smean <- c(0,0) for (i in 1:2) smean[i] <- mean(toy$diam[1:nsel[i]]) Thanks Pedro [[alternative HTML version deleted]]

Here are several ways to get there, but your original loop is fine once it is corrected: > for (i in 1:2) smean[i] <- mean(toy$diam[toy$group==i][1:nsel[i]]) > smean [1] 0.271489 1.117015 Using sapply() to hide the loop: > smean <- sapply(1:2, function(x) mean((toy$diam[toy$group==x])[1:nsel[x]])) > smean [1] 0.271489 1.117015 Or use head() > smean <- sapply(1:2, function(x) mean(head(toy$diam[toy$group==x],...

Dear all, A question related to the following has been asked on R-help before, but I could not find any answer to it. Input will be much appreciated. I got an unexpected sign of the "slope" parameter associated with a covariate (diam) using zeroinfl(). It led me to compare the estimates given by zeroinfl() and hurdle(): The (significant) negative estimate here is surprising, given the biology of the species: > summary(zeroinfl(bnl ~ 1| diam, dist = "poisson", data = valdaekar, EM = TRUE)) Count model coeffici...

I have an ANOVA model with 2 factors "Environment" and "Site", "Diameter" is the response variable. Site should be nested within Environment. Site is also a random factor while Environment is fixed. I can do this analysis using the "aov" function by using these commands: >model<-aov(Diam~Env+Error(Env%in%Site),data=environ) >summary(mode...

Hi All. I am a new R user. Trying to do scatterplot. Not sure how to resolve this error message A<-subset (ErablesGatineau, station=="A") > B<-subset (ErablesGatineau, station=="B") > > plot(diam ~ biom) > abline(lm(diam ~ biom), col = "red") > > goodcases <- !(is.na(diam) | is.na(biom)) > lines(lowess(diam[goodcases] ~ biom[goodcases])) > > library(car) > scatterplot(diam ~ biom, reg.line = lm, smooth = TRUE, + labels = FALSE, boxplots = FA...

Hi Pedro, This may not be much of an improvement, but it was a challenge. selvec<-as.vector(matrix(c(nsel,unlist(by(toy$diam,toy$group,length))-nsel), ncol=2,byrow=TRUE)) TFvec<-rep(c(TRUE,FALSE),length.out=length(selvec)) toynsel<-rep(TFvec,selvec) by(toy[toynsel,]$diam,toy[toynsel,]$group,mean) Jim On 4/3/16, Pedro Mardones <mardones.p at gmail.com> wrote: > Dear all; > > This must have a rather...

...) 1: 1 2 4 3 5: 4 3 1 2 9: 3 4 2 1 13: 2 1 3 4 17: Read 16 items > coin <- factor(treat,label=c("5c","10c","20c","50c")) > time <- scan() 1: 12.0 13.3 13.0 8.9 5: 9.7 7.2 7.3 8.6 9: 9.2 9.0 7.6 7.4 13: 10.1 5.7 7.5 8.9 17: Read 16 items > diam <- c(19,24,28,33) > treat [1] 1 2 4 3 4 3 1 2 3 4 2 1 2 1 3 4 > diam[treat] [1] 19 24 33 28 33 28 19 24 28 33 24 19 24 19 28 33 > d <- diam[treat] > > lsmod <- lm(time ~ spinner + order + poly(d,3)) > anova(lsmod) Analysis of Variance Table Response: time D...

The list of your interest is R-help not R-sig-mac stefano Il giorno 26/gen/06, alle ore 01:20, Sylvain Charlat ha scritto: > Hi, > > Is there any simple way to get histogram for different levels of > factor? > > Say you have the following data set: > > Island Sp.diam > Moorea 1.21 > Moorea 1.27 > Moorea 1.28 > Moorea 1.22 > Moorea 1.28 > Rurutu 1.5 > Rurutu 1.67 > Rurutu 1.75 > Rurutu 1.55 > Rurutu 1.7 > Rurutu 1.55 > Rurutu 1.59 > Rurutu 1.66 > Rurutu 1.7 > > Is there anything better than: > >...

...***************** The program : ----- # Packages # library(spatstat) library(RandomFields) library(adapt) library(MarkedPointProcess) # Program # data(BITOEK) win <- ripras(steigerwald$coord) PointProcess <- ppp(x=steigerwald$coord[,1],y=steigerwald$coord[,2],window=win,marks=steigerwald$diam) plot(PointProcess) rfm.test(coord=steigerwald$coord,steigerwald$diam,MCrepetitions=19) ----- regards. Florent Bonneu Laboratoire de Statistique et Probabilit??s bureau 148 b??t. 1R2 Universit?? Toulouse 3 118 route de Narbonne - 31062 Toulouse cedex 9 bonneu at cict.fr

...now if there is a way of directly calculating the F-ratio of a random effect using the "aov" function. I have 2 factors in my model, "population" which is random and "length" which is the length of female fish within each population. The dependent variable is "diam" which is the average diameter of eggs produced by each female. At present I set up the model like this: >model <- aov(diam~population*length, data=data) >anova(model) then using the output: >ratio=MS-population/MS-interaction >1-pf(ratio, df-population,df-interaction)...

Hello, I am trying to learn how to create a simulated dataset of a forest stand: I must simulate 10 stands, for each of the stands I have been creating a random number of trees. For each tree I should create a random diameter. The problem I have is that I cannot use a matrix because the length of the three items is different. I was trying with list but I have no clue on how to make a dataset out of them... tried "merge".... Anyone has some clever suggestion? THANK you !! Giovanna Giovanna Ottaviani Aalmo...

...and ... The premise is: a DF like so > loremIpsum <- "Lorem ipsum dolor sit amet, consectetur adipiscing elit. Quisque leo ipsum, ultricies scelerisque volutpat non, volutpat et nulla. Curabitur consequat ullamcorper tellus id imperdiet. Duis semper malesuada nulla, blandit lobortis diam fringilla at. Vestibulum nec tellus orci, eu sollicitudin quam. Phasellus sit amet enim diam. Phasellus mattis hendrerit varius. Curabitur ut tristique enim. Lorem ipsum dolor sit amet, consectetur adipiscing elit. Sed convallis, tortor id vehicula facilisis, nunc justo facilisis tellus, sed el...

Dear R-experts---is there a relatively low-pain way to get unicode characters into a plot to a pdf device? pdf(file="cardsymbols.pdf") plot( 0, xlim=c(0,5), ylim=c(0,5), type="n") text(1,1, "&spades;") text(2,2, "&hearts;") text(3,3, "&diams;") text(4,4, "&clubs;") dev.off() (these are the characters that I need the most NOW, but this is a more generic question.) sincerely, /iaw ---- Ivo Welch (ivo.welch at gmail.com)

...of each plant is indicated in the column "species", which is a factor column with many different levels (species). There are multiple records for each species, and there is no pattern to the order in which the species names appear in the data frame. e.g. uniqueID species elev ht diam 1 1 sp2 3.5 1.3 55 2 2 sp2 4.2 0.5 15 3 3 sp3 3.2 1.0 13 4 4 sp65 2.2 2.0 14 5 5 sp43 5.4 5.7 20 6 6 sp2 2.5 4.1 32 7 7 sp12 1.1 0.9 5 8 8 sp3 3.4 3.6 2 I wo...

Is there a easy way to add components like: http://www.dynarch.com/projects/calendar/ http://tinymce.moxiecode.com/example_full.php?example=true in ruby forms, by plugins or whatever? Thank you Rodrigo Dom?nguez Iplan Networs www.iplan.com.ar <http://www.iplan.com.ar/> ? rdominguez@iplan.com.ar www.rorra.com.ar <http://www.rorra.com.ar/> ? rorra@rorra.com.ar --------------

I just realised I hadn't updated the version of PHP Markdown and Extra on the dingus as I pretended I did in my two announcements. It is now fixed: PHP Markdown 1.0.2b7 and PHP Markdown Extra 1.1b1 can now be tested on the dingus. <http://www.michelf.com/projects/php-markdown/dingus/> Sorry for any inconvenience. Michel Fortin michel.fortin@michelf.com http://www.michelf.com/

...(dia,tapply(ALTURA0[,2],ALTURA0[,1],median, na.rm=TRUE)) colnames(ALTURA_0)=c("ddt","reg") #PSEUDO PSEUDO0=subset(d,d$TRAT==f[1],select=c(1,REG[5])) PSEUDO_0=data.frame(dia,tapply(PSEUDO0[,2],PSEUDO0[,1],median, na.rm=TRUE)) colnames(PSEUDO_0)=c("ddt","reg") #DIAM DIAM0=subset(d,d$TRAT==f[1],select=c(1,REG[6])) DIAM_0=data.frame(dia,tapply(DIAM0[,2],DIAM0[,1],median, na.rm=TRUE)) colnames(DIAM_0)=c("ddt","reg") #AF AF0=subset(d,d$TRAT==f[1],select=c(1,REG[7])) AF_0=data.frame(dia,tapply(AF0[,2],AF0[,1],median, na.rm=TRUE)) colnames(AF_0)=...

...ilinglist but it seems I can't access them it the moment, altough I did subscribe and login. Is there anything wrong ? Thanks for your reply ! -- With best regards, Kees Lemmens. ----------------------------------------------------------------------- Delft Institute of Applied Mathematics (DIAM), Faculty Electrical Engineering, Mathematics and Computer Science, Delft University of Technology, Delft, The Netherlands,