similar to: read.table and factor

Greetings R Help Group, How does one effect a multiphasic logistic growth model with 4 phases (e.g. Koops 1986; Weigel, Craig, Bidwell and Bates 1992; Grossman and Koops 2003) with R. Before writing to the group, the R help archives were searched, the web was searched with Google, Venables and Ripley 2002 was consulted, Pinheiro and Bates 2000 was consulted, Bates and Watts 2007 was bought and

I have the following problem. Suppose I save the data from an experiment in the following way. There are 21 pairs of std and cf. Each pairing is an experimental condition, and 4 data values are saved. The conditions are run (and data saved) in random order, though blocked such that all conditions are run once, then all in random order a second time, and so on (let's say 3 or 4 repetitions)

I use AfterStep on Linux. I made an R.xpm file (attached) and put this icon in the "wharf". To invoke R I click on the R icon. I did not find an offical R.xpm, and so I am submitting this one for your collective consideration.... The idea behind the icon is that it looks like a scatterplot, which is "statistical". My main beefs with my R.xpm are that it doesn't have a

Greetings, I am trying the get the standard errors of multiple comparisons using Tukey's HSD. These are not reported by the function TukeyHSD. When I apply the following code to the data, which I store as PROLE4.TXT, several unexpected things happen. First, the function TukeyHSD works for all the comparisons but the function simint doesn't. Second, after the application of na.omit

Full_Name: Version: 1.0.0 OS: linux Submission from: (NULL) (193.62.250.209) ?par lists and describes xlog and ylog, though these are obsolete. It does not list log. I found out about log from ?plot.defaults--so that help is OK. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

Hi, I had a look at the MASS4 scripts in the MASS package, in Ch 11.3 Factor Analysis, there is a section of codes like: data(ability.cov) ability.FA <- factanal(covmat = ability.cov, factors = 1) ability.FA (ability.FA <- update(ability.FA, factors = 2)) #summary(ability.FA) round(loadings(ability.FA) %*% t(loadings(ability.FA)) + diag(ability.FA$uniq), 3)

I was wondering if anyone out there has created a series of images (e.g. jpegs) using R, then animated them on a web page. I have looked around a bit and it seems fairly complicated. One approach that seems sensible is to use a java script that displays the list of images. I have only used html on web pages before so I don't know how to use java... Ideally the final result is easily ported to

B??dnie zakodowany tekst zosta? usuni?ty... Plik: nie znany Url: https://stat.ethz.ch/pipermail/r-help/attachments/20060719/3a883e4f/attachment.pl

I have not really used lm before and I was hoping for some help on a simple problem. Here is a toy version of the problem I want to solve. y x grp -.9 1 a -.8 2 a -.7 3 a -.7 1.5 b -.5 2.5 b -.3 3.5 b -.19 2.7 c -.11 3.7 c -.41 4.7 c I want to fit a model that has one y-intercept and three slopes, one for

Full_Name: Bill Simpson Version: 65.1 OS: linux Submission from: (NULL) (193.62.250.209) Try this code as separate lines entered interactively, then try doit() doit<-function() { x<-seq(1,10) y1<-x+rnorm(10,1,1.5) y2<-x+rnorm(10,1,1) t.test(x,y1,paired=TRUE) t.test(x,y2,paired=TRUE) } doit() apparently executes only the last of a series of t.test()s. Maybe this is no bug,

There were some questions about hist() a couple of days ago which triggered this post. My question/suggestion is about the y-axis in hist. There are reasons to prefer making the y-axis density=relative frequency/bin width. One reason is that the height of the plot does not depend on the bin width; another is that if your histogram is in density then you can easily superimpose a smooth theoretical

Dear expe-R-ts, I try some test data for a factorAnalysis (resp. pca) in the sense of Prof. Ripley's MASS ? 11.1, p. 330 ff., just to prepare myself for an analysis of my own empirical data using R (instead of SPSS). 1. the data. ## The test data is (from the book of Backhaus et al.: Multivariate ## Analysemethoden. Springer 2000 [9th ed.], p. 300 ff):

Hi, are there readily available R packages that are able to perform FA on ordinal and/or nominal data? If not, what other approaches and helpful packages would you suggest? BR, Jay

I'm struggling to explain an ASAN report I'm now getting that I didn't get previously on the same code. In fact the report only happens with -O2 and not when I remove the -O flags which makes it hard to debug and makes me suspect it's dependent on exactly which instructions the code generation decides to access the bytes involved. Afaict the C code shouldn't be accessing the

Full_Name: Bill Simpson Version: 0.64.1 OS: linux Submission from: (NULL) (193.62.250.209) Two problems with dev.print help. 1. Under Examples: dev.print(width=6, height=6, horizontal=F) ^ FALSE 2. dev.print help does not describe any of its parameters (e.g. height or width). Not described on the other dev.xxx pages either

On Wednesday, 26 February 2020 10:43:27 CET Richard W.M. Jones wrote: > On Wed, Feb 26, 2020 at 11:21:18AM +0200, Veselin Kozhuharski wrote: > > Hallo Rich, > > > > Here is the fd list and total number just before collectd application > > crashes. Before that the number of used fd's is constantly increasing. It > > looks like a fd leak inside libguestfs to me.

Hi I have used fa() to perform a factor analysis of a psychological battery which is thought to have 11 factors. I can identify which factors the loadings relate to easily enough because I can see which items are loading onto each of the columns in the $loading output. However, how can I identify which items or loadings are being used to create each of the columns in the $scores output? I have

Hi R users: I don?t know if anybody have had the same problem with 'step()' funtion. If I type the commands ------------------------------------------------------------------------------- testData<-read.table("testData.dat",header=T) model1J<-glm(MCHNV~offset(Offset1),data=testData,family="poisson") step(model1J,direction="forward",

I just installed 99.0a, and I can't set mar and mgp: > par(mar(2,2,.1,.1)) Error in par(mar(2, 2, .1, .1)) : couldn't find function "mar" I had the same happen with mgp. Are these bugs or a problem in my installation (everything else seems OK)? Is R core going to kill mar and mgp? If so, how should I fix the huge margins and huge space between axis title and axis? Try this

Full_Name: Bill Simpson Version: 65.1 OS: Linux Redhat 6.1 Submission from: (NULL) (193.62.250.209) Try this: #default plot symbols and lettering are too small, need to scale up par(cex=2,mex=.7) par(mar=c(5,5,1,1)) spdiff<-c(1,2,3,4,5) dpdet<-c(1,2,3,4,5)/10 dpsp<-c(2,3,4,5,6)/10 dpdir<-c(2,4,6,8,10)/10 plot(spdiff,dpdet,pch=15,ylim=c(0,max(dpdet,dpsp,dpdir)),xlab="Speed