[R] What is the difference between "probe level data" and "summarized data"

2012-08-07 Thread ali_protocol
Hi all,

What is the difference between "probe level data" and "summarized data".
Thanks in advance.



--
View this message in context: 
http://r.789695.n4.nabble.com/What-is-the-difference-between-probe-level-data-and-summarized-data-tp4639376.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Weired location of points on MA plot?

2012-05-05 Thread ali_protocol
Dear all,

 I have a MA plot for "replicate" single channnel arrays, the loess curve is
aligned to M=0. However,  I expect that M values be as close as possible to
M=0 for replicate arrays, which is not the case. How may I obtain a measure
that indicates the average distance of points to M=0?

Thanks indeed.

--
View this message in context: 
http://r.789695.n4.nabble.com/Weired-location-of-points-on-MA-plot-tp4610624.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Obtaining translated, rotated and scaled matrics in procrustes analysis

2012-04-26 Thread ali_protocol
Hi all,

I wondered how I can get the translated, rotated and scaledmatrix in a
procrustes analysis.  This does not help:


##
library (vegan)
#defining the target matrix:
a= cbind (c (1,2,3,4), c( 10,12,14,16))
#defining the matrix to be rotated:
b= a; b[,2]= c(9,11,15,16)
c=procrustes (a,b); g=fitted (c)
# trying to get to the procrustes matrices:
d= t (apply(b,1,'-',c$translation))
#points ( d, col=4)
 e= d %*% c$rotation
#points ( e, col=6)
f= e*  c$scale
#points ( f, col=7)
# Should be  a matrix(0,3,2) but is not:
f-g

#(!) Thanks in Advance.


--
View this message in context: 
http://r.789695.n4.nabble.com/Obtaining-translated-rotated-and-scaled-matrics-in-procrustes-analysis-tp4589030p4589030.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Comparing 2 means. pls help

2012-04-09 Thread ali_protocol
I am interested in the difference of 2 data:

mat1= c(2.2, 2.3, 2.2,2.5)
mat2= c(2.6, 2.8, 2.7,2.4)

mat= mat2-mat1

I perform an action on both mat1 and mat2, and I get
mat1prime and mat2prime:

mat1prime= c(2.5, 2.5, 2.3,2.5)
mat2prime= c(2.6, 2.8, 2.7,2.6)

matprime= mat2prime-mat1prime


I want to check wether mat and matprime have the same mean.
How should I do this in R?

--
View this message in context: 
http://r.789695.n4.nabble.com/Comparing-2-means-pls-help-tp4542237p4542237.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] Sum of sd between matrix cols vs spearman correlation between them

2012-04-07 Thread ali_protocol
Hi Petr,

Thank you. Yes, normalized.matrix" and "normalized" are the same, my mistakes.

Operators: like Cyclic loess normalization and quantile Normalization,
from limma or affy packages.

Would cv= sd/mean be a better replacement for sd?

On 4/5/12, Petr Savicky [via R]
 wrote:
>
>
> On Thu, Apr 05, 2012 at 07:46:52AM -0700, ali_protocol wrote:
>> Hi all,
>>
>> I have a matrix  (n*2), I want to compare 2 operators (2 normalization for
>> array results) on these matrix.
>> The 2 columns should ideally become the same after operations
>> (normalization). So to compare operations,
>> I do this for each normalization:
>>
>> s= sum (apply (normalized.matrix, 2,sd))
>> c= cor (normalized[,1],normalized [,2], method='pearson')
>>
>>
>> I expect that if normalization 1 is superior, s should be less and c
>> greater
>> than normalization2, but both s and c change in 1 direstion. Is this
>> possible or am I doing something wrong?
>
> Hi.
>
> Is "normalized.matrix" and "normalized" the same matrix?
>
> Can you specify, which operators you use for normalization?
>
> Without having this information, i guess that comparing the
> correlations alone can be used, since it does not depend on
> the scaling the numbers. By an appropriate scaling factor,
> the sd may be changed to any value, but this does not say
> much about the amount of information in the data. On the other
> hand, the correlation does not change by scaling, so it may
> be a more reliable measure.
>
> Note that
>
>   apply(normalized.matrix, 2, sd)
>
> is the same as
>
>   sd(normalized.matrix)
>
> Hope this helps.
>
> Petr Savicky.
>
> __
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
>
> ___
> If you reply to this email, your message will be added to the discussion
> below:
> http://r.789695.n4.nabble.com/Sum-of-sd-between-matrix-cols-vs-spearman-correlation-between-them-tp4535057p4535407.html
>
> To unsubscribe from Sum of sd between matrix cols vs spearman correlation
> between them, visit
> http://r.789695.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=4535057&code=bW9oYW1tYWRpYW5hbGltb2hhbW1hZGlhbkBnbWFpbC5jb218NDUzNTA1N3wtMTkwNDgxNTk3NQ==


