similar to: iterative output to file by row

Hi, I am running R 2.0.1.1. on Windows. It is a Dell Dimension with a 3.2 Ghz Processor and 4Gb RAM. When using the ReadAffy() function to read in 97 arrays, I get the below error messages: Error: cannot allocate vector of size 393529 Reached total allocation of 1024Mb: see help(memory.size) When I use the comman "memory.limit(size=4000)" to increase the memory size to the

Hello, My R skills are somewhere between novice and intermediary, and I am hoping that some of you very helpful forum members, whom I've seen work your magic on other peoples' problems/questions, can help me here. I have a matrix with the following format: (i) individual plants comprising many different genotype groups (i.e., a plant is genotype 1 or genotype 2 or genotype 3, etc). The

Full_Name: Christina Merz Version: R version 2.5.0 (2007-04-23) OS: mingw32 Submission from: (NULL) (213.70.209.132) R> version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status

I'm using R 1.8.1 on Red Hat Linux 9. I've worked out the linear model decomposition by hand. Using "aov" with "Error" to fit a split-plot model, I get very different results depending on whether I use the "data" argument for "aov" or not. When I use the "data" argument, the ANOVA table from "summary" is correct but the

Hi everyone I'm using Tukey HSD as post-hoc test following a lme analysis. I'm measuring hemicelluloses in different species treated with three different CO2 concentrations (l=low, m=medium, h=high). The whole experiment is a split-plot design and the Tukey-function from the package multcomp is suitable for lme-analysis with random factors. The analysis works fine but I get a non

Hi guys, I have a list of 250 numbers as a result of using the ?by function! List of 246  $ 0   : num [1:28] 22 11 31...  $ 1   : num [1:15] 12 14 9 ... .. .. ..  - attr(*, "dim")= int 250  - attr(*, "dimnames")=List of 1 The problem is that each list of 250 has different length! I would like to get the values of each list in a vector like vector_0 = (22,11,31,..), is this

Hi guys, I still did not solve my problem properly! I have to compare the values of two lists of 250 numbers as a result of using the ?by function! List1 of 250  $ 0   : num [1:28] 22 11 31...  $ 1   : num [1:15] 12 14 9 ... .. .. ..  - attr(*, "dim")= int 250  - attr(*, "dimnames")=List of 1 List2 of 250  $ 0   : num [1:24] 20 12 22...  $ 1   : num [1:17] 11 12 19 ... .. ..

Dear all! W2k, R 2.5.1 I am working with an ongoing malting barley variety evaluation within Sweden. The structure is 25 cultivars tested each year at four sites, in field trials with three replicates and 'lattice' structure (the replicates are divided into five sub blocks in a structured way). As we are normally keeping around 15 varieties from each year to the next, and take in 10 new

I'm trying to export a large data frame to a text file for permanent storage. The only thing I could find was the treeglia Package but that didn't work. Any suggestions? Thanks! Christina D Smith PhD Student, GRA Statistics Department Kansas State University

Dear all, I have a time series of daily measurements that starts like this: KWhourly[1:10,] Date Time Value 01/01/13 00:00:00 1.6 01/01/13 01:00:00 1.6 01/01/13 02:00:00 1.6 01/01/13 03:00:00 1.6 01/01/13 04:00:00 1.6 01/01/13 05:00:00 1.6 01/01/13 06:00:00 1.6 01/01/13 07:00:00 1.6 01/01/13 08:00:00 1.6 01/01/13 09:00:00 1.6 I am trying to get average daily,

Hi, I'm installing an Asterisk box with a TDM2400P in Brazil. I can make analog phones work while lines are not working. Since I do not know anything about brazilian lines, is there anybody who can tell me what is wrong/missing in my conf files (below)? TIA Giorgio _zaptel.conf:_ fxoks=9-16 fxsks=17-24 defaultzone=br loadzone=br* * _zapata.conf:_ context = inbound_zap echocancel = 128

OK, that was a stupid error in the le_int prototype. It's fixed in svn now. Jean-Marc Le mercredi 28 septembre 2005 ? 11:53 +1000, Jean-Marc Valin a ?crit : > OK, it seems like I screwed up the wav header code between 1.1.8 and > 1.1.9. Everything seems fine (playing to soundcard also works), except > that the length reported for the wav is too short, which is why it stops >

Dear All ~ I have * setup & running ok (with two Wildcard X100P's to PSTN). I also have two analog phones connected into same through a SIPURA 2000. These work fine, except that when I call out through PSTN & try to send DTMF tones to (say) a remote PBX to dial an extension, the gain seems to go wild (high), and the DTMF tones are not recognized at the other end. I tried setting the

Hi, R core team I am using the function {aidsEst} in package [micEcon] to do an AIDS model now. So far, everything is good. But I want to test the auto correlation and heteroskedasticity of the individual equation from AIDS demand system. How can I return the individual equation? PS: serial correlation test is {bgtest} in package [lmtest] and heteroskedasticity is {bptest} in package

Hi, Does anyone have the data of the rsync transfer speed vs. other tools like ftp? Or where is the place to find such information that compare the transfer speed between different tools? If I were to transfer a single large file, say 500Mega Byte, previously non-existent on the remote host, what would be the fastest tool I can choose for transfer? Thanks, Christina Li

Dear R Help, I wonder the way to show the source code of [vcov] command. Usually, it can show the source code after input the command and enter. But for [vcov], it shows function (object, ...) UseMethod("vcov") <environment: namespace:stats> I appreciate for your help. Best wishes. Christina [[alternative HTML version deleted]]

I'm attempting to run an LME model with random intercepts for location and have received an error message. The model and error message are given below: > mixed.distance <- lme (fixed= log(min5inpnc) ~ i93centerline, random=~1|location, data = nosmokecomp) Error in y[revOrder] - Fitted : non-conformable arrays I believe that my dataset may have dimensions that are not working with

Dear R Core Team, I am using package {urca} to do cointegration and estimate ECM model, but I have the following two problems: (1) I use ca.jo() to do cointegration first and can get the cointegration rank, alpha and beta. The next step is to test some restrictions on beta with blrtest(),bh5lrtest(), and bh6lrtest(). But none of them can add restrictions on all the cointegration

Dear All, I am trying to iterate an iterative process (i know R is not the best place for so much looping but i have to so tough!) and store the resulting data from each iteration to then be graphed. In simple form the script looks like this for(i in tracks) { a <- (subset(data, data$Id==i)) result.dataset <- function(a, other.arguements) ### this function is the other

Hello, I am a novice R user and am having difficulty retrieving the values from 21 iterations of the R function integrate. The only way I have found is to do a write.table and then a read.table as shown in the code below. I would rather capture the 21 values inside the braces ( sapply might work, but I can't set it up without getting an error in function) so I could compute the