similar to: Weighted Principle Components analysis

Hello, I am using function princomp from the package psych. I have my principle component object mypc: mypc <- princomp(covmat=mycor) plot(mypc) # shows me a screeplot Question: how could I actually see the values displayed in the screeplot. I don't mean on the graph - I just want to know the actual value for each component (e.g., 10, 3.2, 1.8, etc.) I need to know how much variance,

Hello! I am running principle components analysis using princomp function in pacakge psych. mypc <- princomp(mydataforpc, cor=TRUE) Question: I'd like to use pairwise deletion of missing cases when correlations are calculated. I.e., I'd like to have a correlation between any 2 variables to be based on all cases that have valid values on both variables. What should my na.action be in

Hello! I have a data frame mysample (sorry for a long way of creating it below - but I need it in this form, and it works). I regress Y onto X1 through X11 - first without weights, then with weights: regtest1<-lm(Y~., data=mysample[-13])) regtest2<-lm(Y~., data=mysample[-13]),weights=mysample$weight) summary(regtest1) summary(regtest2) Then I calculate Durbin-Watson for both regressions

Hi all I'm using Samba version 3.0623 on my rhel3 with kernel 2.4.21 I mounted a windows partition with the command "smbmount". Now I have the problem, that I cannot copying files which are larger than 2GB. The cp command starts and everything works fine until the 2GB are reached. Then the bash gives me the error: "The maximum file size is exceeded" the

Hello list, I am configuring NUT on a Dell workstation that is connected to an APC Back UPS Pro 1300, and I am running into an issue where starting everything manually seems to be successful, but fails when using the init.d scripts. The short story is that I get the following error when trying to check the status of the UPS: [root at mypc:/etc/nut]# upsc apc_pro1300 at localhost Error:

Dear Tomas, Thanks for your email. This can be easily reproduced if there is access to a shared folder. I've reproduced this behaviour on both Windows 7 and Windows 10. Steps to reproduce: 1. Ensure Offline Files is turned on within Windows 2. Using Windows Explorer, browse to a folder shared on a network using a UNC path, e.g. \\mypc\myshare\ 3. Create a test file, e.g. test.txt 4. Within

Hello Rowloand, ok, I just keep "default" for the mounting. I changed the path of my share : [homes] comment = User Home Folder path = /home/%U read only = no follow symlinks = yes wide links = yes browseable = no hide files = /desktop.ini/$RECYCLE.BIN/ I need Wide Links because in the share there are symbolics links to outside the share. Anyway, I still get the same error with all

I am not sure if this is an R-users question, but since most of you here are statisticians, I decided to give it a shot. I am using the lm() function in R to fit a dependent variable to a set of 3 to 5 independent variables. For this, I used the following commands: >model1<-lm(function=PBW~SO4+NO3+NH4) Coefficients: (Intercept) SO4 NO3 NH4 0.01323 0.01968

Hi. I was comparing the components from princomp's loadings and the eigen given the same input. I found that the sign of componenets (+/-) are opposite between the two components (from princmop and eigen) but the magnitudes are identical. Why? Thanks! [[alternative HTML version deleted]]

I am new to R. I just tried to recreate in R (using sem package and the identical input data) a solution for a simple measurment model I have found before in LISREL. LISREL had no problems and converged in just 3 iterations. In sem, I got no solution, just the warning message: "Could not compute QR decomposition of Hessian. Optimization probably did not converge. in: sem.default(ram =

Le 05/01/2020 ? 14:27, Rowland penny via samba a ?crit?: > On 05/01/2020 11:36, Ni Co wrote: >> Hello Rowloand, >> >> ok, I just keep "default" for the mounting. >> >> I changed the path of my share : >> >> [homes] >> comment = User Home Folder >> path = /home/%U >> read only = no >> follow symlinks = yes >> wide

Hi I am struggling with nested random effects and hope someone can help. I have individuals (ID) who are nested within families (FAM). I want to model an outcome variable, and take account of the intercorrelation of individuals within each family. I think this amounts to two random effects, one nested within the other. How can I model this in R? So far I have tried using the

Hello, I finally run my measurement model in sem - successfully. Now, I am trying to print out the path diagram that is based on the results - but for some reason it's not working. Below is my script - but the problem is probably in my very last line: # ANALYSIS OF ANXIETY, DEPRESSION, AND FEAR - LISREL P.31 library(sem) # Creating the ANXIETY, DEPRESSION, AND FEAR intercorrelation matrix

I think you should be more suspicious of yourself, Dimitri. A letter T variable can easily arise in the problem domain when you are not thinking of logical values at all, at which point your cavalier use of T as a synonym for TRUE can suddenly become a bug. -- Sent from my phone. Please excuse my brevity. On July 27, 2017 8:18:03 AM PDT, Dimitri Liakhovitski <dimitri.liakhovitski at

Dear R-list, maybe some of you could point me in the right direction: Are you aware of any FREE Fortran or Java libraries/actual pieces of code that are VERY efficient (time-wise) in running the regular linear least-squares multiple regression? More specifically, I have to run small regression models (between 1 and 15 predictors) on samples of up to N=700 but thousands and thousands of them. I

Hello! I have a column in my data frame that I have to split: I have to distill the numbers from the text. Below is my example and my solution. x<-data.frame(x=c("aaa1_bbb1_ccc3","aaa2_bbb3_ccc2","aaa3_bbb2_ccc1")) x library(stringr) out<-as.data.frame(str_split_fixed(x$x,"aaa",2)) out2<-as.data.frame(str_split_fixed(out$V2,"_bbb",2))

Just a thought: Did you try na.rm = TRUE in case you have an object named "T" in scope? -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Thu, Jul 27, 2017 at 7:49 AM, Dimitri Liakhovitski <dimitri.liakhovitski at

Hello! I am trying to understand how ggplot2's geom_bar treats NAs. The help file says: library(ggplot2) ?geom_bar na.rm: If FALSE, the default, missing values are removed with a warning. If TRUE, missing values are silently removed. I am trying it out: md <- data.frame(a = c(letters[1:5], letters[1:4], letters[1:3], rep(NA, 3))) str(md); levels(md$a) ggplot(data = md, mapping = aes(x =

Hello! s<-"start"; e<-"end" middle<-as.character(c(1,2,3)) I would like to get the following result: "start 123 end" or "start 1 2 3 end" or "start 1,2,3 end" How can I avoide this (undesired) result: paste(s,middle,e,sep=" ") Thank you! -- Dimitri Liakhovitski Ninah Consulting www.ninah.com

?hanks for the advice, Jeff. Will keep it in mind. But I am anal - I shy away from using letters and words that "look familiar" to me in R (such as mean, sd, T, etc.) But still, it's a good advice. On Thu, Jul 27, 2017 at 11:53 AM, Jeff Newmiller <jdnewmil at dcn.davis.ca.us> wrote: > I think you should be more suspicious of yourself, Dimitri. A letter T > variable can