search for: paran

Horn's parallel analysis method of deciding how many factors to retain has been implemented in R in the paran package. This update to paran adds an option to produce a graph of the analysis along the lines of what Horn presented in his seminal 1965 paper. It also makes a few minor improvements to the formatting of the program's output. Alexis Dinno _______________________________________________ R...

...sp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbspurl(http://mailimg.paran.com/skin/dotline3.gif) ; &nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&nbsp&am...

...date as something other than date (character). Let's say I am opening a CSV file, one of the columns of which is a date or time. How do I specify that when opening the file? Thanks for the help, Jim -- ------------------------------------- James J. Roper, Ph.D. Universidade Federal do Paran? Depto. de Zoologia Caixa Postal 19020 81531-990 Curitiba, Paran?, Brasil ===================================== E-mail: jjroper at gmail.com Phone/Fone/Tel?fono: 55 41 33611764 celular: 55 41 99870543 ===================================== Zoologia na UFPR http://zoo.bio....

Hi, I have a table of x rows and y columns. The table is huge and so i'd like to create a subset of the data containing rows where any of the y values are below a threshold, say 1e-4. Is there a simple way of doing this in R? thanks Rich

Hi All, I might misunderstood how cut works. But following behaviour surprises me. vv <- seq(150, 346, by= 4) cc <- cut(vv, 12) cc[vv == 330] Results [1] (330,346] I would have expected 330 to fall into (313,330] category. Can you please advice what do I do wrong? Many Thanks, Jussi Lehto -------------- next part -------------- Visit our website at http://www.ubs.com This message

Dear Sir/Madam, I got some problem about using R. I am a student, and need to use R a lot. But sometimes when I am working in computer lab, we can not install some particular pacakges by ourselves because of the previliges. Can the folder of Librabry of R, move to, for instance the personal folder, my documents... then we can install any specified package? Best regards,

Dear All, to produce output of several columns of a data frame, I tried to use lapply and also l_ply. In both cases, I would like to print a header line containing also the name of the respective column in the data frame. For example, I would like the following lapply(data.frame(a=1:3, b=2:4), function(x) print(deparse(substitute(x)))) to produce: [1] "a" [1] "b" and

Hallo All, I would like to apply a function to all permutations of variables in a dataframe (except the first). What is the best way to achieve this? I produce the permutations using: nvar <- ncol(dat) - 1 perms <- as.matrix( expand.grid(rep( list(1:0) , nvar ))[ , nvar:1] ) Thanks in advance Serguei Test-dataframe, comma-delimited: code,wav,w,area,gdp,def,pop,coast,milspend,agr

...nd it works fine. The firewall is disabled in ws2003, but I'm a linux user and it is possible that I'm missing some configuration to allow the connection. Any suggestion could be usefull. Thanks, Diego G. Pasqualin Brazil -- Diego G. Pasqualin UFPR - Universidade Federal do Paran? C3SL - Centro de Computa??o Cient?fica e Software Livre

...SecondEvent)/2 else DateFound Only the first half of the expression was ever evaluated. I hope I have not been to brief, and will certainly appreciate any help. Thanks, Jim -- James J. Roper Population Dynamics and Conservation of Terrestrial Vertebrates Caixa Postal 19034 81531-990 Curitiba, Paran?, Brasil ======================================= E-mail: jjroper at gmail.com Phone/Fone/Tel?fono: 55 41 33611764 celular: 55 41 99870543 Casa: 55 41 33857249 ======================================= Ecologia...

Is there a way to combine two columns within a data frame? Example data: id snp AL1 AL2 1500 30 A B 1510 30 A A 1520 30 A B This is what I would like: indv snp AL1AL2 1500 30 AB 1510 30 AA 1520 30 AB Any help is greatly appreciated. Alysta

I noticed that when I fit a linear model using 'lm' there is an attribute called "factors" that is added to the "term". It doesn't seem to appear for 'model.matrix', just 'lm'. I have been unable to find where it gets constructed or what it means? It looks like a two dimensional array that I may be able to use so I would just like to get some

I have an HTML form on the screen, that is connected to a database (standard Ruby way). I found that there are people that would rather not fill in this information, but would rather search Amazon for the info. Therefore, I added a link to the right of one of the fields that, when pressed, will search Amazon for what ever you typed into that field. So, how can I do that? At first I tried

Hi I am struggling to fit a non-linear trend using the likfit function in geoR. Specifically I want a sigmoidal function, something like SSfpl in the nls package to fit the trend. But it seems trend.spatial in geoR only works with lm or glm type models. Any ideas how I can specify the model to calculate the kriging parameters using REML, including the parameters of a sigmoidal trend function

I am attempting to develop a multiple regression model using selected model variables that should all be treated as numeric (mostly real) values. However, R considers one specific variable "mass" automatically to be of class "factor", probably because "mass" consists of integer values that are repeated. I now want to force R to treat "mass" as a numeric

...wd=2,lty='dotted') text.cca(dbRDA,dis='sites',col='black',cex=0.8,labels=FatoresSymbols$RepSimb) dev.off() ___________________________________ MSc. Rodrigo Aluizio Centro de Estudos do Mar/UFPR Laborat?rio de Micropaleontologia Avenida Beira Mar s/n - CEP 83255-000 Pontal do Paran? - PR - BRASIL Fone: (0**41) 3455-1496 ramal 217 Fax: (0**41) 3455-1105

...Can you tell me if this is a configuration issue on Dovecot? Thank you all -- Att. Rubens Aurelio Mascari ..................................... _ Linux & Internet Developer .... ICQ# 2730907 ............... ?v? Odara Internet LTDA............ (44) 3226-3662 ............. /(_)\ Maring? - Paran? .............. http://www.odara.com.br .... ^ ^ ______.-""-.__ QuEm TeM o CoNhEcImEnTo TeM o PoDeR __.-""-._______

I have defined a function to compute the value of a beta distribution of the second kind (the existing beta distribution of th stats package is the beta distribution of the first kind). It works perfectly for a single value, but I want to apply it to a vector of 22 000 values. I can use a loop for the calculation of each value but it runs very very slowly. So, what can I change ? Hers's the

Hello, i know that i have to use optim, but i'm confused how its possible maximize the sum over all l[i] and get the optimized max(LL), r and alpha? LL <- function(trans,time){ for(i in 1:length(trans){ l[i] <- log(lgamma(r+trans[i] - gamma(r+1)*(alpha/alpha+t[i]))**r)*(t[i]/alpha+t[i]))**trans[i] } return(sum(l)) } i'm confused how i have to set r and alpha and i

Hi I have recently been analyzing birthweight data from a clinic. The data has obvious defects in that there is digit preference on certain weights making them overrepresented. This shows as spikes in the histogram on certain well rounded weights like 2, 2.5, 3, etc. I would like to show this to government officials but can't figure out how I should present the finding in an easy to