similar to: Format Data Issue??

Hi all, I have a pretty basic question about categorical variables but I can't seem to be able to find answer so I am hoping someone here can help. I found that if the factor names are all in numbers, fitting the model in lm would return labels that are not very recognizable. # Example: let's just assume that we want to fit this model fit <- lm(height ~ age + Seed,

Hallo! I want to plot multiple grouped data in a postscript file using a loop. As I use a loop no plot (or just one empty plot) is generated. Here an example: library(nlme) data(Loblolly) # example data from nlme postscript("PSFile.ps") for (i in 1:1) # just as example { plot(Loblolly) } dev.off() Result: Just an empty PSFile.ps. (Withoput the loop it works. May anybody help? Karl

Hello, I am comparing models made with nlme functions and non-nlme functions, based on Akaike's AIC. The AIC values I get for exactly the same model formulation --for example a linear model with no random effects fit with gls and lm, respectively-- do not fit, although the values of the four model parameters are exactly the same. For example: m1 <- gls(height ~ age, data = Loblolly) m2

I ran into a problem trying to make a plot from a file that's read using source. Basically, I have the following code in a file "plot.R" : library(nlme) data(Loblolly) fm1 <- nlme(height ~ SSasymp(age, Asym, R0, lrc), data = Loblolly, fixed = Asym + R0 + lrc ~ 1, random = Asym ~ 1, start = c(Asym = 103, R0 = -8.5, lrc = -3.3))

R users, I am trying to find some way to find the value of a column that is repeated the most for each StandID of a dataframe. I have research methods online and the help page, but have had no success in finding a solution. I have tried using the table function but it returns items for the whole dataset and not by the StandID. Any help will be appreciated. Thanks in advance. R version 2.11.1

Hey I'm not really sure what I should put on here, but I am having trouble with my R code. I am trying to get the p-values, R^2s etc for a number of different groups of variables that are all in one dataset. This is the code: #Stand counter st<-1 #Collections stands<-numeric(67) slopes<-numeric(67) intercepts<-numeric(67) mses<-numeric(67) rsquares<-numeric(67)

Hola a todos Les consulto por un problema con xtable y nlme, tomando un ejemplo del manual de nlme para obtener los resultados en latex utilizando xtable, se puede utilizar el siguiente código, pero hay un problema y causa error. library(nlme) library(xtable) fm1 <- nlme(height ~ SSasymp(age, Asym, R0, lrc), data = Loblolly, fixed = Asym + R0 + lrc ~ 1,

R-users, Attached is the file (SR_use2.txt) I'd like to include and includes column headers. nat_est is the response variable and is the number of species at a particular point. The other variables are the explanatory vars (vark, var2, var1, UK, U2, U1, GK, G2, G1, PK, P2, P1). Here is Walsh's sample code for hier.part: data(urbanwq) env <- urbanwq[,2,8] hier.part(urbanwq$lec,

can you recommend a good manual for R that starts with a data set and gives demonstrations on what can be done using R? I downloadedR Langauage definition and An introduction to R but haven't found them overly useful. I'd really like to be able to follow some tutorials using a dataset or many datasets. The datasets I have available on R are Data sets in package 'datasets':

Hey, i have a question about how to reorganize a data frame in the easiest way. my example: what would be the easiest way to bring a data.frame such like this: ---- nr height age Seed 1 1 4.51 3 301 2 15 10.89 5 301 3 29 28.72 10 301 4 43 41.74 15 301 5 57 52.70 20 301 6 71 60.92 25 301 7 2 4.55 3 303 8 16 10.92 5 303 9 30 29.07 10 303 10 44

Hello, How do I specify the formula and random effects without a startup object ? I thought it would be a mixture of nls and lme. after trying very hard, I ask for help on using nlme. Can someone hint me to some examples? I constructed a try using the example from nls: #variables are density, conc and Run #all works fine with nls DNase1 <- subset(DNase, Run == 1 ) fm2DNase1 <- nls(

I have a need to call and pass arguments to nlme() from within another function. I use R version 1.8. I have found an apparent way to make this work, but I would appreciate some comments on whether this fix is really appropriate, or there is another way to do it that does not involve changing the source code. I don't have enough experience to start changing the sorurce code of a library

Hello all. I have been working on a (fairly simple) function for way too long. I’ve really hit a wall and I was hoping someone might be able to point me in the right direction. I have (attempted) to create a function that has an embedded while loop. The while loop works fine by itself, however, when the while loop is embedded in the function, the function fails. I’m not sure why this

Hello all, I have created a barplot that shows change in hardwood/softwood density from 1965 to 2005 in 5 year periods (1965,1970, etc). I would like to have an X-axis where the labels for the years line up after every two bars in the plot (there is one bar for hardwood, and another for softwood). Below is my script: density<-read.table("F:\\Megan\\Vtest.csv", header=TRUE,

Hi to all on the Rsync mailing list. I noticed some posts back on the 25th September about my Rsync page (http://www.gaztronics.net/rsync.php). There were questions regarding it's suitability, why it was there, etc, and how I did not conform to the GPL. I have since sought advice from the Rsync administrators, who have advised me with regards to the GPL, etc. The page has been updated

Dear all I try to refine my nlme models and with partial success. The model is refined and fitted (using Pinheiro/Bates book as a tutorial) but when I try to plot plot(augPred(fit4)) I obtain Error in predict.nlme(object, value[1:(nrow(value)/nL), , drop = FALSE], : Levels (0,3.5],(3.5,5],(5,7],(7,Inf] not allowed for vykon.fac > Is it due to the fact that I have unbalanced

I'm trying to use tapply to output means and SD or SE for my data but seem to be limited by how many times I can subset it. Here's a snippet of my data > stems353[1:10,] Time DataSource Plot Elevation Aspect Slope Type Species SizeClass Stems 1 Modern Cameron 70F221 1730 ESE 20 Conifer ABCO Class1 3 2 Modern Cameron 70F221 1730

Ahhh yes silly me. Toying around with it I came to that conclusion. Specifying the secrets file as the source of the password was my problem. I created a file on the local machine and put the password in there just as plain as could be. I made it readable only for root. Then I ran the command and it worked like a charm. Now I am off to figure out a way to make each module log to a different

Good day to all. I would like to congradulate you all on a brilliant product. I have been watching this list for some time as I find the topics most informative. I have a carefully crafted question that I would like to ask you guys. But first some info to help give you a better picture. I am using a RH 9.0 box and the latest rsync to backup different servers to an offsite data facility. The

I am moving from an anonymous rsync setup to a password protected rsync setup. I have my module involved in the test configured to auth only for one user. The secrets file is working just fine for I have tested it. Before I decided to go with a password protected setup, I had a script that cron ran every night at 2 am. worked like a charm. I added the --password-file=/etc/rsyncd.secrets argument