similar to: from 3 numeric variables to a string

Howdee, My question appears at #6 below: 1. I want to model the growth of each of a large number of individuals using a 4-parameter logistic growth curve. 2. nlme does not converge with the random structure that I want to use. 3. nlsList does not converge for some individuals. 4. I decided to go around nlsList using: t(sapply(split(data, list(data$id)), function(subd){coef(nls(mass ~

Howdee, *** I know that the lmList() function exists, yet I don't want to use it. *** Would anyone be kind enough to tell how I can apply the function lm() to each level of a given factor so to obtain the intercept and slope for each factor level within a matrix? For instance, suppose a dataframe containing 3 variables: id, x and y. I want to compute the function lm() for each level

Howdee, I am able to fit a 4-parameter logistic growth curve to a dataset which comprise many individuals (using R v. 2.3.1). Yet, if I want to obtain the parameters for each individual (i.e., for each 'id') using nlsList, then I obtain an Error message which I have trouble interpreting. Any advice as to how I can solve this problem? Thanks for your time, Marc > reg <-nls(mass ~

Is there a way to apply a function, say cor(), to each combination of some number of variables, and this, without using loops? For example, I have day, hour, var1 and var2. How could I compute cor(var1,var2) for each day*hour combination and obtain a matrix with day, hour and the cor value for each combination? Thanks for your time, Marc =================== Marc Bélisle Professeur adjoint

Howdee, One of my student spotted something I can't explain: a probability >1 vs a normal probability density function. > dnorm(x=1, mean=1, sd=0.4) [1] 0.9973557 > dnorm(x=1, mean=1, sd=0.39) [1] 1.022929 > dnorm(x=1, mean=1, sd=0.3) [1] 1.329808 > dnorm(x=1, mean=1, sd=0.1) [1] 3.989423 > dnorm(x=1, mean=1, sd=0.01) [1] 39.89423 > dnorm(x=1, mean=1, sd=0.001) [1]

How do one dertermine the optimal number of rescaling cycles to use when performing DCA using the decorana function (library vegan)? Thanks for your time, Marc =================== Marc Bélisle Professeur adjoint Chaire de recherche du Canada en écologie spatiale et en écologie du paysage Département de biologie Université de Sherbrooke 2500 Boul. de l'Université Sherbrooke, Québec J1K 2R1

Hello, I apologize if this seems like an obvious question, but I have been looking everywhere and have yet to find an answer. I am doing a multinomial regression with multinom() in the nnet package. I have a 3 level ordered response (ordered()) variable and 4 predictors, 3 of which are numerical and one which is an ordered factor (also ordered()) with 5 levels (a, b, c, d, e). My question is in

Dear R experts, For quite some time I have been trying to solve a mistery of reading a seemingly trouble-free text file. The data is temperature reconstruction arranged as a huge grid, preceded by seven "header lines" (which you see better if file is opened in Firefox or Chrome). This is the data (gridded temperature reconstruction)

I can't seem to be able to change the font size in an effect display. I've tried the following: > par(cex.lab=4) > plot(effect ("alti",reg8), ylab="detection probability") and > plot(effect ("alti",reg8), ylab="detection probability", cex=4) but nothing changes. Can anyone help me? thanks. Emilie Berthiaume graduate student

Hi, I'm just getting started with R and I'm having problems with some simple operations: I want to get the the sum of the column "SStot" for each year using the function by. The data set is named "SS". I've tried this: by (SS, year, sum(SStot)) and it's not working. Is it because there's a different number of rows for each year? How else can I do

Hi list! I did the following regression: reg1 <- glm(alti~sp + ovent + vivent + nuage, family=gaussian, data=meteo1) I was interested in knowing the effect of the species (sp) in reg1 and so I used the function «effect»: effect.sp <- effect ("sp", reg1, se=TRUE) with this output: sp AK BW NH OS RT SS 2.730101 2.885363 2.753774 2.750311

One possibility ist boxplot(sapply(ListOfNames, get, env = .GlobalEnv)) Hope that this helps Thomas -----Urspr??ngliche Nachricht----- Von: Lefebure Tristan [mailto:Tristan.Lefebure at univ-lyon1.fr] Gesendet: Mittwoch, 7. Juli 2004 10:34 An: r-help at stat.math.ethz.ch Betreff: [R] boxplot a list of objects Hi list, #Imagine we have vectors of different length (in practice 100 vectors):

Hi I use RSvgDevice to output plot, and modify them using svg editor (inkscape or sodipodi on Linux). Some month ago, results were perfect. I did exactly the same analysis today on the same data, and unfortunatly the results are different. While looking to the svg file, it seems that all information concerning "fill" and "stroke" of objects are lost. The consequence is that

Greetings, We wanted to announce v. 1.4.0 of the R package 'KScorrect', which carries out the Lilliefors correction to the Kolmogorov-Smirnoff (K-S) test for use in (one-sample) goodness-of-fit tests. Aside from several minor changes, the biggest change is that the Monte Carlo algorithm now supports parallel implementation, using the platform-independent 'doParallel' and

Greetings, We wanted to announce v. 1.4.0 of the R package 'KScorrect', which carries out the Lilliefors correction to the Kolmogorov-Smirnoff (K-S) test for use in (one-sample) goodness-of-fit tests. Aside from several minor changes, the biggest change is that the Monte Carlo algorithm now supports parallel implementation, using the platform-independent 'doParallel' and

Hello, i' m looking for a function to compute geary statistic for spatial correlation. Thanks. -- St?phane DRAY --------------------------------------------------------------- Biom?trie et Biologie ?volutive - Equipe "?cologie Statistique" Universite Lyon 1 - Bat 711 - 69622 Villeurbanne CEDEX - France Tel : 04 72 43 27 56 Fax : 04 78 89 27 19 04 72 43 27 57 E-mail

Hi, Just something I don't understand: data <- data.frame(V1=c(1:12),F1=c(rep("a",4),rep("b",4),rep("c",4))) data_ac <- data[which(data$F1 !="b"), ] levels(data_ac$F1) Why the level "b" is always present ? thanks Tristan, R 2.0.1 for Linux Fedora 3 -- ------------------------------------------------------------ Tristan

Hi, can anyone help me with this problem i''ve got! I have a problem with routing on my Mandrake 7.2 Linux Box i was wondering if you could help me with. I live in Australia and as broadband access isn''t very available at the moment, i have 2 modems to dialup to my isp (dynamic ips assigned) to increase my bandwidth at an affordable level. What currently happens is i have a

Hi, I tried to install R packages required for the GRASS-R interface by using the following command (copied from http://www.geog.uni-hannover.de/grass/statsgrass/grass6_r_install.html): /install.packages (c("sp", "spgrass6","rgdal","maptools"), dependencies =TRUE)/ rgdal package was installed successfully; There were problems to find 'sp' and

Hello, I am trying to play with RODBC library and Excel Files. In my file (doubs.xls) there are 2 spreadsheets: > library(RODBC) > connection<-odbcConnect("Excel Files") > sqlTables(connection) TABLE_CAT TABLE_SCHEM TABLE_NAME TABLE_TYPE REMARKS 1 F:\\Th?se\\R\\Doubs NA Faune$ SYSTEM TABLE NA 2 F:\\Th?se\\R\\Doubs NA Milieu$