Re: [R] prblem with NA

Hi r-help-boun...@r-project.org napsal dne 06.01.2009 02:29:18: Hello - kayj wrote: Hi all I have a data set with the total number of columns =ncol, and the total number of rows=nrow. I am trying to loop over the values and id the value is less than or equal to 100 to be

Re: [R] prblem with NA

Wacek Kusnierczyk wrote: Prof Brian Ripley wrote: On Tue, 6 Jan 2009, Petr PIKAL wrote: [...] Logical vectors can be treated like numeric with TRUE=1 and FALSE=0, More accurately, 'like integer' ('numeric' often means 'double'). why would this be more accurate?

[R] hi I have a problem with implementing R-project

Hi, This is Kiran I tried to work on r-project through C#.Net While calling the r-project I got an error that *System.Runtime.InteropServices.COMException(0x80040013):Exception from HRESULT:0x80040013 at STATCONNECTORSRVLib.StatConnectorClass.Init(String bstrConnectorName) at

[R] hi i have a problem with implementing R-Project in C#

Hi, This is Kiran I tried to work on r-project through C#.Net While calling the r-project I got an error that *System.Runtime.InteropServices.COMException(0x80040013):Exception from HRESULT:0x80040013 at STATCONNECTORSRVLib.StatConnectorClass.Init(String bstrConnectorName) at

[R] Changing Matrix Header

