similar to: Iterations of random sampling

I'm using .C() and .External() and have no problems sending integers, reals or strings from R to C. Nor do I have problems sending integers or reals back from C to R. But I'm pulling my hair out trying to set a string value in a C function and then sending it back from C to to R. I've searched the usual sources and tried various casts, macros and allocation schemes, but I'm

Hello there. Happy new year for everyone! I need help with a table. This table contains 300 rows and 192 columns. Being the first column the ID of my samples that can have several observations. I need to generate e NEW table that contains a single ID with the sum of the observations by columns: For example: Example ID?? A??? B? ? C? ? D? ? E? ? F? ? G.... 191 columns a1?? 0??? 0??? 0? ? 1???

Hello everybody, I'm a new user of R and I'm working right now with the ACE function from the acepack library. I Have a question: Is there a way to predict new responses using ACE? What I mean is doing something similar to the following code that uses PPR (Projection Pursuit Regression): library(MASS) x <- runif(20, 0, 1) xnew <- runif(2000, 0, 1) y <- sin(x) a <- ppr(x, y,

Hi, I am a long time SPSS user but new to R, so please bear with me if my questions seem to be too basic for you guys. I am trying to figure out how to analyze survey data using logistic regression with multiple imputation. I have a survey data of about 200,000 cases and I am trying to predict the odds ratio of a dependent variable using 6 categorical independent variables (dummy-coded).

I have a data frame of ~200 columns and ~20,000 rows where each column consists of binary responses (0,1) and a 9 for missing data. I am interested in finding the columns for which there are fewer than 100 individuals with responses of 0. I can use an apply function to generate a table for each column, but I'm not certain whether I can subset a list based on some criterion as subset() is

I have built and uploaded to CRAN the Carbon version of R for MacOS/MacOSX. You'll find the base set (rm151.sit) and the recommended packages (recommended.sit) separately or all in one single file (rm151_FULL.sit) on CRAN and their mirrors soon. I have also built the Bioconductor suite of packages (corresponding to the Bioconductor_1.0.tgz file) as it seems to be one of the most

I have built and uploaded to CRAN the Carbon version of R for MacOS/MacOSX. You'll find the base set (rm151.sit) and the recommended packages (recommended.sit) separately or all in one single file (rm151_FULL.sit) on CRAN and their mirrors soon. I have also built the Bioconductor suite of packages (corresponding to the Bioconductor_1.0.tgz file) as it seems to be one of the most

Hi!!! I'm trying to implement a SUS server (local windows update) I-ve already solved my problem to modify registry in order to establish windows update configuration. Now i've a quiestion all my client logon my samba PDC enusing a netconfig.pol isend all configuration everithing works ok. But i've a trouble, my users don't have privileges to install anything and i don't

dear all, i'm trying to reproduce an spss-anova in R. It is an 2x3x3 repeated measures desingn with repeated contrasts. In R i've coded a contrast matrix for all factors and made a split in the aov summary - but I can't get the repeated interaction contrasts. The output from SPSS looks like this: TaskSw * CongNow * CongBefore: SS df Mean Square F Sig. 1 vs. 2 1 vs. 2 1 vs. 2

Dear useRs, I have two different length vectors: one column (1...m) and one row vector (1...n): 20 40 20 60 5 4 2 Now I have to calculate ratios between column vector elements and each row vector elements: 4 5 10 8 10 20 4 5 20 15 12 30 Thank's in advance for any suggestions, Andrej

Dear mailing list I have a data that have 20,000 rows and 20 columns. Io wonted to extract the 10th row only. Example the 10th, 20th, 30th 40th…..20000 th. can you please help me how do I do that.Than kyou. Example is below. Inpute: AG GG GG AG CC CC CC CC CT CC CT CT GG GG GG GG CC CC CC CC GG GG GG GG CC CC CC CC GG CG CG GG GG GG GG GG *CC CC CC CC* AA AG AG AA AA AA AA AA GG AG AG GG GG AG AG

myDF: d1 d2 d3 d4 d5 -0.166910351 0.022304377 -0.00825924 0.008330689 -0.000925938 -0.166910351 0.022304377 -0.00825924 0.008330689 -0.000925938 -0.166910351 0.022304377 -0.00825924 0.008330689 -0.168225938 -0.166910351 0.022304377 -0.00825924 0.008330689 -0.168225938 -0.166910351 0.022304377 -0.00825924 0.008330689 -0.168225938 -0.166910351

Greetings, I am recoding a dummy variable (coded 1,0) so that 0 = 2. I am using the line sciach$dummyba[sciach$ba==0] <- 2 I notice that it creates a new column dummyba, with 0 coded as 2 but with 1's now coded as NA. Is there a simple way around this in the line I am using, or do I need to have an additional line sciach$dummyba[sciach$ba==1] <- 1 Thanks in advance. David

When I try citypro <- sas.get('c:/ndri/cvar/data', member = 'citypro2') I get the following errors and warnings: 'sas' is not recognized as an internal or external command, operable program or batch file. 'less' is not recognized as an internal or external command, operable program or batch file. Error in sas.get("c:/ndri/cvar/data", member =

Hi, Im' trying to put a SUS server on my network, iv'e already installed it, i've a PDC using samba, and my clients get registry configuration to get the updates using a Netconfig.pol file this parts works ok. So I'm having troubles with my clients, they don't get the updates from the server when they are into the samba domain it looks like they don't see the sus server, i

Hi, I'm trying to make a boxplot with the data at the end of the message, and when I try to execute the command >boxplot(Diatoms) (or for any other field instead of "Diatoms") I get the following error message: Error in oldClass(stats) <- cl : adding class "factor" to an invalid object Any advice would be much appreciated. Thanks a lot, Miriam Date

So I am very new to R. Have been using python for a project and need to calculate the correlation coefficient matrix for my data set. the data is in the range of 10-15 observations of 230,000 variables. ie the correlation matrix would be 230,000X230,000 Using python and the numpy.corrcoef() I run out of memory if I try to do this with more than ~30,000 variables. I was able to load the data into

Hello, How ls() only functions or only data objects (basically anything other than functions) such as data.frame, numeric ...? John [[alternative HTML version deleted]]

Hi, I am trying to create a power curve to show how the power of a t-test varies depending on the mean. Any ideas how I should go about this? Louisa _________________________________________________________________ [[alternative HTML version deleted]]

Hi, I'm trying to write a function that will divide a given range of numbers into 3 sets using sample(), without repetition. Currently I'm trying this approach: r <- 1:10 s1 <- sample(r,size=3) Next, I want to remove the selected elements from r and sample() from the remainder. r <- r[ -(r=s1) ] s2 <- sample(r,size=3) When I go to remove the elements contained in s2 from