similar to: How-To construct a cov list to use a covariance matrix in factanal?

I have a data frame that looks like this: > print(df) V2 V3 V4 V5 V6 V7 V8 V9 V10 V11 V12 1 FN 8.637 28.890 31.430 31.052 29.878 33.215 32.728 32.187 29.305 31.462 2 FP 19.936 30.284 33.001 35.100 30.238 34.452 35.849 34.185 31.242 35.635 3 TN 0.000 17.190 16.460 21.100 17.960 15.120 17.200 17.190 15.270 15.310 4 TP 22.831 31.246 33.600 35.439 32.073

I try to calculate descriptive statistics for one of the variables in the data frame, however command sapply calculates these statistics for every value of the variable separately. How to make it calculate range (as well as other statistics) for all column? Here are commands and results: > as1$trust [1] 5.957510 5.888664 6.168135 6.419472 5.668796 6.026923 6.456721 7.017946 5.294411

Is that an error that prevents anything from working? It's probably expected... On Thu, Jul 7, 2016 at 2:07 PM, abcd <abcdeluxe at web.de> wrote: > Thanks for the quick help! > > Indeed a kernel update resolved the issue, there one problem left > though: I get an error: > > Xorg.0.log: Couldn't load sub module "fb". > > I can't find out what

Hi! I have some questions about MARS model's coefficient of determination. I use the MARS method in my master's thesis and I have noticed some problems with the MARS model's R^2. You can see the following example that the MARS model's R^2 is too big when i have used mars() -function for MARS model building, and when I have made MARS-model using a linear regression, it

Hello R users, For the last few days I am struggling with the following task: my data.frame: A1 A2 A3 B1 B2 B3 58.81 53.292 54.501 13.013 17.39 19.407 56.02 56.251 54.033 20.099 13.15 10.411 55.376 53.099 57.625 13.396 21.031 13.22 58.584 53.194 54.218 13.038 16.854 19.289 55.7 55.921 53.847 19.942 13.153 9.828 55.093 52.934

http://biomet.oxfordjournals.org/cgi/reprint/86/3/677 biometrika1999 http://biomet.oxfordjournals.org/cgi/reprint/94/4/1006 biometrika2000 Hi All: I just want to try some luck. I am currenly working on my project,one part of my project is to reanalysis the kenward cattle data by using the method in Mohsen's paper,but I found I really can get the same or close output as he did,so,any

Hello Coming from matlab background, I could use some help on this one. I need to generate a data set based on this equation X(t) = 3.8x(t-1) (1-x(t-1)) + e(t), where e(t) is a N(0,0,001) random variable I need say 100 values. How do I do this? Thanks

Hello, I am a new R user and I am trying to estimate some generalized linear models (glm). I am trying to compare a model with a gaussian distribution and an identity link function, and a poisson model with a log link function. My problem is that while the gaussian model has significantly lower (i.e. "better") AIC (Akaike Information Criterion) most of the coefficients are not

Thanks for the help so far, I have already solved the main problem, but sadly it seems I can't find out the last step. I hope you can help me there... Here ist the important part: It is trying to load some nvidia module, although i have already done apt-get purge nvidia* [ 4.636] (II) LoadModule: "glx" [ 4.638] (II) Loading /usr/lib/xorg/modules/extensions/libglx.so [

Thanks for the help so far, I have already solved the main problem, but sadly it seems I can't find out the last step. I hope you can help me there... Here ist the important part: It is trying to load some nvidia module, although i have already done apt-get purge nvidia* [ 4.636] (II) LoadModule: "glx" [ 4.638] (II) Loading /usr/lib/xorg/modules/extensions/libglx.so [

Hello Folks, I have a database of 2000+ days with 35 observations each. I am trying to modeling a time series by day, but it seems a problem that I don<t have the time of the observation. I achieve something interesting by using the barplot function, but I`d rather working with ggplot2, since I have the book by Hadley Wickham in hands. I start by transforming my data into a dataframe, but

Hello, I've been struggling with this issue for a few days now. My configuration is pretty much a stock Ubuntu 12.04 server, with samba4 and bind9 installed from packages. Samba - 4.0.0alpha18 bind - 9.8.1-P1 I've recompiled the bind package to add dlz support (added flag --with-dlopen), but I'm getting an error whenever a domain client connects and attempts to update its DNS

Hi, Is anybody know the purpose of this method (xen_evtchn_do_upcall)? When I run a user level application involved in TCP receiving and the SoftIRQ for eth0 on the same CPU core, everything is OK. But if I run them on 2 different cores, there will be xen_evtchn_do_upcall() existing (maybe when the local_bh_disable<http://www.cs.fsu.edu/~baker/devices/lxr/http/ident?i=local_bh_disable>() or

Hi, all: Two examples are shown below. I want to use the multiple comparison of Dunnett. It succeeded in upper case "example 1". However, the lower case "example 2" went wrong. In "example 2", the function pmvt return NaN, so I cannot show this simtest result. Is there any solution? (I changed the variable "maxpts" to a large number in front of the

I have the following data frame: > algaedata = year DIV cellsperml 2001 BAC 72.808 2001 CHL 3.273 2002 BAC 14.002 2002 CYA 220.896 2002 UNI 464.699 2003 BAP 0 2003 BAC 1.782 2004 CYA 315.799 2005 UNI 222.532 2005 BAP 0.2 2005 CYA 163.627 2005 BAC 324.949 2006 CHL 1.636 2006 BAC 199.145 2007 CHL 19.635 2007 CYA 134.174 2007 BAC 485.405 2007 CHL 11.454 2007 CYA 104.721 ...which makes a fine

Hello. I have 16 subjects with 1-4 obs per subject. I am using the package "nlme" to fit a simple random effects (variance components model) with 3 parameters: overall mean (fixed effect), between subject variance (random) and within subject variance (random). I need a 3x3 variance-covariance matrix that includes all 3 parameters in order to compute the variance of a specific

Hello R-classicist won't like the following question due to there's no minimal reproducible example and the question was posted already. Anyway I'm not able understand this rather complicated version of "add points to wireframe". what I have is matrix 3x2000 with a dem<-(x,y,z, coordinates) (not provided). with follwing code I want to produce the wireframe and add some

R Help - I'm trying to run a correlation matrix with a covariate of "age" and will at some point will also want to covary other variables concurrently. I'm using the "psych" package and have tried other methods such as writing a loop to extract semi-partial correlations, but it does not seem to be working. How can I accomplish this? library(psych) > set.cor(y =

I am using the package nlme to fit a simple random effects (variance components model) with 3 parameters: overall mean (fixed effect), between subject variance (random) and within subject variance (random). I have 16 subjects with 1-4 obs per subject. I need a 3x3 variance-covariance matrix that includes all 3 parameters in order to compute the variance of a specific linear

Hi, I am trying to find m breakpoints of a linear regression model. I used the segmented package. It works fine for small number of predicators and breakpoints.(3 r.v. 3 points). However, my model has 14 variables it even would not work even for just one breakpoints!. The error message is always estimated breakpoints are out of range. Since my problem is time related problem. So I