Dear all, I have the following matrix. dat A A A A A A A A A A [1,] 0 0 0 0 0 0 0 0 0 0 [2,] 0 0 0 0 0 0 0 0 0 1 [3,] 0 0 0 0 0 0 0 0 0 2 How can I change it into: [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 0 0 0 0 0

Re: [R] Changing Matrix Header

Hi Gundala, try chopping off the top row like newx-as.matrix(x[2:dim(x)[1],]) OR try changing it to a data frame... new x-data.frame(x,row.names=NULL) #pretty sure its not row.names but there is probably an equivalent for col.names OR look into ?read.table and specify header = F Cheers,

Re: [R] Changing Matrix Header

Hello, colnames( dat ) - NULL will do the trick. Carlos J. Gil Bellosta http://www.datanalytics.com On Tue, 2009-01-06 at 17:14 +0900, Gundala Viswanath wrote: Dear all, I have the following matrix. dat A A A A A A A A A A [1,] 0 0 0 0 0 0 0 0 0 0 [2,] 0 0 0 0 0 0 0

[R] problems with ace (ape-package)

Hi everyone I am trying to reconstruct ancestral character states with R. I used the ace function of the package ape. I have discrete characters and tried to use the ER and ARD model. The problem is, that I get negative lik$anc values. Does anybody by chance have an idea what might be wrong?

Re: [R] C dll compilation + S Poetry example

On 05/01/2009 11:26 PM, Pedro Mardones wrote: Dear all; Working with the following code extracted from the document S Poetry by Patrick Burns (from CRAN), I haven't been able to load the resulting dll into R. The code is basically the calculation of the quadratic form x'Qx: static double

Re: [R] Converting data frame to symmetric matrix

Try this: reshape(d, timevar = x, idvar = y, direction = wide) On Tue, Jan 6, 2009 at 6:29 AM, poas...@umich.edu wrote: Dear Sir or Madam, I have the following data frame (which is just a toy example of my larger dataset) L3 - LETTERS[1:3] x=c(1,1,2,2,3,3,4,4,5,5) y=1:10 d -

Re: [R] RODBC connection die - using more than 1 Rgui/Rcmdr

On 05/01/2009 9:26 PM, jaey.ahn wrote: Yes, I am aware of that. RODBC connection is automatically shutting down without any warnings or error message. (That's why I jotted down Oracle error message when trying to reconnect). I don't think it's a matter of Oracle since I've tested many query

[R] Trouble using ohlcPlot

Hi, I'm trying to create a time series that will work with ohlcPlot: mts-ts(data=c(results$OpenPrice,results$HighPrice,results$LowPrice, results$ClosePrice),c=mts,names=c(Open, High, Low,Close)) ohlcPlot(mts) fails with: Error in if ((!is.mts(x)) || (colnames(x)[1] != Open) || (colnames(x)[2]

Re: [R] prblem with NA

Hi r-help-boun...@r-project.org napsal dne 06.01.2009 10:02:12: On Tue, 6 Jan 2009, Petr PIKAL wrote: [...] Logical vectors can be treated like numeric with TRUE=1 and FALSE=0, More accurately, 'like integer' ('numeric' often means 'double'). I learned that somewhere but it is not

[R] Warning message:In pt(q, df, lower.tail, log.p) : NaNs produced

Hi friends, Any idea why do i get this warning?And also why all computed p-values are NaN. Have shown below what i did in Windows r-console.:-- df c1 c2 1 1 50 2 NA NA 3 4 NA 4 7 6 5 NA 7 6 10 10 r-cor(x=df,y=NULL,use=complete.obs,method=c(pearson)) r c1 c2 c1

[R] [R-pkgs] pmml 1.2.0 (predictive modelling markup language)

Version 1.2.0 of pmml has been released and is available from CRAN. The pmml package (http://rattle.togaware.com/pmml.html) is part of the Rattle data mining suite http://rattle.togaware.com. It generates representations of analytic models built in R using the open standard predictive modelling

[R] Large Dataset

Hi alI, I have a 3.1 GB Dataset ( with 11 coloumns and lots data in int and string). If I use read.table; it takes very long. It seems that my RAM is not big enough (overload) I have 3.2 RAM and 7GB SWAP, 64 Bit Ubuntu. Is there a best sultion to read a large data R? I have seen, that

Re: [R] Large Dataset

type ?memory into R and that will explain what to do... S - Original Message - From: Edwin Sendjaja edw...@web.de To: r-help@r-project.org Sent: Tuesday, January 06, 2009 11:41 AM Subject: [R] Large Dataset Hi alI, I have a 3.1 GB Dataset ( with 11 coloumns and lots data in int

Re: [R] Warning message:In pt(q, df, lower.tail, log.p) : NaNs produced

Obviously your df=0 in pt(STATISTIC, df). Regards, Yihui -- Yihui Xie xieyi...@gmail.com Phone: +86-(0)10-82509086 Fax: +86-(0)10-82509086 Mobile: +86-15810805877 Homepage: http://www.yihui.name School of Statistics, Room 1037, Mingde Main Building, Renmin University of China, Beijing, 100872,

Re: [R] Large Dataset

Hi Simon, Thank for your reply. I have read ?Memory but I dont understand how to use. I am not sure if that can solve my problem. Can you tell me more detail? Thanks, Edwin type ?memory into R and that will explain what to do... S - Original Message - From: Edwin Sendjaja

Re: [R] Warning message:In pt(q, df, lower.tail, log.p) : NaNs produced

Yihui Xie wrote: Obviously your df=0 in pt(STATISTIC, df). However, there are 3 complete observation pairs, not 2, and he wants r as the off-diagonal element and not the whole correlation matrix. (and df denotes two different things, but that is just aesthetics.) Regards, Yihui -- Yihui

Re: [R] Large Dataset

When I do it on a Mac installation I get: Help for the topic memory was not found. Is that a Linux-specific function? Or perhaps you meant to type: ?Memory Which does produce useful information. -- David Winsemius sessionInfo() R version 2.8.0 Patched (2008-11-14 r46932)

Re: [R] Large Dataset

I think he meant: ?Memory edwin When I do it on a Mac installation I get: Help for the topic memory was not found. Is that a Linux-specific function? Or perhaps you meant to type: ?Memory Which does produce useful information. -- David Winsemius sessionInfo() R version 2.8.0

Re: [R] adding a curve with xaxs=i

Unfortunately, I need add=T because there are many other things in the graph. I could swear that it worked 2 years ago. Maybe there's something wrong with an update. Original-Nachricht See if this gets you closer to what you want. curve (hm,lty=1,xaxs=i, xlim=c(0,1),

Re: [R] adding a curve with xaxs=i

Valentina Kraus wrote: Unfortunately, I need add=T because there are many other things in the graph. I could swear that it worked 2 years ago. It can't have worked. Axis settings are fixed before you start adding things. Presumably you used xlim, xaxs et al. on one of the previous commands.

[R] Selecting variables from a list of data.frames

I have a simulation program that generates a data frame for each run. I aggregate the data.frames into a list (df.list). The structure of all data frames is the same, only the values are different. I then want to aggregate the various runs. Currently I use the following method (for three

Re: [R] Large Dataset

Only a couple of weeks ago I had to deal with this. adjust the memory limit as follows, although you might not want 4000, that is quite high memory.limit(size = 4000) Simon. - Original Message - From: Edwin Sendjaja edw...@web.de To: Simon Pickett simon.pick...@bto.org Cc:

Re: [R] Changing Matrix Header

Try this: unname(dat) On Tue, Jan 6, 2009 at 6:14 AM, Gundala Viswanath gunda...@gmail.comwrote: Dear all, I have the following matrix. dat A A A A A A A A A A [1,] 0 0 0 0 0 0 0 0 0 0 [2,] 0 0 0 0 0 0 0 0 0 1 [3,] 0 0 0 0 0 0 0 0 0 2 How can I change it into:

Re: [R] Large Dataset

Yup, it was a typo. But I always try capital if lower case doesnt work, Sorry. - Original Message - From: David Winsemius dwinsem...@comcast.net To: Simon Pickett simon.pick...@bto.org Cc: Edwin Sendjaja edw...@web.de; r-help@r-project.org Sent: Tuesday, January 06, 2009 12:40 PM

Re: [R] Sweave in LATEX

On 1/6/2009 6:44 AM, Mr Derik wrote: Hello I have been setting up my computer to run Sweave. I have got the whole thing working on example files, except that my MikTex returns an Undefined Control Sequence error for \Sexpr and my output file contains verbatim code sequences at the apropriate

Re: [R] Selecting variables from a list of data.frames

Magnus account at zulutime.net writes: I have a simulation program that generates a data frame for each run. I aggregate the data.frames into a list (df.list). The structure of all data frames is the same, only the values are different. I then want to aggregate the various runs.

[R] for loop and if problem

Hi, I'm heaving difficulties with a dataset containing gene names and positions of those genes. Not such a big problem, but each gene has multiple exons so it's hard to say where de gene starts and where it ends. I want the starting and ending position of each gene in my dataset. Attached is the

Re: [R] Error : unused arguments in pairs()

Dear Sarah, Thank you a lot, It does exactly what I need. By the way, I tried doing what Prof. Ripley suggested I just was not able to get it right - I am pretty new to this after all. Thank again, Herwig Sarah Goslee wrote: You didn't do what Prof. Ripley suggested - adding a ...

[R] history: recording

Hi, I'm using windows xp and R 2.8.0 - I wonder what is the command to put in a script that has the same effect as when in a plot you choose menu History and Recording. Best wishes Troels -- Troels Ring - - Department of nephrology - - Aalborg Hospital 9100 Aalborg, Denmark - - +45 99326629

Re: [R] Sweave in LATEX

Thanks for replying. I'm working on a windows XP sp3 PC. I am trying to run the test file provided by R at the moment: % -*- mode: noweb; noweb-default-code-mode: R-mode; -*- \documentclass[a4paper]{article} \title{A Test File}

[R] Semi-definite quadratic programming

Hi, I am currently using solve.QP from the quadprog package to solve some quadratic optimization problems of the form: min[ -d'b + (1/2) b'Db ] under constraints A'b = b_0 solve.QP appears to use an implementation of the Goldfarb and Idnani algorithm. I now have a problem of this form where the

[R] Alternate looping

R-help, I'm using the for control flow to graph plots continuously: par(mfrow=c(3,5), mar=c(0.6,.2,1.2,.2),yaxt=n,xaxt=n) for(j in 1:11) { for(i in 1:15) { species - spAldur[spAldur$ar == 1993+i spAldur$aldur == j,] plot(spec...@coords[,1], spec...@coords[,2], xlim=c(-10,-3.5),

Re: [R] for loop and if problem

I'm heaving difficulties with a dataset containing gene names and positions of those genes. Not such a big problem, but each gene has multiple exons so it's hard to say where de gene starts and where it ends. I want the starting and ending position of each gene in my dataset. Attached is

Re: [R] for loop and if problem

On Tue, Jan 06, 2009 at 07:21:48AM -0800, Sake wrote: I'm heaving difficulties with a dataset containing gene names and positions of those genes. Not such a big problem, but each gene has multiple exons so it's hard to say where de gene starts and where it ends. I want the starting and ending

Re: [R] Sweave in LATEX

Your example works for me. I'd guess there's a problem with the way you handled Sweave.sty, but I don't really know what it would be. Can't you tell MikTex to use the Sweave file from its original location, using -include_directory? Duncan Murdoch On 1/6/2009 10:48 AM, Mr Derik wrote:

Re: [R] Interval censored Data in survreg() with zero values!

Terry Therneau schrieb: Apologies -- you are being more subtle than I thought. Nevertheless, I think that the censoring language isn't quite right. You are thinking of a hierarchical model: z ~ N(Xb, sigma), where Xb is the linear predictor, whatever covariates you think

Re: [R] Large Dataset

Hi, I am not very knowledgeable about this kind of stuff but my guess is that if you have a fairly slow computer and massive data sets there isnt alot you can do except get a better computer, buy more RAM or use something like SAS instead? Hopefully someone else will chip in Edwin, best of

Re: [R] Sweave in LATEX

One other suggestion: use the SweavePDFMikTex function from patchDVI (available on http://r-forge.r-project.org/projects/sweavesearch) to handle both the Sweave call, and the pdflatex call. Duncan Murdoch __ R-help@r-project.org mailing list

Re: [R] Sweave in LATEX

Another thing to try is placing Sweave.sty in the same directory as your .Rnw file and experiment with and without \usepackage{Sweave} in your .Rnw file. On Tue, Jan 6, 2009 at 11:39 AM, Duncan Murdoch murd...@stats.uwo.ca wrote: Your example works for me. I'd guess there's a problem with the

[R] smoothed contour lines

Hi all, I'm trying to draw a contour plot with rounded (smoothed) contour lines instead of sharp angles; when the grid consists of only a few points in x- and y- axis, the resulting contour is in facts rather ugly since very sharp angles may appear. I did not find any way to do it, by using

Re: [R] history: recording

Thanks a lot - windows(record=TRUE) seems to do the job. Sorry not to have found it. Troels jim holtman skrev: ?windows On Tue, Jan 6, 2009 at 10:44 AM, Troels Ring tr...@gvdnet.dk wrote: Hi, I'm using windows xp and R 2.8.0 - I wonder what is the command to put in a script that has the

Re: [R] the first and last observation for each subject

Just in case anyone is still interested, here are some comparisons of the time it says to compute grouped medians via sapply(split(x,group),median) and gm(x,group), which uses the trick used by rle() to find the first and last entries in each group. Which method is fastest depends on the nature

Re: [R] Large Dataset

Hi Ben, Using colClasses doensnt improve the performace much. With the data, I will calculate the mean, min, max, and standard deviance. I have also failed to import the data in a Mysql Database. I dont have much knowledge in Mysql. Edwin Edwin Sendjaja edwin7 at web.de writes: Hi

Re: [R] Large Dataset

The sqldf R package can import a file into an sqlite database and extract a portion of it. You basically need two statements: one to specify the name and format of the file and one to specify what you want to extract. See home page at: http://sqldf.googlecode.com On Tue, Jan 6, 2009 at 12:10 PM,

Re: [R] Large Dataset

For the mean, min, max and standard deviance (deviation I suppose) you don't need to store all data in the memory, you can calculate them incrementally. Read the file line by line (if it is a text file). G. On Tue, Jan 6, 2009 at 6:10 PM, Edwin Sendjaja edw...@web.de wrote: Hi Ben, Using

Re: [R] Large Dataset

Increase the memory as much as you can, read in the data, (however long it takes) then aggregate the data into smaller chunks, selecting only the bits you want. Remove the big original data set from memory (using rm()) and keep (or save the smaller aggregated data using wite.table()) If

[R] Using apply for two datasets

I can run one-sample t-test on an array, for example a matrix myData1, with the following apply(myData1, 2, t.test) Is there a similar fashion using apply() or something else to run 2-sample t-test with datasets from two groups, myData1 and myData2, without looping? TIA, Gang

Re: [R] Using apply for two datasets

I think that you can use mapply for this. On Tue, Jan 6, 2009 at 3:24 PM, Gang Chen gangch...@gmail.com wrote: I can run one-sample t-test on an array, for example a matrix myData1, with the following apply(myData1, 2, t.test) Is there a similar fashion using apply() or something else to

Re: [R] Large Dataset

Bellow, you can see the R data. But this stucks even in first line (read.table..). I dont know how to calculate this and write the result into a new table. Edwin data - read.table(test.data) data - subset(data, (data$Zusatz!=60) (data$Zusatz!=0))

[R] Duplicated messages

Folks: I am getting duplicate messages on posts. Please correct my details if I'm wrong, but I believe it's because folks are posting to **both** the addresses, r-help@r-project.org and r-h...@stat.math.eth.ch . I believe the first is just an alias for the second, and that's why the suplicate

Re: [R] Using apply for two datasets

Hi Gang, Perhaps this post might be useful in this case. Please take a special lookat Gábor Csárdi's reply. http://www.nabble.com/apply,-t-test-and-p-values-to20012292.html#a20012292 HTH, Jorge On Tue, Jan 6, 2009 at 1:10 PM, Gang Chen gangch...@gmail.com wrote: Thanks a lot for the quick

Re: [R] history: recording

On Tue, 6 Jan 2009, Troels Ring wrote: Thanks a lot - windows(record=TRUE) seems to do the job. If you want to control an already open graphics device you will need the R-devel version, which has a function msgWindow() to change a lot of settings, including 'record'. Sorry not to have

[R] Two Noobie questions

1. I have a list of lm (linear model) objects. Is it possible to select, through subscripts, a particular element (say, the intercept) from all the models? I've tried something like this: List[[1:length(list)]][1] All members of the list are similar. My goal is to have a list of the intercepts

Re: [R] Two Noobie questions

On Jan 6, 2009, at 1:50 PM, AllenL wrote: 1. I have a list of lm (linear model) objects. Is it possible to select, through subscripts, a particular element (say, the intercept) from all the models? I've tried something like this: ?coef if your list of models is ml, then perhaps

Re: [R] Alternate looping

On Tue, 06-Jan-2009 at 04:27PM +, Luis Ridao Cruz wrote: | R-help, | | I'm using the for control flow to graph plots continuously: | | par(mfrow=c(3,5), mar=c(0.6,.2,1.2,.2),yaxt=n,xaxt=n) | for(j in 1:11) | { | for(i in 1:15) | { | | species - spAldur[spAldur$ar == 1993+i

[R] sqldf with date (class date) variables ?

Hello, I am having trouble with retrieving some data in queries involve with variables with data type date. See the enclosed example: ll-c(21DEC2006,15DEC2006) ss-data.frame(ll) ss-transform(ss,ll=as.date(as.character(ll))) (ss) ll 1 17156 2 17150 tt-sqldf(select ll from ss) (tt) ll 1

Re: [R] for loop and if problem

On Tue, 6 Jan 2009, Sake wrote: Hi, I'm heaving difficulties with a dataset containing gene names and positions of those genes. Not such a big problem, but each gene has multiple exons so it's hard to say where de gene starts and where it ends. I want the starting and ending position of each

Re: [R] Two Noobie questions

Thanks for your help! I combined the above two to get the following, which seems to work (if somewhat inelegant): int.List-unlist(lapply(lmList, function(x) {coef(x)[1]}),use.names=FALSE) lmList is my list of lm objects. -Allen David Winsemius wrote: On Jan 6, 2009, at 1:50 PM,

Re: [R] Using apply for two datasets

Perhaps you can convert your matrices to data frames as in: mapply(t.test,as.data.frame(myData1),as.data.frame(myData2)) to test by column and mapply(t.test,as.data.frame(t(myData1)),as.data.frame(t(myData2))) to test by row? - Original Message From: Gang Chen gangch...@gmail.com

[R] R SEM package

Does anyone know if the sem package in R can implement a stacked model comparison, for example as in LISREL or AMOS? Thanks, Anthony -- Anthony Steven Dick, Ph.D. Post-Doctoral Fellow Human Neuroscience Laboratory Department of Neurology The University of Chicago 5841 S. Maryland Ave. MC-2030

[R] Overlaying several qqnorm curves in same frame

Hi all, I want to create a rather standard overlaid qqnorm plot on a single variable, with different subgroups of the same dataset plotted using different colors/symbols/etc. (I don't want side-by-side, rather different-colored curves on the same graph) I managed to do it rather tediously using

Re: [R] R SEM package

Dear Anthony, I'm not sure what you mean by a stacked model comparison. If you mean a model fit to multiple groups, then the answer is no. If you mean comparison of nested models, then anova() can be used to get a LR test. See ?anova.sem. I hope this helps, John --

Re: [R] Overlaying several qqnorm curves in same frame

On 1/6/09, Assaf oron assaf.o...@gmail.com wrote: Hi all, I want to create a rather standard overlaid qqnorm plot on a single variable, with different subgroups of the same dataset plotted using different colors/symbols/etc. (I don't want side-by-side, rather different-colored curves on

Re: [R] Using apply for two datasets

Thanks a lot for the suggestions, Jorge and Satoshi Takahama! Both approaches work well... Gang On Tue, Jan 6, 2009 at 2:12 PM, Satoshi Takahama s.takah...@yahoo.com wrote: Perhaps you can convert your matrices to data frames as in: mapply(t.test,as.data.frame(myData1),as.data.frame(myData2))

[R] Contributed Documentation

Hi the list, I wrote a tutorial about S4. Is it possible to have a link to it in the page Contributed Documentation or R Documentation on the CRAN web site ? Who shall I contact ? Christophe __ R-help@r-project.org mailing list

Re: [R] R SEM package

Dear Anthony, As I said, the sem package won't do multiple-group analysis. I would like at some point to add multiple-group analysis to the capabilities of the package but I can't say when I might find time to do that. Sorry, John -- John Fox, Professor Department

[R] Social Networks - how to get the same network layout every time when I plot the network?

Hi, I am using the function plot.network from the network library. However, I noticed that every time when I use the function I get a different layout (the network is rotated?). I would like to get the exact same picture every time when I use the function. Is that possible (maybe set a seed or

[R] Drawing from an empirical distribution

Hi All, Does anybody know if there is a simple way to draw numbers from an empirical distribution? I know that I can plot the empirical cumulative distribution function this easy: plot(ecdf(x)) Now I want to pick a number between 0 and 1 and go back to domain of x. Sounds simple to me. Any

Re: [R] Drawing from an empirical distribution

If the ecdf is 'ecdf(x)', do just: sample(x, size=whatever, replace=TRUE) HTH, Antonio. 2009/1/6 culpritNr1 ig2ar-s...@yahoo.co.uk: Hi All, Does anybody know if there is a simple way to draw numbers from an empirical distribution? I know that I can plot the empirical cumulative

Re: [R] Drawing from an empirical distribution

the empirical distribution gives probability 1/n to each of n observations. rather than sampling the unit interval, just resample the dataset. If x is your dataset, and you want an independent sample of size k, sample(x,size=k,replace=TRUE) albyn On Tue, Jan 06, 2009 at 02:39:17PM

Re: [R] Drawing from an empirical distribution

Thank you. That's exactly what I was looking for. Antonio, Fabio Di Narzo wrote: If the ecdf is 'ecdf(x)', do just: sample(x, size=whatever, replace=TRUE) HTH, Antonio. 2009/1/6 culpritNr1 ig2ar-s...@yahoo.co.uk: Hi All, Does anybody know if there is a simple way to draw numbers

Re: [R] Drawing from an empirical distribution

Sure, but it would be more 'fun' to modify ecdf() slightly to produce an ecqf() function -- essentially reversing the arguments to approxfun()-- and then use ecqf(runif(whatever)) no nit-picking about efficiency, please. url:www.econ.uiuc.edu/~rogerRoger Koenker

Re: [R] Drawing from an empirical distribution

Nit-picking about syntax does seem needed, mea culpa, I intended something more like: Qn - ecqf(x) Qn(runif(whatever)) On Jan 6, 2009, at 5:06 PM, roger koenker wrote: Sure, but it would be more 'fun' to modify ecdf() slightly to produce an ecqf() function -- essentially

[R] Memory Efficiency of Symmetric Matrix

-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 I'm generating a symmetric correlation matrix using a data matrix as input: mat - cor(data.mat) My question is: Is there a more memory efficient way to store this data? For instance, since: all(mat == t(mat)) every value is duplicated, and I should

Re: [R] Social Networks - how to get the same network layout every time when I plot the network?

If you want the same network each time, then set the random seed to the same value: n-6 dat - rbinom(n*(n-1)/2,1,.6) net-diag(n) net[lower.tri(net)] - dat net[upper.tri(net)] - t( net )[upper.tri(net)] net #the network library(network) g-network(net,directed=FALSE) set.seed(1) plot.network(g)

Re: [R] Memory Efficiency of Symmetric Matrix

You can do mat[lower.tri(mat, diag=F)] Søren Fra: r-help-boun...@r-project.org på vegne af Nathan S. Watson-Haigh Sendt: on 07-01-2009 01:28 Til: r-help@r-project.org Emne: [R] Memory Efficiency of Symmetric Matrix -BEGIN PGP SIGNED MESSAGE- Hash:

[R] How-To extract data from a matrix of lists subject to a set of given constraints

My goal is to store the DWT coefficients from a number of time series in such a way to be able to extract subsets satisfying giben conditions. I have defined the following 11 matrices whose x-axis represent time intervals and y-axis represent the number of time series: MS - 11 d1.mat -

Re: [R] Contributed Documentation

Christophe Genolini cgenolin at u-paris10.fr writes: Hi the list, I wrote a tutorial about S4. Is it possible to have a link to it in the page Contributed Documentation or R Documentation on the CRAN web site ? Who shall I contact ? Christophe Try c...@r-project.org . Posting to

Re: [R] Memory Efficiency of Symmetric Matrix

See also the dist() function documentation. If you use indexing as described in ?dist it is straightforward to maintain and use a vector of the distances. Steven McKinney Statistician Molecular Oncology and Breast Cancer Program British Columbia Cancer Research Centre email: smckinney +at+

Re: [R] Generating GUI for r-scripts

Daren Tan daren76 at hotmail.com writes: I have developed some scripts that basically ask for input tab-limited format files, do some processing, and output several pictures or csv. Now I need to have some gui to wrap on top of the scripts, so that end-users can select their input files,

Re: [R] Memory Efficiency of Symmetric Matrix

the SparseM package might be what you are looking for http://www.econ.uiuc.edu/~roger/research/sparse/SparseM.pdf On Jan 7, 11:36 am, Søren Højsgaard soren.hojsga...@agrsci.dk wrote: You can do mat[lower.tri(mat, diag=F)] Søren Fra:

Re: [R] RODBC connection die - using more than 1 Rgui/Rcmdr

First of all, I apologize for misleading you by mentioning irrelevant error message. Forget about the error message I stated in the first place. I should have clearly stated in the first place. No error message was shown at the moment of the problem. Without any warnings and error message shown,

[R] breaking a loop in R

Hi All, I was wondering if there is anything that breaks a loop in R. For example, For (k in 1:100) For ( i in 1:10){ If ( condition ){ Break the inner loop } } } In the above case, if the program runs the if statement, then I want the inner loop for (i in 1:10) to stop looping and skip

Re: [R] breaking a loop in R

something like this should do it: breakFlag - FALSE For (k in 1:100) For ( i in 1:10){ If ( condition ){ breakFlag - TRUE break } } if (breakFlag) break } On Tue, Jan 6, 2009 at 7:58 PM, kayj kjaj...@yahoo.com wrote: Hi All, I was

Re: [R] Memory Efficiency of Symmetric Matrix

-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 andrew wrote: the SparseM package might be what you are looking for http://www.econ.uiuc.edu/~roger/research/sparse/SparseM.pdf On Jan 7, 11:36 am, Søren Højsgaard soren.hojsga...@agrsci.dk wrote: You can do mat[lower.tri(mat, diag=F)] Søren

Re: [R] Converting data frame to symmetric matrix

try this: L3 - LETTERS[1:3] x=c(1,1,2,2,3,3,4,4,5,5) y=1:10 d - data.frame(cbind(x,y), fac=I(sample(L3, 10, replace=TRUE))) m.out - matrix(ncol=10, nrow=10) m.out[cbind(d$x, d$y)] - d$fac On Tue, Jan 6, 2009 at 3:29 AM, poas...@umich.edu wrote: Dear Sir or Madam, I have the following data

Re: [R] breaking a loop in R

? `break` On Tue, Jan 6, 2009 at 7:58 PM, kayj kjaj...@yahoo.com wrote: I was wondering if there is anything that breaks a loop in R [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] Understanding dsyrk_ in C code

-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 I'm trying to understand some C code in an R package I'm using. I'm address this question here as it's matrix algebra...and I'm no pro at that! the C command reads: double alpha = 1.0, beta = 0.0; dsyrk_(L, N, nGenes, nGenes, alpha, mat1, nGenes,

Re: [R] Contributed Documentation

I would too. On Tue, Jan 6, 2009 at 8:17 PM, Ben Bolker bol...@ufl.edu wrote: Christophe Genolini cgenolin at u-paris10.fr writes: Hi the list, I wrote a tutorial about S4. Is it possible to have a link to it in the page Contributed Documentation or R Documentation on the CRAN web site ?

[R] Bar Plot with Connected Points on 1 Y-Axis

Hi Everyone, Have created a bar plot of the data below using the following code: barplot(TACC,space=0,names.arg=Year). I now want to add a series of connected points to represent the catch. I tried to do this using line(Catch) or points(Catch), however both of these commands result in each

[R] Destroying Attributes of a Matrix

Dear all, I have the following structure of a matrix str(mdat) num [1:32268, 1:10] 0 0 0 0 0 0 0 0 0 0 ... - attr(*, dimnames)=List of 2 ..$ : NULL ..$ : NULL How can I destroy the attributes such that it simply gives: str(mdat) num [1:32268, 1:10] 0 0 0 0 0 0 0 0 0 0 ... - Gundala

Re: [R] Bar Plot with Connected Points on 1 Y-Axis

on 01/06/2009 09:07 PM jimdare wrote: Hi Everyone, Have created a bar plot of the data below using the following code: barplot(TACC,space=0,names.arg=Year). I now want to add a series of connected points to represent the catch. I tried to do this using line(Catch) or points(Catch),

Re: [R] Destroying Attributes of a Matrix

In this case, how about? dimnames(mdat) - NULL Gundala Viswanath wrote: Dear all, I have the following structure of a matrix str(mdat) num [1:32268, 1:10] 0 0 0 0 0 0 0 0 0 0 ... - attr(*, dimnames)=List of 2 ..$ : NULL ..$ : NULL How can I destroy the attributes such that it simply

[R] heatmap.2 and three colours for specific r anges‏

Hi, I hope to show a heatmap with thre colours, no gradation. How to specify heatmap.2 to map green for values less than -1, gray for values between -1 and 1, and red for values greater than 1 ? Thanks [[alternative HTML version deleted]] __

[R] directory of files in a zip file?

Hi, All: How can I get the list of files contained in a zip file? zip.file.extract will extract a specific 'file' from 'zipname', but how can I get the names of the files in 'zipname'? Thanks, Spencer __ R-help@r-project.org

[R] R Threatens SAS According to The NYT

On Wed, Jan 7, 2009 at 11:05 AM, ajay ohri ajayo...@yahoo.com wrote: FYI..not a R -Help Topic, buy I dont know which list to post discussions like this. Regards, Ajay -- Forwarded message -- From: ajay ohri ajayo...@yahoo.com Date: Wed, Jan 7, 2009 at 10:46 AM Subject:

  1   2   >