similar to: RCMD Warnings on src directory.

Hi all, I am trying to find sd of my rows in a matrix and i get column sd inspite of extracting rows. I tried to do the sqrt(var(x)) but that did'nt work as well, Here is my data genes 15 24 63 40 25 42 46 35 23 53 37 45 30 37 50 55 40 51 30 48 x<-sd(genes[1:5,]) y<-sqrt(var(genes[1:5,])) I get 4 sds for the 4 columns instead of 5 sds for my 5 rows. Thanks you in advance.

Hi, my vector V<- c(1,0.5,0.06,0.056,0.01,0.04,0.4,0.9,0.82,0.1) if( V > 0.5) { V <- 1 - V } I get a warning saying only the first element will be used in comparing (if V > 0.5). However, my results tell me vis-versa ,it actually compares every element of the vector V with 0.5 and that is waht i want it to do. Using a for loop is expensive and time consuming since my actual

Hi all, My data genes [,1] [,2] [,3] [,4] [1,] 25 72 23 55 [2,] 34 53 41 33 [3,] 26 43 26 44 [4,] 36 64 64 22 [5,] 47 72 67 34 stu<-t.test(genes[,1:2],genes[,3:4]) > stu$p.value [1] 0.4198002 i get 1 pvalue for the entire col1:col2 Vs col3:col4. I am trying to get 5 p values for the 5 rows i have. I am trying to avoid a for loop coz my

Dear all, I would like to fit a linear regression with replication (on each year, observation is replicated, e.g 4 times). The independent variable ranges for instance 1-5 year, so I expect to have a linear fit of 5 points. For that purpose I do these (with dummy variables x and y): x<-rep(seq(1:5),4) y<-rnorm(20) linreg<-lm(y~x) fitted.values(linreg) # why produce 20 points of

Hi All, I get the following error when i perform a t.test. studentt<-apply(tlr, 1, function(x) t.test(x[2:6], x[7:11],var.equal=TRUE)$p.value) # tlr is a table of 11 columns and 22500 rows. I am not able to post tlr due to its size. Error in if (stderr < 10 * .Machine$double.eps * max(abs(mx), abs(my))) stop("data are essentially constant") : missing value where

I have a simple table of data: Result Var1 Var2 Var3 1 0.10 0.78 0.12 0.38 2 0.20 0.66 0.39 0.12 3 0.10 0.83 0.09 0.52 4 0.15 0.41 0.63 0.95 5 0.60 0.88 0.91 0.86 6 -0.02 0.14 0.69 0.94 I am trying to achieve two things: 1) Manipulate this data so that I have the "Result" data unchanged, and all the other data

I want to use the grplasso package on a data set where I want to fit a linear model.? My interest is in identifying significant?beta coefficients.? The documentation is a bit cryptic so I'd appreciate some help. ? I know this is a strategy for large numbers of variables but consider a simple case for pedagogical puposes.? Say I have?two 3 category predictors (2 dummies each), a binary

I'm trying to overlay a density plot on a previously plotted histogram. However, i need to use the same axis as of the 1st(histogram) plot to plot the second. My second plot is creating its own axis and causing my plot to extend the entire histogram instead of getting a subplot on a portion of the histogram. I tried 'fig' and 'new' parameters with no luck. Thanks in advance.

Hi all, I think the last time i posted this topic i started on the wrong foot. Thnaks alot to everyone who responded. i'm coding in R(first time) for a paper my colleague is publishing. i plotted a histogram for 6000 values. I am told to plot experimental vs theoretical vlaues from the histogram and do a non linear least square curve minimization and compute the mean and sd of the new x

Trying to do a meta-analysis of correlations in R using the meta package; have tried several things and keep getting a similar error. Can anyone help explain the error? > cor<-c(-0.3018, 0.667, -3.8002, -0.607, -0.4885, -3.8002, -0.0701, 0.1348, -0.9505, -0.5709, -0.6127, -1.2419, -0.1511, -0.1054)> n<-c(3,4,3,3,3,3,16,36,30,9,3,3,30,4)> library(meta)> metacor(cor, n, data=NULL,

Hi all, I'm trying to use the nls function to do a nls curve fit. I have a histogram for my data. I need to get athe best curve fit for this. i am plotting experimental vs the theoretical values. my experimental values are got from my dnorm function, mean=mean of my mid values(hist) and sd=sd of my mid values. For the theoretical i'm using mean and sd of all 20000 raw values as the start

Hi all, I'm trying to customize the nls function code to spit out the sd and mean of the new fit. All it gives me are the fitted values. THis function seems to work fine for my project otherwise. All i need is the SD and MEAN. Can someone help me to customise. I frankly did not understand the nls code at all. I have 10 x and 10 y values from a histogram. I need to use nls to do the best fit

How can I receive parameter estimates for a given curve that has been fit with the NLS function. The problem is that I have a "n parameter" curve and I want the optimal fit. The NLS procedure gives me final function values and not the individual parameter estimates that were used to define this "best" fit. What function can I use to get these parameters? Thanks

Good Day All ! I have a 2-D vector of mode numeric and a parallel 1-D vector of mode numeric. Here are my values. 2-D vector 1-D vector 80 75 85 80 80 70 80 90 80 80 60 70 80 70 70 85 75 95 85 85 70 60 90 60 70 My 1-D vector is the average across the rows of my 2-D vector. I process each column of the 2-D vector in a loop. When I trim my 2-D vector column

Hi all, I downloaded all the required r-tools from Dr.Brian Ripley's website. RCMD build dAnal It built the tar file for me. After which i tried RCMD check dAnal and i get the following error. Making package dAnal adding build stamp to DESCRIPTION file installing R files installing data files installing man source files installing indices Error in

Hi, I wonder if it is possible to create an DLL or EXE file performing R procedures. Instead of running R, reading data and calling some procedures, I would like to use R functions in the following way: "C:\linearRegression.exe data.txt" which would produce let's say file "output.txt" with the results. Is there some way how to do it? Thanks a lot. Pavel Vanecek

i'm coding in R(first time) for a paper my colleague is publishing. i plotted a histogram for 6000 values. Taking the mean and sd of the midpoints i did a dnorm and got the densities. pl<-dnorm(trimmedvals,mean=midsmean,sd=midsSD) (now in a loop of 5 times) i plotted these experimental values against theoretical values using the nls function the following way.

Doesn't look like nonparametric fit to me, since nls() is used to fit to a gaussian density, so the result is a gaussian density (with estimated parameter). What I do not understand is why people would do this. This is not the first time I've seen people doing this, on both R-help and S-news (if my memory is still any good). If the objective is to fit a Gaussian distribution to the

Hi All, Running R CMD javareconf -e (or R CMD javareconf as root I am rot on my machines) Fails. / [root at penguins1lanalt etc]# R CMD javareconf Java interpreter : /usr/bin/java Java version : 1.6.0_30 Java home path : /usr/java/jre1.6.0_30 Java compiler : not present Java headers gen.: Java archive tool: Java library path:

Building a package using R.2.8.0 on Windows XP gives problems when compiling C code: Trying to run 'rcmd build' on R.2.8.0 on Windows XP gives ... ---------- Making package gRbase ------------ adding build stamp to DESCRIPTION making DLL ... gcc -std=gnu99 -Ic:/Programs/R/current/include ------------------------------------ gcc.exe: no input files make[3]: *** [setops.o] Error 1