[R] rev matrix

Hello I am looking for a function which performs the same thing as rev but on a matrix rather then a vector, based on one column. I'd be happy to read the man page if such a function exists. Thanks, Vera Lindsay __ [EMAIL PROTECTED] mailing list

[R] Gumbel distribution

Does R have built in Gumbel distribution (pdf, ecdf, hazard, parameters estimation) for the minimum case? Thanks Anne [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

Re: [R] rev matrix

On Thu, 14 Oct 2004 [EMAIL PROTECTED] wrote: Hello I am looking for a function which performs the same thing as rev but on a matrix rather then a vector, based on one column. I'd be happy to read the man page if such a function exists. If I understand you, it is easy to do by indexing, but

Re: [R] Gumbel distribution

please, take a look at package evd Matthias Does R have built in Gumbel distribution (pdf, ecdf, hazard, parameters estimation) for the minimum case? Thanks Anne [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list

Re: [R] maximum likelihood estimation

Roger == Roger D Peng [EMAIL PROTECTED] on Thu, 14 Oct 2004 17:06:25 -0400 writes: Roger What lead you to believe that mle() is defunct? It's Roger still in the `stats4' package in my installation of Roger R. yes indeed. Just to clarify possible confusions: The 'mle'

Re: [R] Problem with number characters

Spencer == Spencer Graves [EMAIL PROTECTED] on Thu, 14 Oct 2004 13:41:24 -0700 writes: Spencer It looks like you have several non-printing Spencer characters. nchar will give you the total number Spencer of characters in each character string. Spencer strsplit can break

Re: [R] C/C++

Hi, Unfortunately I have only some beginners hints to this very interesting topic. Writing R extensions provides examples how to execute r functions from within C code. There is also somewhere a doc how to use R as a shared library (?). But I dont remember in which doc exactly. Furthermore

[R] R and Java

Hi Bobby, You are right - functions is a StringBuffer object. I've compiled without the syntax error. The error message is: Error in dev.print( can only read screen device Do you think this means that RTerm doesn't actually produce a Screen plot when

Re: [R] cluster analysis

Dear James, sorry, this is not really an answer. I use cutree to obtain clusters from an hclust object. I do not get from the identify help page that identify should do anything like what you expect it to do... I tried it out and to my surprise it behaved as you said, i.e., it indeed does

Re: [R] cluster analysis

ChrisH == Christian Hennig [EMAIL PROTECTED] on Fri, 15 Oct 2004 11:43:53 +0200 (MEST) writes: ChrisH Dear James, ChrisH sorry, this is not really an answer. nor this. I'm answering Christian... ChrisH I use cutree to obtain clusters from an hclust ChrisH object. I do

Re: [R] cluster analysis

On Friday 15 Oct 2004 10:43 am, you wrote: PS: It seems that each value is typed twice because classi is named, and each value is also a name. Try as.vector(classi). (Perhaps a little useful help in the end?) Indeed. I have tried, for example: as.vector(classi[[1]]) and

[R] Announcing books on R - policy {was Statistical analys...}

Vadim == Vadim Ogranovich [EMAIL PROTECTED] on Thu, 14 Oct 2004 10:50:08 -0700 writes: Vadim I thought that maybe authors of books on R should be Vadim allowed (encouraged ?) to announce Vadim availability/revisions of their books via the Vadim R-packages list? Well, most

[R] Testing for normality of residuals in a regression model

Hi all, Is it possible to have a test value for assessing the normality of residuals from a linear regression model, instead of simply relying on qqplots? I've tried to use fitdistr to try and fit the residuals with a normal distribution, but fitdsitr only returns the parameters of the

Re: [R] cluster analysis

On Friday 15 Oct 2004 11:02 am, you wrote: or unname(classi) -- which is slightly more expressive in this case and possibly more desirable in other situations. Martin Maechler, ETH Zurich Thanks, Martin. I've tried, like you suggested: un_classi - unname(classi) but

[R] Re: Testing for normality of residuals in a regression model

Dear Federico, see: ? shapiro.test(stats) Shapiro-Wilk Normality Test and ? jarque.bera.test(tseries) Jarque-Bera Test They are the most common tests used for normality testing. Ciao Vito Federico Gherardini wrote on Fri Oct 15 14:44:18 CEST 2004: Hi all, Is it

Re: [R] Testing for normality of residuals in a regression model

What about shapiro.test(resid(fit.object)) Stefano On Fri, Oct 15, 2004 at 02:44:18PM +0200, Federico Gherardini wrote: Hi all, Is it possible to have a test value for assessing the normality of residuals from a linear regression model, instead of simply relying on qqplots? I've tried

[R] length with missing values

R-help I have a martix with missing values( in which I want the sample size by column) When I : apply(matrix,2,length) I get the length of the vector regardless of missing values. I can't pass an argument to length in apply. Alternatively I could ifelse ( is.na ( matrix [, columns in matrix

Re: [R] Testing for normality of residuals in a regression model

Hi Frederico, take also a look at the package nortest: help(package=nortest) Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/16/396887 Fax: +32/16/337015 Web:

Re: [R] Filling polygons with points

[EMAIL PROTECTED] wrote : Are there any possibilities to fill a polygon with a point pattern or with a symbol pattern like '+' oder '-' instead of shading lines? In `spatstat' there is code for testing whether points belong to a spatial region. The region can consist of any number of

Re: [R] length with missing values

Hi Luis, just use: apply( mat, 2, function(x) length(x[!is.na(x)]) ) I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/16/396887 Fax: +32/16/337015

RE: [R] length with missing values

Hi Luis, Are you looking for the number of missings in each column? If so, try this:- bob-matrix(c(2:4, NA, 5:7, NA), ncol=2) colSums(is.na(bob)) [1] 1 1 Hope this helps, Ian. mangosolutions R Consulting and Training Services Tel +44 118 902 6620 Fax +44 118 902 6401 -Original

Re: [R] length with missing values

Dimitris Rizopoulos [EMAIL PROTECTED] writes: Hi Luis, just use: apply( mat, 2, function(x) length(x[!is.na(x)]) ) or: apply(!is.na(mat),2,sum) -- O__ Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of Biostatistics 2200 Cph. N (*) \(*) -- University of

Re: [R] length with missing values

or even: colSums(!is.na(mat)) Best, Dimitris - Original Message - From: Peter Dalgaard [EMAIL PROTECTED] To: Dimitris Rizopoulos [EMAIL PROTECTED] Cc: Luis Rideau Cruz [EMAIL PROTECTED]; [EMAIL PROTECTED] Sent: Friday, October 15, 2004 1:43 PM Subject: Re: [R] length with missing values

Re: [R] length with missing values

On 15 Oct 2004 at 12:01, Luis Rideau Cruz wrote: R-help I have a martix with missing values( in which I want the sample size by column) When I : apply(matrix,2,length) Hi Try apply(matrix, 2, function(x) sum(!is.na(x))) matrix is reserved word, so it is preferable to use other name

Re: [R] length with missing values

On Fri, 15 Oct 2004, Sundar Dorai-Raj wrote: Luis Rideau Cruz wrote: I have a martix with missing values( in which I want the sample size by column) When I : apply(matrix,2,length) I get the length of the vector regardless of missing values. I can't pass an argument to length

Re: [R] R plot problems

At 10:51 2004-10-15 +0100, you wrote: If you are using base grafics, you use the argument srt (see help(par)), How do you get a, say, 45 degree rotation of the axis labels using `srt'? The help page for par's `las' argument says Note that other string/character rotation (via argument srt to par)

RE: [R] Testing for normality of residuals in a regression model

Dear Federico, A problem with applying a standard test of normality to LS residuals is that the residuals are correlated and heterskedastic even if the standard assumptions of the model hold. In a large sample, this is unlikely to be problematic (unless there's an unusual data configuration), but

Re: [R] R plot problems

On Fri, 2004-10-15 at 07:18, Henric Nilsson wrote: At 10:51 2004-10-15 +0100, you wrote: If you are using base grafics, you use the argument srt (see help(par)), How do you get a, say, 45 degree rotation of the axis labels using `srt'? The help page for par's `las' argument says Note

Re: [R] cluster analysis

James == James Foadi [EMAIL PROTECTED] on Fri, 15 Oct 2004 11:36:14 +0100 writes: James On Friday 15 Oct 2004 11:02 am, you wrote: or unname(classi) -- which is slightly more expressive in this case and possibly more desirable in other situations. Martin

RE: [R] R plot problems

Use text(), for which srt will work. As an example: plot(rnorm(20), rnorm(20), xaxt=n) text(-2:2, rep(par(usr)[3], 5), c(negative two, negative one, zero, one, two), srt=45, xpd=NA, adj=1) HTH, Andy From: Henric Nilsson At 10:51 2004-10-15 +0100, you wrote: If you are using

Re: [R] C/C++

Have you checked Section 5.14 (page 64) of R-extension titled Using these functions in your own C code. Section 5.7 tells you about some of the available subroutines. On Fri, 2004-10-15 at 01:31, doktora v wrote: Hey everyone, I have been looking for a while for ways to integrate R's

[R] edit plots from the ADE4 package

Dear R-Help i have a newbie question, how edit size fonts and colors in the plots (scatter.acm or dudi.acm function) of multiple correspondence analysis in the ADE4 Package? thanks? Rafael Gutierrez Estadístico Unidad de Tecnología Cerro Matoso S.A. Tel. 4-7723350 Fax. 4-7723236

[R] constrained splines in GLM

Hi, I would like to use constrained splines in a GLM model (Poisson link) to take into account the nonlinear effect of some covariate. The constraints I need are described below. I have several variables that I need concurrently in the same model. I looked at package mgcv but I do not know if/how

Re: [R] Testing for normality of residuals in a regression model

Thank you very much for your suggestions! The residuals come from a gls model, because I had to correct for heteroscedasticity using a weighted regression... can I simply apply one of these tests (like shapiro.test) to the standardized residuals from my gls model? Cheers, Federico

[R] pdf device --- axes labels text disappeared?

Dear R Wizards: Running R 1.9.1. on amd64. Promise- c(0,20,40); Expect- c(0, 20, 0.2*20+.8*40 ); # this omits printing numbers on the axes labels. pdf(file = bug.pdf ); plot(Promise, Expect, type=b, ylim=c(0,60)); dev.off(); # this works postscript(file = bug.eps ); plot(Promise, Expect, type=b,

Re: [R] Testing for normality of residuals in a regression model

John Fox wrote: Dear Federico, A problem with applying a standard test of normality to LS residuals is that the residuals are correlated and heterskedastic even if the standard assumptions of the model hold. In a large sample, this is unlikely to be problematic (unless there's an unusual data

Re: [R] Testing for normality of residuals in a regression model

John Fox wrote: Dear Federico, A problem with applying a standard test of normality to LS residuals is that the residuals are correlated and heterskedastic even if the standard assumptions of the model hold. In a large sample, this is unlikely to be problematic (unless there's an unusual data

[R] Problems saving a graphic

Hi! I'm with a little problem when saving a graphic into pdf. The graphic has the label: Milhões de toneladas. The problem is with the accent. After I save the graphic into pdf the õ doesn't appear. Somebody could help me? Thanks, Talita Perciano Costa Leite Graduanda em Ciência da Computação

RE: [R] edit plots from the ADE4 package

Hi Try to change the parameters in the s.class function corresponding to col (which is 1 by default) and using the parameter font as a small positive integer determining a text font for characters and hence an interline spacing. Leonardo Trujillo Southamtpton Statistical Sciences Research

[R] does rpart has similar function as edit.tree?

Hi, list, The 'tree' package has edit.tree to change split variables by user. Does 'rpart' package have something similar? I couldn't find either in the manual or archive. Thanks in advance. Best, Auston __ [EMAIL PROTECTED] mailing list

RE: [R] pdf device --- axes labels text disappeared?

Dear R Wizards: Running R 1.9.1. on amd64. Promise- c(0,20,40); Expect- c(0, 20, 0.2*20+.8*40 ); # this omits printing numbers on the axes labels. pdf(file = bug.pdf ); plot(Promise, Expect, type=b, ylim=c(0,60)); dev.off(); # this works postscript(file = bug.eps ); plot(Promise,

Re: [R] pdf device --- axes labels text disappeared?

On Fri, 2004-10-15 at 09:13, [EMAIL PROTECTED] wrote: Dear R Wizards: Running R 1.9.1. on amd64. Promise- c(0,20,40); Expect- c(0, 20, 0.2*20+.8*40 ); # this omits printing numbers on the axes labels. pdf(file = bug.pdf ); plot(Promise, Expect, type=b, ylim=c(0,60)); dev.off(); #

Re: [R] Problem with number characters

Note that there are also regexp classes that define certain character sets, most notably [:graph:] , which can make it easy to create appropriate regexps. More is in ?regex . Martin Maechler maechler at stat.math.ethz.ch writes: : : Spencer == Spencer Graves spencer.graves at pdf.com :

[R] Within-Subjects ANOVA TukeyHSD

Is there a problem with using TukeyHSD for individual comparisons of means coming from a within-subjects ANOVA? I always get: no applicable method for TukeyHSD Thank you, Steffen __ [EMAIL PROTECTED] mailing list

Re: [R] pdf device --- axes labels text disappeared?

No problems under R 2.0.0 or 1.9.1 on that platform. Are you *sure* this is a bug in the PDF and not in your pdf reader? Please read the PDF source file to find out. Mine has lines like /F2 1 Tf 12.00 0.00 -0.00 12.00 68.40 47.52 Tm (0) Tj /F2 1 Tf 12.00 0.00 -0.00 12.00 144.39 47.52 Tm (10) Tj

Re: [R] pdf device --- axes labels text disappeared?

/iaw, * [EMAIL PROTECTED] [2004-Oct-15 06:13 AKDT]: Dear R Wizards: Running R 1.9.1. on amd64. Promise- c(0,20,40); Expect- c(0, 20, 0.2*20+.8*40 ); # this omits printing numbers on the axes labels. pdf(file = bug.pdf ); plot(Promise, Expect, type=b, ylim=c(0,60)); dev.off(); Version

Re: [R] one more Rcmdr problem

I did experience the same problem. After installing R 2.0.0 patched and downloading the source for Rcmdr_0.99-12 from John Fox's Web page http://socserv.socsci.mcmaster.ca/jfox/Misc/Rcmdr/ and recompiling the package it now works. Ctrl-C, Ctrl-X, and Ctrl-V, in the text boxes of Rcmdr now work.

Re: [R] Testing for normality of residuals in a regression model

Berton Gunter wrote: Quite right, John! I have 2 additional questions: 1) Why test for normality of residuals? Suppose you reject -- then what? (residual plots may give information on skewness, multi-modality, data anomalies that can affect the data analysis). Because I want to know if my model

Re: [R] pdf device --- axes labels text disappeared?

Hi: Thank you everyone. You were all correct. This turns out to be a viewer bug---and even more likely a viewer installation bug, not an xpdf bug, much less an R bug. Apologies for having wasted everyone's time. Regards, /ivo --- ivo welch __

Re: [R] Problem with number characters

Scott, Has anyone suggested yet that some options might need adjusting in the Windows program that wrote the file? In my experience, CSV files from Windows applications are typically pure ASCII (though yours clearly isn't). Another possibility is that the program can export in plain ASCII with

[R] add-on package/windows installation probLem

Hi, I have tried to download the bootstrap.tar.gz from the website below: http://cran.r-project.org/src/contrib/Old/0.50/INDEX.html I have extracted the folder to a zip file...and then use the load packages from local zip files. I use R2.0.0.0 and windows NT, i get the error: 1: error -1 in

RE: [R] edit plots from the ADE4 package

Hi again, You could change the size font using the parameter clabel in both functions s.label and s.class. For the typical example in Fine, 1996 you could notice that slabel(acm$co,clabel=1) produce different results than

Re: [R] Testing for normality of residuals in a regression model

Berton Gunter wrote: Exactly! My point is that normality tests are useless for this purpose for reasons that are beyond what I can take up here. Thanks for your suggestions, I undesrtand that! Could you possibly give me some (not too complicated!) links so that I can investigate this matter

[R] Building package compatible w/ R v1.9.1 and R v2.0.0?

Hi, just in the process of updating my packages for R v2.0.0 and I have not had time to followed the R v2.0.0 discussions so maybe my questions have already been answered. A concern I have is that when creating packages they should be backward compatible with R v1.9.x for a while until all users

RE: [R] Testing for normality of residuals in a regression model

Berton Gunter wrote: Exactly! My point is that normality tests are useless for this purpose for reasons that are beyond what I can take up here. Thanks for your suggestions, I undesrtand that! Could you possibly give me some (not too complicated!) links so that I can investigate

RE: [R] one more Rcmdr problem

Brilliant. With R 2.0.0 patched and Rcmdr 0.9-12 all problems are gone. Thank you all, and a special thank you for the R Commander; it saves a lot of effort for teaching purposes. Best regards, Adrian -Original Message- From: Erich Neuwirth [mailto:[EMAIL PROTECTED] Sent: 15 octombrie

[R] SJava 0.65 and R 2.0.0 (slightly long)

Hi, I had some Java code that worked with SJava 0.65 under R 1.8.*. I'm trying to get it to work with R 2.0.0. My JVM is Sun 1.5.0 (running on Fedora Core 2) I downloaded and installed SJava 0.65 in my personal directory using R CMD INSTALL SJava-0.65.tar.gz -l ~/src/Rlibrary and before

Re: [R] add-on package/windows installation probLem

You downloaded a .tar.gz file and treated it as a zip file, for some reason not in any R documentation. On Fri, 15 Oct 2004, Nebahat Noyan wrote: Hi, I have tried to download the bootstrap.tar.gz from the website below: http://cran.r-project.org/src/contrib/Old/0.50/INDEX.html I have

Re: [R] add-on package/windows installation probLem

Nebahat Noyan wrote: Hi, I have tried to download the bootstrap.tar.gz from the website below: http://cran.r-project.org/src/contrib/Old/0.50/INDEX.html 0. For sure you are joking ... The package is 7.5 years old (16-Apr-1997) - too old to be installable by recent versions

Re: [R] Building package compatible w/ R v1.9.1 and R v2.0.0?

Henrik Bengtsson wrote: Hi, just in the process of updating my packages for R v2.0.0 and I have not had time to followed the R v2.0.0 discussions so maybe my questions have already been answered. A concern I have is that when creating packages they should be backward compatible with R v1.9.x for a

RE: [R] one more Rcmdr problem

Dear Adrian, Since I moved the discussion of this problem to the r-devel list, I should have reported back to r-help (1) that the problem was general to tlctk (not just the Rcmdr package), and (2) that Duncan Murdoch made a change to the patched Windows version of R 2.0.0 that made the problem go

[R] Mandrake rpm's for R 2.0.0

A binary RPM for Mandrake 10.0 official has been uploaded and will soon be available on CRAN, as well as the SRPM with which users of other flavours of Mandrake will be able to build (and hopefully contribute) their own RPM's. Please read the readme file if you have problems installing.

RE: [R] Testing for normality of residuals in a regression model

Let's see if I can get my stat 101 straight: We learned that linear regression has a set of assumptions: 1. Linearity of the relationship between X and y. 2. Independence of errors. 3. Homoscedasticity (equal error variance). 4. Normality of errors. Now, we should ask: Why are they needed?

RE: [R] Testing for normality of residuals in a regression model

Dear Kjetil, I don't believe that these are BLUS residuals, but since the last n - r effects are projections onto an orthogonal basis for the residual subspace, they should do just fine (as long as the basis vectors have the same length, which I think is the case, but perhaps someone can

RE: [R] Testing for normality of residuals in a regression model

Dear Federico, The problem is the same with GLS residuals -- even if the GLS transformation produces homoskedastic errors, the residuals will be correlated and heteroskedastic (with this problem tending to disappear in most instances as n grows). The central point is that residuals don't behave

RE: [R] Testing for normality of residuals in a regression model

Dear Andy, At the risk of muddying the waters (and certainly without wanting to advocate the use of normality tests for residuals), I believe that your point #4 is subject to misinterpretation: That is, while it is true that t- and F-tests for regression coefficients in large sample retain their

Re: [R] fidelity of generated raster images (R and perl)

Wrap-up: Different solutions were proposed for turning a matrix into a bitmap. Henrik Bengtsson [EMAIL PROTECTED] R.classes bundle from (http://www.maths.lth.se/help/R/R.classes/) MonochromeImage R.classes is still not compatible with 2.0.0 (minor modification to DESCRIPTION etc are needed) so

RE: [R] Testing for normality of residuals in a regression model

Hi John, Your point is well taken. I was only thinking about the shape of the distribution, and neglected the cases of, say, symmetric long tailed distributions. However, I think I'd still argue that other tools are probably more useful than normality tests (e.g., robust methods, as you

Re: [R] Testing for normality of residuals in a regression model

OK, I'll expose myself: I tend to do normal probability plots of residuals (usely deletion / studentized residuals as described by Venables and Ripley in Modern Applied Statistics with S, 4th ed, MASS4). If the plots look strange, I do something. I'll check apparent outliers for

[R] categorical varibles in coxph

Hello, I wonder when I do coxph in R: coxph( Surv(start, stop, event) ~ x, data=test) If x is a categorical varible (1,2,3,4,5), should I creat four dummy varibles for it? if yes, how can I get the overall p value on x other than for each dummy variable? Thanks Lisa Wang Princess Margaret

[R] Using accrualReport.s in Hmisc/Design

Hello, I apologize for the post if this is not the appropriate place. I have searched the R List archives and the Vanderbilt.edu website, but am not able to understand how to use the accrualReport.s and associated scripts available at the Vanderbilt site. It must be very simple, but I am at a

Re: [R] Using accrualReport.s in Hmisc/Design

Aric Gregson wrote: Hello, I apologize for the post if this is not the appropriate place. I have searched the R List archives and the Vanderbilt.edu website, but am not able to understand how to use the accrualReport.s and associated scripts available at the Vanderbilt site. It must be very

Re: [R] categorical varibles in coxph

On Fri, 15 Oct 2004, Lisa Wang wrote: Hello, I wonder when I do coxph in R: coxph( Surv(start, stop, event) ~ x, data=test) If x is a categorical varible (1,2,3,4,5), should I creat four dummy varibles for it? if yes, how can I get the overall p value on x other than for each dummy variable? No.

Re: [R] tree version 1.0-16

Kalaylioglu, Zeynep (IMS) wrote: Hi There is something weird going on with the tree package version 1.0-16. So I want to download an updated tree package. I found that R website has version 1.0-18. I want to try this one and see how the tree algorithm is performing. What is the best way to

RE: [R] combine many .csv files into a single file/data frame

On 15-Oct-04 bogdan romocea wrote: Dear R users, I have a few hundred .csv files which I need to put together (full outer joins on a common variable) to do a factor analysis. Each file may contain anywhere from a few hundred to a few thousand rows. What would be the most efficient way to

[R] power in a specific frequency band

Dear R users I have a really simple question (hoping for a really simple answer :-): Having estimated the spectral density of a time series x (heart rate data) with: x.pgram - spectrum(x,method=pgram) I would like to compute the power in a specific energy band. Assuming that frequency(x)=4

Re: [R] Testing for normality of residuals in a regression model

Liaw, Andy wrote: . . . . Also, I was told by someone very smart that fitting OLS to data with heteroscedastic errors can make the residuals look `more normal' than they really are... Don't know how true that is, though. Certainly true, since the residuals will be a kind of average, so the

[R] Problem Compiling R-2.0.0 on Linux Alpha

Hi folks: I'm trying to compile R-2.0.0 on my Alpha (164SX; alphapca56-unknown-linux-gnu) running RedHat Linux 7.2 with gcc3.4.1. Here is where the compile terminates: building package 'tcltk' mkdir -p -- ../../../library/tcltk/R mkdir -p -- ../../../library/tcltk/demo mkdir -p --

[R] 20,000 * 6 data values

Hello, Rusers: What is the maximum number of data R can handle? Or I have to use SAS? I am trying to do some microarray data analysis. But I am totally new. Did anyone use R to do microarray analysis? Many thanks, Sun __ [EMAIL PROTECTED] mailing

[R] sapply and loop

Dear all, I am doing 200 times simulation. For each time, I generate a matrix and define some function on this matrix to get a 6 dimension vector as my results. As the loop should be slow, I generate 200 matrice first, and save them into a list named ma, then I define zz-sapply(ma, myfunction)