similar to: 3-dimensional looping Q.

Hi everyone- I have used several methods to calculate principal components rotated using the varimax procedure. This is simple enough. But I would like to calculate the % of variance explained associated with each PC before and after rotation. factanal returns the % of variance explained associated with each PC but I cannot seem to get it to change after rotation. Many thanks for your

Hello again- Thanks to Prof. Ripley for responding to my previous question. I would like to clarify my question using sample code. I will use some sample code taken from ?prcomp Again, I would like to compare the % variance explained by each PC before and after rotation. < code follows > data(USArrests) pca = prcomp(USArrests, scale = TRUE) # proportion variance explained by each

Hi. is there a way to transform a table that has numeric and categorical variables into a matrix? thanks Maria Mercedes Gavilanez Department of Biological Sciences 107 Life Sciences Building Louisiana State University Baton Rouge, LA 70803 (225)578-4284

Dear fellow R users, I am now investing time in learning how to use compiled C code to produces functions that can be used in R. I am just starting, and there is much that I need to learn, so I have a question that might be straight forward. I am learning how to use function in the C library GSL (gnu scientific library), to write C code, that I then plant to use in R. Is there any problem in

Does anyone know of a literature reference, or a piece of code that can help me calculate the amount of variation explained (R2 value), in a regression constrained to have a slope of 1 and an intercept of 0? Thanks! Sebastian J. Sebastián Tello Department of Biological Sciences 285 Life Sciences Building Louisiana State University Baton Rouge, LA, 70803 (225) 578-4284 (office and lab.)

Dear R-List, My questions concerns missing values. Specifically, is is possible to use different "types" of missingness in a dataset and not a one-size-fits-all NA? For example, data may be missing because of an outright refusal by a respondent to answer a question, or because she didn't know an answer, or because the item simply did not apply. In later analysis it is sometimes

I have a data.frame as follows: a 3 a 2 a 1 b 3 b 2 c 2 c 3 c 1 c 1 Each factor (a, b, c) should be monotonically decreasing, notice that factor 'c' is not. I could use some help to figure out how to form a logical structure (mostly just syntax), that will check each 'next value' for each factor to see if it is less than the previous value. If it is less than the

I have a really long functions, and at the end of the function, I am using a if statement to tag certain keywords based on whether they have certain values contained in them. However, the if statement doesn't seem to work. When I had split up the commands into various functions, it worked fine, but I'm not sure what going on now that it's combined into a single function. myfunc

Dear R users, I am looking for more efficient way to compute the followings -------------------------------------------------------------------------- a <- matrix(c(1,1,1,1,2,2,2,2),4,2) b <- matrix(c(1,2,3,4),4,1) Eventually, I want to get this matrix, `c`. c <- matrix(c(1/1,1/2,1/3,1/4,2/1,2/2,2/3,2/4),4,2) --------------------------------------------------------------------------

Hi all, I just svn-updated the 2.5 branch on my machine and I noticed this warning during the build. *** Warning: Linking the shared library /home/maurice/installation/llvm/Debug/lib/libLTO.la against the non-libtool *** objects /home/maurice/installation/llvm/Debug/lib/LLVMCppBackend.o /home/maurice/installation/llvm/Debug/lib/LLVMMSIL.o /home/maurice/installation/llvm/Debug/lib/LLVMCBackend.o

Maurice Gittens wrote: > Hi all, > > I just svn-updated the 2.5 branch on my machine and I noticed this > warning during the build. > > *** Warning: Linking the shared library > /home/maurice/installation/llvm/Debug/lib/libLTO.la against the non-libtool > *** objects /home/maurice/installation/llvm/Debug/lib/LLVMCppBackend.o >

Hi all, On my Fedora 10 x86-64 the current svn dies in the following way: /home/maurice/installation/llvm-dev/utils/TableGen/Debug/TGLexer.o: In function `llvm::TGLexer::getNextChar()': /home/maurice/installation/llvm-dev/utils/TableGen/TGLexer.cpp:71: undefined reference to `llvm::TGSourceMgr::FindBufferContainingLoc(llvm::TGLoc) const'

Dear list, I have two vectors: x <- c("one","two") y <- paste(rep(x,2),"blah") I want to replace all occurrences of each element of x in y with something else, so that y looks like this: y [1] "something else blah" "something else blah" "something else blah" [4] "something else blah" I can do this using a loop: for (

Dear R Listers- I am attempting to analyse the survival of seeds in cages (exclosures) that differ in their permeability to rainforest mammals. Because I did not observe the moment of seed disappearance, my data is interval censored. This limits my options for analysis (as I understand it) to survreg, in the survival package. Because I repeated the experiment in 8 sites, I have a random

Dear R Gurus - Please help me to make a multiple-panel plot using xYplot in which the y-axis of each panel is dependent only on the data of that plot. Here is the dataset: response date mean Lower Upper 1 density 1 271.0000000 249.2500000 289.0000000 2 density 4 296.7500000 280.5000000 313.0000000 3 recruits 1 271.0000000 249.2500000

Dear all, I have been using the function gls in the package nlme in R to fit some spatial regressions (as described in Dormann et al.). However, I have been struggling trying to find a way to calculate a measure of effect size from these models, so I wanted to know if any of you had an idea on how to do this. More precisely, I am producing a multiple model with an exponential correlation

Hello R-listers, I am using xtable to create some LaTeX tables for me to Sweave into my document. I am now trying to set the number of digits that will be reported in several _rows_ in my table. The 'digits' option in xtable() sets the number of digits to display for specified columns. Am I missing another obvious option to do the same thing for rows? A quick RSiteSearch didn't turn

On 2009-03-25 15:48, Maurice Gittens wrote: > Hi all, > > On my Fedora 10 x86-64 the current svn dies in the following way: > > /home/maurice/installation/llvm-dev/utils/TableGen/Debug/TGLexer.o: In > function `llvm::TGLexer::getNextChar()': > /home/maurice/installation/llvm-dev/utils/TableGen/TGLexer.cpp:71: > undefined reference to >

Rsync 3.2.7 is running on the Gentoo computer, which doesn't have a version, other than it's "current". I'm running the script from this computer. Rsync 3.1.2 is on the source computer, where the files come from, which is Ubuntu 18.0.4.6. I'm copying to a CIFS share mounted on the Gentoo computer. The rsync scripts are all similar to this one: /usr/bin/rsync -v -a

Hi, See below I reply your message for <https://stat.ethz.ch/pipermail/r-help/2008-April/160966.html>[R] predict.glm & newdata posted on Fri Apr 4 21:02:24 CEST 2008 You say it ##works fine but it does not: if you look at the length of yhat2, you will find 100 and not 200 as expected. In fact predict(reg1, data=x2) gives the same results as predict(reg1). So I am still looking for