similar to: For loops in R

Hi, I am trying to call a C function, that I wrote to parse a flat file, into R. The argument that will go into this function is an input file that I need to parse and write the desired output in an output file. I used some hit and trial approach but i keep on getting the "file not found" or "segmentation fault" error. I know that the error is in passing the argument but I

Hello all, I have two matrices, pop and pop2, each the same number of rows and columns that I want to compare for equality. I am concerned about efficiency in this operation. I've tried a few things without success so far. Doing something simple like: if (pop==pop2) { cat('equal') } else { cat('NOT equal') } results in the warning: 1: In if (pop == pop2) { : the

Dear all, I'm quite sure that this is a stupid question, but I'll ask anyway. I want to perform an ANCOVA with two continuous factors and three categorical factors. Plant population growth rate (GR) = dependent variable Seed reduction due to herbivory (SR) = continuous explanatory variable Herbivore species (HS, 2 levels) = categorical explanatory variable Population (Pop, 24 levels) =

Hello all, I have the following function call to create a matrix of POP_SIZE rows and fill it with bit strings of size LEN: pop=create_pop_2(POP_SIZE, LEN) I have 3 questions: (1) If I did keep_pop[1:POP_SIZE] == pop[1:POP_SIZE] to keep a copy of the original data structure before manipulating 'pop' potentially, would this make a deep copy or just shallow? Ie

I would like to run a for loop with an index going from 0 to 499 but the following seems to miss out the first value: C <- 499 for (i in 0:C) The alternative is: C <- 500 for (i in 1:C) { #Then every time I use i, I replace it with i-1 } Is this a good way to do it or is tere a better way? Thank you, ThomasThis message and any attachment are intended solely for the addressee and may

I am using *Genetic Algorithm* to maximize some function which use data. I use GA package in R for this ( https://cran.r-project.org/web/packages/GA/index.html) Below is my code library(GA) set.seed(1) Dat = data.frame(rnorm(1000), matrix(rnorm(1000 * 30), nc = 30)) Fitness_Fn = function(x) { return(cor(Dat[, 1], as.matrix(Dat[, -1]) %*% matrix(x, nc = 1))[1,1]) } ga(type =

Dear Daniel As far as I can see you have re-generated the data before calling ga() so the data is not just a permutation of the first set. Michael On 07/05/2025 15:36, Daniel Lobo wrote: > I am using *Genetic Algorithm* to maximize some function which use data. > > I use GA package in R for this ( > https://cran.r-project.org/web/packages/GA/index.html) > > Below is my code

Dear All I new to using R and am struggling with some matrix multiplication. I have two matrices, one containing random numbers, these are multiplied together to get another matrix which is different each time. When I put in another for loop to repeat this process a multiple times the matrices are all the same. I?m sure there is a way to keep the randomness of the different matrices but I think

Dear R Users, I am a novice learner of R software. I am working with Genetic Matching - GenMatch(), but I am getting an Error message as follows:  Error in GenMatch(Tr = Tr, X = X.binarynp, BalanceMatrix = BalanceMatrix.binarynp,  :    GenMatch(): input includes NAs Could you please suggest me correcting the above problem? My GenMatch command is, > gen1 <- GenMatch(Tr = Tr, X =

Dear Sir: When I try to save large and very complex recursive objects with some components containing models (such as the output from arima or lm ), the resulting file sizes increase by 4 meg per save. example directory: amex 8 meg argentina 12 meg australia 16 meg ... Moreover, I am unable to read these file objects back into R. I note that readBin and writeBin

-------------------------------------------------- From: "Ahoussou Sylvie" <sylvie.ahoussou at antilles.inra.fr> Sent: Friday, September 12, 2008 9:48 AM To: "Thomas Lumley" <tlumley at u.washington.edu> Subject: Re: [R] Complex sampling survey _ Use of survey package > Thanks for your answer > > I think I made a mistake when I recopied the 5 first rows of

Version 3.9 of the survey package is now on CRAN. Since the last announcement (version 3.6-11, about a year ago) the main changes are - Database-backed survey objects: the data can live in a SQLite (or other DBI-compatible) database and be loaded as needed. - Ordinal logistic regression - Support for the 'mitools' package and multiply-imputed data - Conditioning plots,

Version 3.9 of the survey package is now on CRAN. Since the last announcement (version 3.6-11, about a year ago) the main changes are - Database-backed survey objects: the data can live in a SQLite (or other DBI-compatible) database and be loaded as needed. - Ordinal logistic regression - Support for the 'mitools' package and multiply-imputed data - Conditioning plots,

I have made the 3-d case do the same as the vector case, which is what the C code clearly intended (a goto label was in the wrong place). This leaves the bigger question of the right thing to do. I note that data frames give an error when any indices are NA. -thomas On Fri, 3 Sep 2004 ripley@stats.ox.ac.uk wrote: > Apart from the inconsistencies, there are two clear bugs here: > > 1)

https://bugzilla.mindrot.org/show_bug.cgi?id=3722 Bug ID: 3722 Summary: check ssh key (and other parameters) before running the proxy command Product: Portable OpenSSH Version: 9.8p1 Hardware: Other OS: Linux Status: NEW Severity: enhancement Priority: P5 Component:

Dear all, I had a R segmentation fault, and then invoked debug mode and ran step by step. When I reached "terms(Y~X1*X2*...*X16)", I would then have "segmentation" fault. However, if I just ran this under regular "R interactive" mode, it would be fine though taking long time. My questions are: 1. Is there a known limit of terms for a formula? 2. Why does the

Dear all, I had a R segmentation fault, and then invoked debug mode and ran step by step. When I reached "terms(Y~X1*X2*...*X16)", I would then have "segmentation" fault. However, if I just ran this under regular "R interactive" mode, it would be fine though taking long time. My questions are: 1. Is there a known limit of terms for a formula? 2. Why does the

Using CVS from 7/12/04 and trying to get the voicemail broadcast feature to work. Voicemail.conf has [mycontext] 3722 => 1234,BroadCast Test,,,cc=*@mycontext . then many other voicemail boxes. ----- whenever I leave voicemail at box 3722, only box 3722 gets the voicemail. It is not expanding it to other voicemail boxes in the [mycontext] context. Even if I replace the cc= line with

I was experiencing a strange pattern of slowdowns when using sample.int(), where sampling from a one population would sometimes take 1000x longer than taking the same number of samples from a slightly larger population. For my application, this resulted in a runtime of several hours rather than a few seconds. Looking into it, I saw that sample.int() is hardcoded to switch algorithms when the

Dear R Users, I am a novice learner of R software. I am working with Genetic Matching - GenMatch(), but I am getting an Error message as follows:  Error in GenMatch(Tr = Tr, X = X.binarynp, BalanceMatrix = BalanceMatrix.binarynp,  :    GenMatch(): input includes NAs Could you please suggest me correcting the above problem? My GenMatch command is, > gen1 <- GenMatch(Tr = Tr, X =