Re: [R] SSPIR problem

2006-06-17 Thread Spencer Graves
Did you try the 'examples' in the SS help file? If no, please copy them from the help file into a script file and work through them line by line, trying different modifications. If you've tried this, have you similarly worked through the 'sspir' 'demo' files? (These can be found in

Re: [R] Question concerning mle

2006-06-17 Thread Spencer Graves
It's difficult to say much at this level of generality, but I have four suggestions: 1. Have you tried creating a reasonable grid of starting values using "expand.grid" and then plotting the resulting likelihood surface? If you have more than 2 parameters, you may want to

Re: [R] problem with legend on other graphics devices

2006-06-17 Thread Gavin Simpson
On Fri, 2006-06-16 at 19:42 +0100, Gavin Simpson wrote: > Hi, > > I have a bit of a problem with a legend in png, eps and pdf plots > produced from a custom plotting function. I was writing a little > function to produce some stock plots for a routine analysis conducted in > our lab. I have a wrap

Re: [R] Having a go at Rpad

2006-06-17 Thread tshort
That means that Apache doesn't know to execute the perl files. Something's wrong with your Apache configuration or directory permissions. Some things to check: * Make sure the *.pl files in ../Rpad/server/ are set to enable execution. * Check your apache http.conf. You need something like the fo

Re: [R] substr or split help needed

2006-06-17 Thread Gabor Grothendieck
Try this: # read in lines and replace all occurences of Simula, all occurrences # of .txt, all commas and all underscores with spaces Lines <- gsub(".*Simula|.txt|[,_]", " ", readLines("/myfile.dat")) # reread what is left skipping over column headings and setting it ourself # If you want factor

[R] Panel Bar plots

2006-06-17 Thread Nealy, David
Hello All, I have been trying to use the par(mfrow(c(m,n)) function to produce a panel of m rows and n columns of bar plots. The different columns were to designate n different categories, as the m rows are to designate different categories as well. However, I cannot seem to get rid of the y-

Re: [R] substr or split help needed

2006-06-17 Thread jim holtman
I used your sample of data. When reading in I used 'as.is=TRUE' to prevent the conversion to 'factor'. I then did an explicit conversion to number on the PLAND column. You should use 'str' to look at the structure of your data frame. You have backsplashes and your 'print' will not show them. Y

Re: [R] lmer binomial model overestimating data?

2006-06-17 Thread Spencer Graves
Martin Henry H. Stevens wrote: > Hi folks, > Warning: I don't know if the result I am getting makes sense, so this > may be a statistics question. > > The fitted values from my binomial lmer mixed model seem to > consistently overestimate the cell means, and I don't know why. I > assum

Re: [R] lmer and mixed effects logistic regression

2006-06-17 Thread Spencer Graves
'lmer' RETURNS AN ERROR WHEN SAS NLMIXED RETURNS AN ANSWER Like you, I would expect lmer to return an answer when SAS NLMIXED does, and I'm concerned that it returns an error message instead. Your example is not self contained, and I've been unable to get the result you got.

Re: [R] managing data

2006-06-17 Thread Neuro LeSuperHéros
Hi Change V1 for the name of the column to be matched on X and Y cbind(X,Y[match(X$V1,Y$V1),]) Neurorox >From: "yohannes alazar" <[EMAIL PROTECTED]> >To: r-help@stat.math.ethz.ch >Subject: [R] managing data >Date: Sat, 17 Jun 2006 14:40:47 +0100 > >Dear mailing list, may some one be kind to he

Re: [R] managing data

2006-06-17 Thread Philipp Pagel
On Sat, Jun 17, 2006 at 02:40:47PM +0100, yohannes alazar wrote: > Dear mailing list, may some one be kind to help me solve following problem. > > I am trying to write a code that will combine two tables "x" and "y". The > first columns of both tables are unique identification for the rows. The >

[R] substr or split help needed

2006-06-17 Thread Milton Cezar
Dear R-friends I have several data files with about 1,900 lines (records) each. I´m using read.table command to read the files. The files looks like LID , TYPE , PLAND D:\Bijou-MC\Simula_P005_H100_R001.txt , Forest , NA

[R] Having a go at Rpad

2006-06-17 Thread vittorio
Context: FreeBSD 6.1on a pentium 4; apache 2 server up & running; mod_perl installed; java & php 4 working. Willing to use Rpad on an office intranet, I installed the Rpad package (R2HTML was still there) and carefully read the instructions to set up a server. Now, on the server I can start t

[R] managing data

2006-06-17 Thread yohannes alazar
Dear mailing list, may some one be kind to help me solve following problem. I am trying to write a code that will combine two tables "x" and "y". The first columns of both tables are unique identification for the rows. The first column of table "X" is a sub set of the first column of "Y". I need t

[R] Projection pursuit regression method

2006-06-17 Thread anil kumar rohilla
Hi, list i was trying to do( Projection pursuit regression ) ppr function in Mass library . It is running very good ,no problem ,but when i refered the references Friedman et al 1981 ,,and SMART guide,, i was not able to find how parametes they r geeting using least square method , if any bo

[R] Getting forcasting equation from nnet results

2006-06-17 Thread juliacel
I'm trying to build forecasting equation from weights of 2-2-1 neural net. Running the nnet function gives me a vecto of 9 weights, but I don't know how to build the equation form these values. Can anyone advice? Or at least tell me where the nnet output is described in details (the manual only gi

Re: [R] modeling logit(y/n) using lrm

2006-06-17 Thread Frank E Harrell Jr
Hamilton, Cody wrote: > After a little digging, it turns out that fit.mult.impute will allow > fitter = glm, so previous suggestions regarding modeling cbind(y,n) as > an outcome will work fine. Thanks! Also lrm can easily handle your setup using the weights argument. Frank > > > > Cody Hami