--
View this message in context: 
http://r.789695.n4.nabble.com/Sum-of-sd-between-matrix-cols-vs-spearman-correlation-between-them-tp4535057p4538971.html
Sent from the R help mailing list archive at Nabble.com.
[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Sum of sd between matrix cols vs spearman correlation between them

2012-04-05 Thread ali_protocol
Hi all,

I have a matrix  (n*2), I want to compare 2 operators (2 normalization for
array results) on these matrix. 
The 2 columns should ideally become the same after operations
(normalization). So to compare operations, 
I do this for each normalization: 

s= sum (apply (normalized.matrix, 2,sd))
c= cor (normalized[,1],normalized [,2], method='pearson')


I expect that if normalization 1 is superior, s should be less and c greater
than normalization2, but both s and c change in 1 direstion. Is this
possible or am I doing something wrong?

Thank you in advance.




--
View this message in context: 
http://r.789695.n4.nabble.com/Sum-of-sd-between-matrix-cols-vs-spearman-correlation-between-them-tp4535057p4535057.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] A introductory question about Zips law (Newbie to statistics)

2012-03-31 Thread ali_protocol
Hi everyone.

Newbie to statistics.

I have 40 matrices of ~400 values. how may I determine whether the
distribution follows zips law?

response <-sample (1:20,400*4, replace= TRUE)

Thank you vry much. 

--
View this message in context: 
http://r.789695.n4.nabble.com/A-introductory-question-about-Zips-law-Newbie-to-statistics-tp4521190p4521190.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] How to test for the difference of means in population, please help

2012-03-29 Thread ali_protocol
Dear Wu,

Thank you indeed very  much  for your response. Is it possible that I
use this method to compare the distribution of original population
that  this observations come from?

Thank you ideed very much.

On 3/28/12, alan wu [via R]  wrote:
>
>
> try anova if they have same observe points.
>
> conditions<-rep(c(1,2),each=1600)
> response<-sample (1:20,1600*2, replace= TRUE)
> points<-rep(rep(1:400, each=4),2)
> replicates<-rep(1:4,800)
> obs<-data.frame(response,conditions,points,replicates)
> fit <- lm(response~conditions+points,data=obs)
> anova(fit)
>
>
> On Mar 27, 3:13 am, ali_protocol 
> wrote:
>> Dear all,
>>
>> Novice in statistics.
>>
>> I have 2 experimental conditions. Each condition has ~400 points as its
>> response. Each condition is done in 4 repereats (so I have 2 x 400 x 4
>> points).
>>
>> I want to compare the means of two conditions and test whether they are
>> same
>> or not. Which test should I use?
>>
>> #populations
>> c = matrix (sample (1:20,1600, replace= TRUE), 400 ,4)
>> b = matrix (sample (1:20,1600, replace= TRUE), 400 ,4)
>>
>> #means of repeats
>> c.mean= apply (c,2, mean)
>> b.mean= apply (b,2,mean)
>>
>> #mean of experiment
>> c.mean.all= mean (c)
>> b.mean.all= mean (b)
>>
>> --
>> View this message in
>> context:http://r.789695.n4.nabble.com/How-to-test-for-the-difference-of-means...
>> Sent from the R help mailing list archive at Nabble.com.
>>
>> __
>> r-h...@r-project.org mailing
>> listhttps://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting
>> guidehttp://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>
> __
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
>
> ___
> If you reply to this email, your message will be added to the discussion
> below:
> http://r.789695.n4.nabble.com/How-to-test-for-the-difference-of-means-in-population-please-help-tp4508089p4511310.html
>
> To unsubscribe from How to test for the difference of means in population,
> please help, visit
> http://r.789695.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=4508089&code=bW9oYW1tYWRpYW5hbGltb2hhbW1hZGlhbkBnbWFpbC5jb218NDUwODA4OXwtMTkwNDgxNTk3NQ==


--
View this message in context: 
http://r.789695.n4.nabble.com/How-to-test-for-the-difference-of-means-in-population-please-help-tp4508089p4516385.html
Sent from the R help mailing list archive at Nabble.com.
[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Testing for difference of distribution of two population(newbie)

2012-03-27 Thread ali_protocol
Dear all,

I am novice in statistics.
I have two matrices (results from my experiments), each having ~400 point. I
want to test whether the points come from a same distribution or not.
Further, I have the results for each matrix (each experiment) in 4
replicates. What should I do?

# Repeats of experiments are in different columns
exp.1 = matrix (sample (1:20,1600, replace= TRUE), 400 ,4)
exp.2 = matrix (sample (1:20,1600, replace= TRUE), 400 ,4)


Thanks in advance

--
View this message in context: 
http://r.789695.n4.nabble.com/Testing-for-difference-of-distribution-of-two-population-newbie-tp4508092p4508092.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] How to test for the difference of means in population, please help

2012-03-27 Thread ali_protocol
Dear all,

Novice in statistics.

I have 2 experimental conditions. Each condition has ~400 points as its
response. Each condition is done in 4 repereats (so I have 2 x 400 x 4
points). 

I want to compare the means of two conditions and test whether they are same
or not. Which test should I use?

#populations
c = matrix (sample (1:20,1600, replace= TRUE), 400 ,4)
b = matrix (sample (1:20,1600, replace= TRUE), 400 ,4)

#means of repeats
c.mean= apply (c,2, mean)
b.mean= apply (b,2,mean)

#mean of experiment
c.mean.all= mean (c)
b.mean.all= mean (b)

--
View this message in context: 
http://r.789695.n4.nabble.com/How-to-test-for-the-difference-of-means-in-population-please-help-tp4508089p4508089.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] How to choose the most 2 different samples in an array experiment

2012-03-15 Thread ali_protocol
Dear all,

What test should be done on ~40 samples (matrixes of 400x1), each in
triplicate, to determine which sample is the most different to all samples?

and What test should be done to determine the two samples which are most
different? 

--
View this message in context: 
http://r.789695.n4.nabble.com/How-to-choose-the-most-2-different-samples-in-an-array-experiment-tp4474424p4474424.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Creat a new list with names of other list

2012-03-15 Thread ali_protocol
Hi all,

I want to create a new list with names of another list. Eeach sublist from
both lists is a matrix, but the matrices with the same name (eg.  list.1
[[1]] and list.2 [[1]]) have different dimensions. How can I create the 2nd
list?
 
tSE = list ()
Norm <- names(Normal)
names(tSE) <- Norm 


does not do it.


--
View this message in context: 
http://r.789695.n4.nabble.com/Creat-a-new-list-with-names-of-other-list-tp4474250p4474250.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Apply a loop containing a function on a list

2012-03-14 Thread ali_protocol
Hi all,
I want to do this:

B.list$aa= (a loop containing My.fun acting on the reults of second
function on a A.list$aa))
or, overally
B.list$aa = function (A.list$aa)
B.list and A.list has many sublists aa, ab and Is there a way I can
apply the function and loop on all sublists of A.list and get B.list?


Thanks in advance.

--
View this message in context: 
http://r.789695.n4.nabble.com/Apply-a-loop-containing-a-function-on-a-list-tp4471188p4471188.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Changing the size of a plot

2012-03-07 Thread ali_protocol
Dear all,

Newbie.
How can I make R create a bigger window for my plot? Thank you in advance.

--
View this message in context: 
http://r.789695.n4.nabble.com/Changing-the-size-of-a-plot-tp4452761p4452761.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] What is the fastest way to combine all columns of a matrix to one column?

2012-02-27 Thread ali_protocol
Dear all,

Newbie: What is the fastest way to combine all columns of a matrix to one
column?

--
View this message in context: 
http://r.789695.n4.nabble.com/What-is-the-fastest-way-to-combine-all-columns-of-a-matrix-to-one-column-tp4423882p4423882.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] What is the function for "smoothing splines with the smoothing parameter selected by generalized maximum likelihood?

2012-01-09 Thread ali_protocol
Dear all, 

I am new to R, and I am a biotechnologist, I want to fit a smoothing spline
with smoothing parameter selected by generalized maximum likelihood. I was
wondering what function implement this,  and, if possible how I can find the
fitted results for a certain point (or predict from the fitted spline if
this is the correct  language)

--
View this message in context: 
http://r.789695.n4.nabble.com/What-is-the-function-for-smoothing-splines-with-the-smoothing-parameter-selected-by-generalized-maxi-tp4278275p4278275.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] ssanova/ ssanova0 and adding the fitted line to a plot.

2012-01-05 Thread ali_protocol
Hi everyone,

I use this code to add the fitted line to a plot, however, I think I cannot
access the fitted equation. How my I find  the equation  in the object ksi?
thanks alot.

library (gss)
generate (simple linear) x and y
x= 2* (1:40)
y= x-1
plot (x,y)
## fit and draw
ksi = ssanova0 (x~y, method= "m")
lines( ksi$qwk, col = "red")

--
View this message in context: 
http://r.789695.n4.nabble.com/ssanova-ssanova0-and-adding-the-fitted-line-to-a-plot-tp4264679p4264679.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] How to create a matrix with 3 dimensions from several 2 dimensional matrice?

2011-12-27 Thread ali_protocol
Hi every  one,

How is it possible  to create a matrix with 3 dimensions from several 2
dimensional matrice? 
Is it possible that each of "elementary/building block" matrices could be
called by its corresponding  original name?  

Thanks alot.

--
View this message in context: 
http://r.789695.n4.nabble.com/How-to-create-a-matrix-with-3-dimensions-from-several-2-dimensional-matrice-tp4237360p4237360.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Cumulative distribution of Spearman's Rank rho values.

2011-12-27 Thread ali_protocol
Hi every one,

newbie here.

I have a matrix and I want to compute the cumulative distribution of
Spearman's Rank rho values, for variables (which are in rows)

 m= matrix ( rnorm (30), 10,3)

how is that done? 

--
View this message in context: 
http://r.789695.n4.nabble.com/Cumulative-distribution-of-Spearman-s-Rank-rho-values-tp4237267p4237267.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Cumulative distribution of coefficient of variation.

2011-12-27 Thread ali_protocol
Hi every one,

newbie here.

I have a matrix and I want to compute the cumulative distribution of
coefficient of variation, for variables (which are in rows)

 m= matrix ( rnorm (30), 10,3)

how is that done?


--
View this message in context: 
http://r.789695.n4.nabble.com/Cumulative-distribution-of-coefficient-of-variation-tp4237243p4237243.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] Estimate of smooth.spline

2011-12-24 Thread ali_protocol
thank you vey much Jean.

On 12/21/11, Jean V Adams [via R]
 wrote:
>
>
> ali_protocol wrote on 12/21/2011 05:39:19 AM:
>
>> Hi everyone,
>> I have:
>>
>> s= smooth.spline (cbind(1,2,3,4,3,3),cbind (4,2,4,6,5,6))
>>
>> how may I obtain s hat (s^)?
>>
>> Thanks a lot.
>
>
> Read the help file on the smooth.spline() function,
>
> ?smooth.spline
>
> The section titled "Value" tells you what the function returns ...
>
>An object of class "smooth.spline" with components
>xthe distinct x values in increasing order, see the ?Details?
> above.
>y the fitted values corresponding to x.
>...
>
> The y component is what you're interested in,
>
> s$y
>
>
> Jean
>   [[alternative HTML version deleted]]
>
> __
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
>
> ___
> If you reply to this email, your message will be added to the discussion
> below:
> http://r.789695.n4.nabble.com/Estimate-of-smooth-spline-tp4221253p4221894.html
>
> To unsubscribe from Estimate of smooth.spline, visit
> http://r.789695.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=4221253&code=bW9oYW1tYWRpYW5hbGltb2hhbW1hZGlhbkBnbWFpbC5jb218NDIyMTI1M3wtMTkwNDgxNTk3NQ==


--
View this message in context: 
http://r.789695.n4.nabble.com/Estimate-of-smooth-spline-tp4221253p4230616.html
Sent from the R help mailing list archive at Nabble.com.
[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Estimate of smooth.spline

2011-12-21 Thread ali_protocol
Hi everyone,
I have:

s= smooth.spline (cbind(1,2,3,4,3,3),cbind (4,2,4,6,5,6))

how may I obtain s hat (s^)?

Thanks a lot.

--
View this message in context: 
http://r.789695.n4.nabble.com/Estimate-of-smooth-spline-tp4221253p4221253.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] Bioconductor. MA plot for qPCR array

2011-12-18 Thread ali_protocol
Dear "Michael Weylandt" and "Juliet Hannah". Thank you indeed very much.

--
View this message in context: 
http://r.789695.n4.nabble.com/Bioconductor-MA-plot-for-qPCR-array-tp4182805p4210685.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Smoothing spline with smoothing parameters selected by "generalized maximum likelihood"

2011-12-18 Thread ali_protocol
Hi there,

How may I  smooth spline  two vectors with the smoothing parameter selected
by generalized maximum likelihood (GML) .?

Thanks a lot.

--
View this message in context: 
http://r.789695.n4.nabble.com/Smoothing-spline-with-smoothing-parameters-selected-by-generalized-maximum-likelihood-tp4210679p4210679.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Binning a 2 column matrix by avarages of rows.

2011-12-17 Thread ali_protocol

Newbie here.  Many apologies in advance for using the incorrect lingo. 
I'm new to statistics and VERY new to R.
I have a "nx2" matrix , I want to sort the values based on the average of 2
columns and put k lowest (or highest) values in bin1, second k high/low
values in bin2, and so on (bins would be of the same dimensions). I should
also know what the first index (or position) of the pair has been.


car-1   29  30
car-2   22  24
car-3   16  16
car-4   41  41
car-5   43  45
car-6   34  36

for exmample in bins of 2x2

bin1
car-3   16  16
car-2   22  24
bin2
car-1   29  30
car-6   34  35
bin3
car-4   41  41
car-5   43  45

--
View this message in context: 
http://r.789695.n4.nabble.com/Binning-a-2-column-matrix-by-avarages-of-rows-tp4207606p4207606.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Bioconductor. MA plot for qPCR array

2011-12-11 Thread ali_protocol
Dear all,

Is there anyway too generate MA plot for 2 qPCR assays (an array of 2x 400).


--
View this message in context: 
http://r.789695.n4.nabble.com/Bioconductor-MA-plot-for-qPCR-array-tp4182805p4182805.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] span and epsilon in Cyclic loess

2011-12-08 Thread ali_protocol
Hi everyone,

Is there any automated method to optimize span and epsilon in cyclic loess
normalization? I'm using cyclic loess normalization for an array with <1000
genes. 

Thanks in advance.

--
View this message in context: 
http://r.789695.n4.nabble.com/span-and-epsilon-in-Cyclic-loess-tp4171916p4171916.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] How to it a "loess curve" and obtain the equation in R?

2011-10-21 Thread ali_protocol
Hi!

How can I fit a loess curve to an array (384 x 2).
How can I obtain the equation for thi fi?

Thanks in advance.

--
View this message in context: 
http://r.789695.n4.nabble.com/How-to-it-a-loess-curve-and-obtain-the-equation-in-R-tp3927622p3927622.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] Creating affybatch objects from matrix (data from qPCR array)

2011-10-21 Thread ali_protocol
Thanks Martin.

Actually I want to do a "cyclic loess", so I chose Affy. How may I find
other packages that perform Cyclic loess?

--
View this message in context: 
http://r.789695.n4.nabble.com/Creating-affybatch-objects-from-matrix-data-from-qPCR-array-tp3921559p3927553.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Creating affybatch objects from matrix (data from qPCR array)

2011-10-20 Thread ali_protocol
Hi!

Is There a way to manually create an affybatch object from qPCR array data?

--
View this message in context: 
http://r.789695.n4.nabble.com/Creating-affybatch-objects-from-matrix-data-from-qPCR-array-tp3921559p3921559.html
Sent from the R help mailing list archive at Nabble.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.