Re: [R] AIC function and Step function

Hi Dana, Many thanks to Christos Hatzis who sent me an offline response, pointing out the new functions that make this much easier than my last suggestions: methods() and getAnywhere() > methods("extractAIC") [1] extractAIC.aov* extractAIC.coxph* extractAIC.glm* extractAIC.lm* ext

Re: [R] AIC function and Step function

Hi Dana, Of course the only true way to know what the AIC calculations are is to read the source code. From within R, what with namespaces, that is becoming increasingly difficult. The AIC() function is not too hard to find from R. > AIC function (object, ..., k = 2) UseMethod("AIC") So now

[R] AIC function and Step function

I would like to figure out the equations for calculating "AIC" in both "step() function" and "AIC () function". They are different. Then I just type "step" in the R console, and found the "AIC" used in "step() function" is "extractAIC". I went to the R help, and found: "The criterion used is

Re: [R] Troubles with the format of dates

Perhaps just this: Lines <- "DATUM P1 P2 2006-11-16 425.21 423.99 2006-12-15 425.12 423.97 2007-01-16 425.16 424.06" library(chron) DF <- read.table(textConnection(Lines), header = TRUE, as.is = TRUE) DF$DATUM <- chron(DF$DATUM, format = "y-m-d", out.form

Re: [R] Package Iso

I have just uploaded a revised version of the Iso package to ``incoming'' on CRAN. It should be available soon. I have fixed the bug that was pointed out to me by the two addressees. I have also added a new feature to the package; it returns a step function representation of the isotoni

Re: [R] "xreg" in ARIMA modelling.

The help file states: "The exact likelihood is computed via a state-space representation of the ARIMA process, and the innovations and their variance found by a Kalman filter." It is possible to include exogenous variables (xreg) this way, but one can only assume this is done [only one person kn

[R] s4 object

Hola muy buen dia!! Me podrían explicar de una manera sencilla en que consiste y como funciona la clase s4 object. Se les agredeceria mucho su mas pronta respuesta! [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

Re: [R] if then statement problem

If I understood your problem correctly, you are just missing a couple of things: q = which(apply(p.unique,2,function(x)all(x==r)) == TRUE) Also, you should probably change this line: if(q>0){c=p.unique[,q]};{c=c(0,0,0)} To something like this: if(length(q)>0){c=p.unique[,q]};{c=c(0,0,0)} Rega

Re: [R] Regression Problem for loop

Dear ales, Try this: # Data # set.seed(123) mat=matrix(rnorm(36),ncol=6) colnames(mat)=paste('x',1:6,sep="") res=rnorm(6) # All info for coefficients. First 4 columns # correspond to the intercept and the next for to the slope t(apply(mat,2,function(x){ sm=summary(lm(res~x))$coeff res=matrix(sm,n

Re: [R] Simultaneous Comparison of Factor Structures

Dear Mark, Although I would like, at some point, to add this feature, the sem package doesn't fit multiple-group models. Sorry, John -- John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario, Canada web: socserv.mcmaster.ca/jfox > -Or

[R] Help with "replicate"

Hello everyone, I have a file called "trans" which has 100 records (rows) and 60 columns. I need to make 2 samples (with 5 records in each) out of this file with no replacement, so I used "replicate". The sampling part works fine and I called the outcome "result". my problem is then I need to ca

[R] if then statement problem

I really need help with an if then statement I am having trouble with. A brief explanation of what I am trying to do: I have a p matrix of all permutations of a vector of length 3 filled with either 1s or -1s (these make up the columns). The p.unique matrix is the unique columns of the p matri

[R] How to optimize a trading singal generator (for-loop)?

Hi guys, I've got some performance troubles with my trading signal generator, which indicates when the system goes long or short. I'm playing with some historical data and the for-loop isn't doing his job very efficient. Is there some vectorial solution for this? Here the for-loop: > trade

Re: [R] 1-Pearson's R Distance

At 5:00 PM +0100 11/27/08, Claudia Beleites wrote: Hi Rodrigo, afaik, (1 - r_Pearson)/2 is used rather than 1 - r_Pearson. This gives a distance measure ranging between 0 and 1 rather than 0 and 2. But after all, dies does not change anything substantial. see e.g. Theodoridis & Koutroumbas: Patt

Re: [R] R and SPSS

Thanks all! Tobias Verbeke wrote: > >>There exists a R plug-in for SPSS. You can find it on the SPSS website. > > ... and there is a page on the R wiki: > > http://wiki.r-project.org/rwiki/doku.php?id=tips:callingr:spss > > HTH, > Tobias > >>> I have a code in R. Could anyone give m

Re: [R] 1-Pearson's R Distance

Rodrigo Aluizio wrote: Hi again List, Well this time I?m writing for a friend (really J). He needs to create a distance matrix based on an abundance matrix using the 1-Pearson?s R index. Well I told him to look at the proxy package, but there is only Pearson Index. He needs it to perform a clust

Re: [R] as.numeric in data.frame, but only where it is possible

On Thu, Nov 27, 2008 at 12:53 AM, Kinoko <[EMAIL PROTECTED]> wrote: > Hi, > > I would like to convert my "character" sequences in my matrix/ > data.frame into "numeric" where it is possible. > I would also like to retain my alphabetic character strings in their > original forms. > "5.1" > 5.1

Re: [R] 1-Pearson's R Distance

Hi Rodrigo, afaik, (1 - r_Pearson)/2 is used rather than 1 - r_Pearson. This gives a distance measure ranging between 0 and 1 rather than 0 and 2. But after all, dies does not change anything substantial. see e.g. Theodoridis & Koutroumbas: Pattern Recognition. I didn't know of the proxy pack

Re: [R] Troubles with the format of dates

Thank you _SO_ much, Petr, you helped me again. That worked like a charm, exactly what I needed! Cheers, Kathi On Thu, 27 Nov 2008 15:01:29 +0100, Petr PIKAL wrote > Hi > > [EMAIL PROTECTED] napsal dne 27.11.2008 14:20:17: > > > Dear useRs, > > > > I'm struggling again with date-related stu

Re: [R] as.numeric in data.frame, but only where it is possible

Hello, On 26 Nov 2008 at 22:53, Kinoko wrote: > Hi, > > I would like to convert my "character" sequences in my matrix/ > data.frame into "numeric" where it is possible. > I would also like to retain my alphabetic character strings in their > original forms. "5.1" > 5.1 "hm" > "hm" >

[R] Simultaneous Comparison of Factor Structures

Dear R Users, I have a situation where I want to compare the factor structures of two samples that were administered then same test. AMOS as well as SAS provide the option of “simultaneous comparison of factor structures” like it was proposed by Jöreskog (1971, Simultaneous factor analysis in

Re: [R] Finding Stopping time

Hi Debanjan, It would be more likely that you get a response if your question was more clear. Your code is very difficult to read, and it doesn't help that you don't provide any context, or comment your code with "### This is calculating the average" kind of statements. What are you trying to do?

Re: [R] Error in Comprting Risks Regression

Hi, That will be difficult to help with the little information you gave. Please read the posting guide and what's at the bottom of this email. Best regards, Arthur Allignol kende jan wrote: Dear All, I am trying to run the following function (a CRR=Competing Risks Regressionmodel) and rece

Re: [R] Reshape with var as fun.aggregate

Thanks, that works fine. ___ Regards, Jason Locklin, B.Sc. Hons. University of Waterloo http://artsweb.uwaterloo.ca/~jalockli On Wednesday 26 November 2008 19:49:48 you wrote: > Hi Jason, > > The reason you are getting this error is that the latest version of > reshape uses fun.aggregate(nume

Re: [R] par(ask=TRUE) and devAskNewPage(ask=TRUE) not working

Have you looked at the TeachingDemos package? -- David Winsemius On Nov 27, 2008, at 7:18 AM, mentor_ wrote: Hi, First, I do some calculation, then a plot. Add some lines and texts to the plot. Second, do some further calculation, then a plot. Add again some lines and texts to the pl

Re: [R] Regression Problem for loop

On Nov 27, 2008, at 9:49 AM, David Winsemius wrote: Have you looked at results of str on a regression object? I would not think that there would be a single p.value associated with such a beast, but that there might be if you examined individual coefficients. ? coefficients ?coef That w

Re: [R] Regression Problem for loop

Have you looked at results of str on a regression object? I would not think that there would be a single p.value associated with such a beast, but that there might be if you examined individual coefficients. ? coefficients ?coef -- David Winsemius On Nov 27, 2008, at 4:03 AM, ales grill wro

[R] 1-Pearson's R Distance

Hi again List, Well this time I’m writing for a friend (really J). He needs to create a distance matrix based on an abundance matrix using the 1-Pearson’s R index. Well I told him to look at the proxy package, but there is only Pearson Index. He needs it to perform a clustering. Well, as soon as h

Re: [R] How to create a string containing '\/' to be used with SED?

What do you mean by "doesn't work" and "the \\" are removed? > a <- "..\\/path\\/file" > a [1] "..\\/path\\/file" > cat(a) ..\/path\/file -- David Winsemius On Nov 26, 2008, at 7:46 PM, ikarus wrote: Thanks Sean!! I followed you suggestion to use gsub() and it worked perfectly! I still

[R] Odp: Troubles with the format of dates

Hi [EMAIL PROTECTED] napsal dne 27.11.2008 14:20:17: > Dear useRs, > > I'm struggling again with date-related stuff: I am using R to draw water > levels at certain measuring > stations. My data comes as a tab-delimited text file and looks like this: > > DATUM P1 P2 > ... > 2006-11-16

Re: [R] lmer refuses nested random factors

Dear Martin, How many levels has PLANT? Does it is an ID for the individual plants? In that case it is pointless to add it to a random effect. Because then each group at the higest level would contain only one observation. There is a mailing list dedicated to mixed effects models: R-Sig-mixed-m

[R] Troubles with the format of dates

Dear useRs, I'm struggling again with date-related stuff: I am using R to draw water levels at certain measuring stations. My data comes as a tab-delimited text file and looks like this: DATUM P1 P2 ... 2006-11-16 425.21 423.99 2006-12-15 425.12 423.97

[R] Error in Comprting Risks Regression

Dear All, I am trying to run the following function (a CRR=Competing Risks Regressionmodel) and receive the error in solve.default. Can anyone give me some insights into where the problem is? Thanks > print(z<-crr(J3500,CD3500,cov)) Error in solve.default(v[[1]]) : Lapack routine dgesv

[R] par(ask=TRUE) and devAskNewPage(ask=TRUE) not working

Hi, First, I do some calculation, then a plot. Add some lines and texts to the plot. Second, do some further calculation, then a plot. Add again some lines and texts to the plot. Third, do some further calculation, then Fourth, . After a plot is complete (means the plot itself, lines and

[R] "xreg" in ARIMA modelling.

Hello, Does anyone know how the parameter estimates are calculated for xreg variables when called as part of an arima() command, or know of any literature that provides this info? In particular, I was wondering if there is a quick way to compare different combinations of "xreg" variables in the

[R] Odp: what is there in a numeric (0)?

Hi [EMAIL PROTECTED] napsal dne 27.11.2008 11:48:32: > > Hello! > I will repeat my question :what is there in a numeric(0)? > > In a function, I have several matrices with dimensions 288x1, also matrices > with dimensions > 0x1 (numeric(0)). > But I really do not know in which matrix th

[R] lmer refuses nested random factors

I am trying to run the following model in R > lmer(leaves.eaten~Geocytotype+(1|TEST/ PLANT),data=cyphoplantfeeding,family=poisson) My experimental setup is 41 replicates (TEST) of an experiment in which there are three Geocytotypes of a plant species in each TEST, and two plant pseudoreplic

[R] Odp: dataframe

Hi [EMAIL PROTECTED] napsal dne 25.11.2008 05:30:52: > > hi there > > I have a dataframe > > abc 123 345 > abc 345 456 > lmn 567 345 > hkl 568 535 > lmn 096 456 > lmn 768 094 > > i want the uniques of column 1 and there corresponsing column 2 and 3 > output > abc 123 345 > lmn 567 345 > h

[R] Help processing large data

Dear all, I have one problem to handle a large dataset... It looks like: "read" "no" "length" 2 2 144 7 7 47490 9 9 310944 11 11 10089 14 14 13152 17 17 27363 and so on There are 13 rows >From this table I need to make a table like 2_1 2 100 2_2 2 44 7_1 7 100 7_2 7 100 ... ... 7_474 7 100

Re: [R] R and SPSS

>There exists a R plug-in for SPSS. You can find it on the SPSS website. ... and there is a page on the R wiki: http://wiki.r-project.org/rwiki/doku.php?id=tips:callingr:spss HTH, Tobias >> >>> I have a code in R. Could anyone give me the best possible way (or just >>> ways!) to integrate it

[R] what is there in a numeric (0)?

Hello! I will repeat my question :what is there in a numeric(0)? In a function, I have several matrices with dimensions 288x1, also matrices with dimensions 0x1 (numeric(0)). But I really do not know in which matrix there are elements or not. So, these matrices that do not contain anythi

Re: [R] survreg and pweibull solved for any distribution

Dear all - Following up my own post, having found a Terry Therneau post about the value of predict(model,type="quantile"), http://tolstoy.newcastle.edu.au/R/e4/help/08/03/5335.html the following code ammendment produces what I was intending to see. It clear that the 0.5 quantile is the in

[R] Regression Problem for loop

Dear all, I have wrote a code for a linear regression. I want to write a loop for so, that I can get estimate for pavlues for six predictors. But I am getting for estmate for only last one. How can I get pvalues for all my predictors in a loop?? Anticipating your help Thanks Ales

Re: [R] Welcome to the "R-help" mailing list (Digest mode)

With igraph, do read.graph(filename, format="ncol", directed=TRUE) I just noticed that the 'directed' option is not documented, sorry about that, Gabor ps. please give a proper subject to your emails... On Thu, Nov 27, 2008 at 9:54 AM, Weijia You <[EMAIL PROTECTED]> wrote: > Hi friends, > > Is

Re: [R] Welcome to the "R-help" mailing list (Digest mode)

Hi friends, Is there anyone who happened to import data set in the "DL" format into R for further analysis? In the package of "network", there's only a method named "read.paj()". So now, I have to get the dl file from the original data set, and use UCINET to convert it to .net file. It's too comp

Re: [R] Error in sqlCopy in RODBC

BKMooney wrote: > > I tried your suggestion... > > library(RODBC) > channel = odbcConnectAccess("db.mdb") > sqlCopy(channel,"Select * from tab","newtab",destchannel=channel, > safer=TRUE,append=TRUE,rownames=FALSE,fast=FALSE) > odbcClose(channel) > > however, I am still running into errors,

Re: [R] R and SPSS

Hi, There exists a R plug-in for SPSS. You can find it on the SPSS website. Hope it helps. Alain Liviu Andronic wrote: > Hello, > > On Wed, Nov 26, 2008 at 9:25 PM, Applejus <[EMAIL PROTECTED]> wrote: > >> I have a code in R. Could anyone give me the best possible way (or just >> ways!) to