similar to: R in MAC add many extra <´s

Dear R users, I am trying to draw error bars in a bar plot, I use this code (tried many others which did not work): install.packages() library(gplots) y <-c(39.02, 46.42) se <- c(7.57,7.35) plot <- barplot(y, beside=TRUE, col=0, ylim=c(0,47), axis.lty=1, main="far-red", xlab="latitude", names.arg=c("56N", "68N")) superpose.eb(plot, y, se,

Hello I have a data set with outlier and it is not normally distributed. I would instead like to use a more robust distribution like t-distribution. My question is if the coefficients of the regression are different from zero, but assuming a t-distribution. Could someone hint me what package to use or.... Thanks in advance Rosario

Hello I have a data set summarized like this: File name= Height Group Ind Age Trait 1 1 1 20 1 1 2 21 1 2 1 22 1 2 2 21 1 3 1 24 1 3 2 45 1 4 1 23 1 4 2 26 2 1 1 45 2 1 2 12 2 2 1 25 2 2 2 26 2 3 1 45 2 3 2 43 2 4 1 23 2 4 2 47 . . . I would like to plot Trait ~ Age but a different plot for each Group. I tried: > plot(Height$Trait ~ Height$Age | Group) But does not work. Any

Hello I am trying to get out of an lm model the fstatistics, however after I run the model I write > names(Model) and the fstatistic does not appear only these. names(Model) [1] "coefficients" "residuals" "effects" "rank" "fitted.values" [6] "assign" "qr" "df.residual"

Hello I have a model with several hundred Y variables, and also several 1000 X variables. The model is linear lm(Y ~ X). My questions are: 1.- how to avoid writing all Xs variables? is list() the right function? 2.- about the multiple Ys with dependence among some of them, how to incorporate that information in the linear model? Thank you Rosario

Hello I seem to find only two types of rotation for the factanal function in R, the Varimax and Promax, but is it possible to run a orthogonal and oblique rotations in R? Thanks in advance Rosario

Hello I am using box.cox() and I get this error message: Warning message: 'box.cox' is deprecated. Use 'bcPower' instead. See help("Deprecated") and help("car-deprecated"). I went to help but I did not understand the explanation, I am still wondering what is really happening. Thanks /R

Hello I am trying to draw a basic black and white map of two European countries. After searching some key words in google and reading many pages I arrived to the conclusion that persp() could be used to draw that map. I have prepared three small example files, which are supposed to be the files required for running that function. xvector is a vector with the longitudes yvector is a vector with

Hello I have a problem on keeping the format when I export a matrix file with the write.table() function. When I import the data volcano from rgl package it looks like this in R: > data[1:5,] [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] 100 100 101 101 101 101 101 100 100 100 101 101 102 102 [2,] 101 101 102 102 102 102 102

Hello I have run an anova analysis for the fallowing model: H_obs=mu+REGION+MANAGEMENT + e When I run it in ASRelm I get the p-value for mu, and, of course also for the two dependent variables (REGION and MANAGEMENT) When I run it in R, I do not get the pvalue for mu. Can some one help me to understand why? and if it is possible to estimate the pvalue for mu in anova() in R? I attach the

Hello I have done resampling from a matrix data with two variables (A,B). After resampling a get a set of lists (see below example for resample size 2). In both lists I have the same two variables (A,B). I want to export this data to an excel file where I have only two columns, one per variable. Thus, stack the data for variable A in a single column and the same for variable B. I have tried to

Hello I am trying to compare two models using anova(), however I get a message error (see below). In the net I only found some information on certain library(car) for which one should use anova with A capital letter (Anova instead of anova), but I could not find car library as it says it does not exist. > Model <- lm(interceptG ~ SW + TSC + FSC + PF + SlopeG + K, data=AllTrait) >

Hi I have 1800 response variables to regress on two factors (latitude and age), what is the script to run all response variables at once instead of writing 1800 models? Thanks R. ________________________________________ L?hett?j?: r-help-bounces at r-project.org [r-help-bounces at r-project.org] k&#228;ytt&#228;j&#228;n Dennis Murphy [djmuser at gmail.com] puolesta L?hetetty: 25.

Dear colleges, I would like to use the lmList function. I have installed the lme4 library but when I try to use the lmList function I get this error message. Error: could not find function "lmList" Here you can see what kind of messages I am getting when installing lme4 package. Does anyone know how to solve this problem? Thanks Rosario > install.packages() --- Please select a

Hello I have a data set with fixed and random effects, therefore I am using the lme function: lm(y ~ xfixed, random=~1|xrandom, data) After this I want to get the F-values for both the fixed and random predictors. I can easily get the F-value and df for the xfixed predictors (anova()), but how to get the F-value for the xrandom predictors? Thanks in advance. /R

Dear colleges I need to get the p values for a table with 15000 entries of t values. Does any of you know how to do it? I can, of course, get one by one but that is not sensible. Thanks Rosario

Dear list, I' m studying FLAC performance, and I'd like to know how much compression can be achieved for different audio files. 1) It seems that for nontonal sound (wideband noise), the compression factor is better than for compound sound (tones + nontonal components), which is typically 2. The reason for this result could be the following: the LPC filter is more suitable for

Jon, I need to create a database of basics blocks and their execution time. The only thing I'm concerned is if a block A is more expensive than a block B. Do you think that even with the overhead I would be able to get the A > B information? Like: overhead + time(A) > overhead + time(B) => A > B. If so, I'm not too much concerned about the accuracy. Not sure if I was clear,

Hi all, I'm installing Samba 3.0.5 in Mac OSX 10.3.4 server as root. First, ./configure ----with-krb5=/usr/sbin Second, when I do 'make' show this error: Compiling libsmb/clikrb5.c libsmb/clikrb5.c: In function `krb5_locate_kdc': libsmb/clikrb5.c:188: error: `krb5_krbhst_handle' undeclared (first use in this function) libsmb/clikrb5.c:188: error: (Each undeclared

Hello I am trying to run a PCA on the attached file, but I get this error message: pc<-prcomp(data[,-(1:2)],scale=T)$x Error in svd(x, nu = 0) : infinite or missing values in 'x' Thanks in advance /R -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: Romania_PCA_Chlor1.txt URL: